Aggregation and complex formation of amyloid beta (Aβ) peptides on a neuronal cell membrane is a hallmark of neuro-disturbance diseases. In this work, we performed molecular dynamics (MD) simulations to investigate the initial stage of interactions of multiple Aβ₄₂ peptides on a GM1 ganglioside-containing membrane that mimics a micro-domain on the neuronal cell surface. Conformational changes of Aβs due to adhesion on the membrane and subsequent molecular interactions among the Aβs were monitored. It was suggested from results of the two 1.0 µs simulation trials that stable complexes of Aβ peptides were not rapidly generated but that a steady binding of two Aβs was gradually formed. Observation of two Aβs that will be a complex with steady binding revealed that one Aβ was bound to the membrane surface, while the other was attached to the first one without strong contact with the membrane. The motion of the first one was restricted and its conformational change was limited, with the basic side-chains of Arg5 and Lys28 working as anchors to hold the Aβ helix region on the membrane. In contrast, the second one had high flexibility and showed diversity in its conformation. The second Aβ can search for an energetically favorable binding position on the first one. A parallel β-sheet structure was formed between the C-terminal sides of the two Aβs. Ala30 was critically important to lead the stable β-sheet conformation at the C-terminal hydrophobic domains of Aβs. In the N-terminal sides, helix structures were kept in both Aβs.