Next Article in Journal
Melatonin and Structurally-Related Compounds Protect Synaptosomal Membranes from Free Radical Damage
Next Article in Special Issue
Lipid Vesicle Aggregation Induced by Cooling
Previous Article in Journal
Molecular Mechanisms of Microcystin Toxicity in Animal Cells
Previous Article in Special Issue
Characterization and Purification of Polydisperse Reconstituted Lipoproteins and Nanolipoprotein Particles
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Self-Assembly of Diamondoid Molecules and Derivatives (MD Simulations and DFT Calculations)

1
Department of Physics, University of Illinois at Chicago, Chicago, IL 60607-7052, USA
2
Departments of BioEngineering, Chemical Engineering and Physics, University of Illinois at Chicago, Chicago, IL 60607-7052, USA
*
Author to whom correspondence should be addressed.
Int. J. Mol. Sci. 2010, 11(1), 288-303; https://0-doi-org.brum.beds.ac.uk/10.3390/ijms11010288
Submission received: 7 December 2009 / Revised: 10 January 2010 / Accepted: 17 January 2010 / Published: 21 January 2010
(This article belongs to the Special Issue Molecular Self-Assembly)

Abstract

:
We report self-assembly and phase transition behavior of lower diamondoid molecules and their primary derivatives using molecular dynamics (MD) simulation and density functional theory (DFT) calculations. Two lower diamondoids (adamantane and diamantane), three adamantane derivatives (amantadine, memantine and rimantadine) and two artificial molecules (ADM•Na and DIM•Na) are studied separately in 125-molecule simulation systems. We performed DFT calculations to optimize their molecular geometries and obtained atomic electronic charges for the corresponding MD simulation, by which we predicted self-assembly structures and simulation trajectories for the seven different diamondoids and derivatives. Our radial distribution function and structure factor studies showed clear phase transitions and self-assemblies for the seven diamondoids and derivatives.

1. Introduction

Diamondoid molecules (which are cage-hydrocarbons) and their derivatives have been recognized as molecular building blocks in nanotechnology [110]. They have been drawn more and more researchers’ attentions to their highly symmetrical and strain free structures, controllable nanostructure characteristics, non-toxicity and their applications in producing variety of nanostructure shapes, in molecular manufacturing, in nanotechnology and in MEMS [6,8]. It is important and necessary to study self-assembly of these molecules in order to obtain reference data, such as temperature, pressure, structure factor, bonding properties, etc. for application in nanotechnology e.g., building molecular electronic devices.
Two lower diamondoids (adamantane and diamantane), three adamantane derivatives (amantadine, memantine and rimantadine) and two artificial molecules (substituting one hydrogen ion in adamantane and diamantine with one sodium ion: ADM•Na and DIM•Na) are studied in this report. We classified them into three groups, as shown in Table 1.
Group 1: Adamantane and Diamantane, the lowest two diamondoids. Due to their six or more linking groups, have found major applications as templates and as molecular building blocks in nanotechnology, polymer synthesis, drug delivery, drug targeting, DNA-directed assembly, DNA-amino acid nanostructure formation, and host-guest chemistry [110]. However these diamondoids do not have good electronic properties which are necessary for building molecular electronics, but some of their derivatives do.
Group 2: Amantadine, Memantine and Rimantadine, the three derivatives of adamantane, have medical applications as antiviral agents and due to their amino groups, they could be treated as molecular semiconductors [8,9].
Group 3: ADM•Na, DIM•Na, the two artificial molecules, substituting one hydrogen ion in adamantane and diamantane with a sodium ion could have potential applications in NEMS and MEMS [810]. This report is aimed at studying the self-assembly and phase transition properties of these seven diamondoids and derivatives for further study of their structures and possibly building molecular electronic devices with them.
We first performed density functional theory (DFT) calculations to optimize initial geometry structures of these seven diamondoids and derivatives and obtained their atomic electronic charges [1116]. Then we performed molecular dynamics (MD) simulation in the study of their self-assembly behaviors. MD simulation methods have been broadly used in studying dynamics of molecules, in spite of its classic approximation [1721]. As typical plastic crystal, or say molecular crystal, the phase transition behaviors of adamantane molecules have been studied by other researchers using MD simulation method [22,23]. Those studies mainly have focused on the transition of the adamantane from FCC to BCC crystal structure or vice versa, and both of these are in the solid phase. Self-assembly, however, is a transition process from a rather random condition of molecules to an ordered state [1]. In the present case of our interest, this process is similar to condensation transition from gas state to liquid state, and then freezing transition from liquid state to solid state, which is a border range of phase transitions. For this reason we chose the Optimized Potentials for Liquid Simulations All Atom (OPLS-AA) force field in our MD simulations. OPLS-AA force field has shown good results not just in the gas and liquid states but also in crystalline phases in a large range of temperatures [17]; moreover, OPLS-AA is a reliable force field for carbohydrates [18] and as a result, obviously, for hydrocarbons. In the OPLS-AA force field, the non-bonded interactions are represented by the general Lennard-Jones plus Columbic form as shown by Equation 1, where φij is the interaction energy of two sites for both intermolecular and intramolecular non-bonded cases:
φ i j = i , j { 4 ε i j [ ( σ i j r i j ) 12 ( σ i j r i j ) 6 ] + q i q j e 2 r i j } k i j
The combining rules used along with this equation are σij = (σii σjj)1/2 and ɛij = (ɛii ɛjj)1/2. For intermolecular interactions (when i, j corresponds to different molecules) kij = 1.0. While for intramolecular non-bonded interactions between all pairs of atoms (i < j) separated by three or more bonds, the 1, 4-interactions are scaled down by kij = 0.5.
Our MD simulations are mainly performed with the use of the Gromacs package [1921] and in order to perform simulations involving small molecules, we added to Gromacs the necessary ITFs “included topology files” for the seven different diamondoids and derivatives we studied. Moreover, Gromacs package is capable enough to accept and apply OPLS-AA force field for MD simulations as well as the fact that it is a fast package with flexible characteristics.
We report a brief summary of our MD simulation strategy: vacuum simulation strategy was used in the first step to determine how many molecules were suitable for our simulation and to obtain equilibrated stable structures. Then we built periodic boundary simulation box; and we used the fast Particle-Mesh Ewald (PME) electrostatics method for the preparation of gas state while we chose the cut-off technique for electrostatics in the self-assembly/freezing procedure.
To simulate the self-assembly (i.e., freezing transition), we applied the simulation annealing technique to gradually decease the system temperature. V-rescale temperature coupling was applied throughout and in all the simulations. This is a temperature-coupling method using velocity rescaling with a stochastic term [18]. The van der Waals force cut-off and neighbor-searching list distance were, both, set to 2.0 nm, which was large enough since above 1.4 nm the van der Waals force was negligible.
Furthermore, in order to obtain the geometries and atomic charges for the MD simulations of all the seven different diamondoids and derivatives we applied the DFT calculation which is a widely used method in quantum chemistry and electronic structure theory.

2. Density Functional Theory (DFT) Calculations

By applying the DFT calculations through Gaussian 03 package [15], we obtained the optimized initial structures of all the seven diamondoids and derivatives and also their atomic electronic charges for the MD simulations. The B3YLP exchange-correlation functional [12] method was chosen with cc-pVDZ basis [14], for all the molecules except for ADM•Na and DIM•Na for which 6-311 + G(d, p) basis [13] was used since sodium was not included in cc-pVDZ basis; NBO (Natural Bond Orbital) analysis [16] was added to calculate the atomic electronic charges which were used as reference to set the atomic charges of nitrogen and sodium atoms, i.e., −0.76 and 0.65, respectively (in electronic charge units: Coulomb). The atomic charges of carbon and hydrogen atoms were mainly opted from default OPLS-AA force field.

3. Molecular Dynamic (MD) Simulation Procedures and Results

The first step was to determine the number of molecules in the simulation box which was suitable to resemble an NVT ensemble. We performed short (20–40 picoseconds) MD simulations of different numbers of adamantane molecules, (i.e., 8, 27, 64, 125, 216, 343, 512, 729 molecules) in vacuum with the temperature set at 100 K, and we successfully obtained all the stable structures as shown in Figure 1.
Although these structures are stable, they may not be used for structural analysis, since the simulations are rather too short. We concluded that a 125-molecule system was big enough to represent the self-assembly behavior of larger systems, which also could save us calculation cost. Therefore, we chose 125 molecules in what is reported in the rest of this paper.
From this point on, all the MD simulation systems involved 125 molecules for the seven different diamondoids and derivatives. In order to perform self-assembly simulations from gas state to liquid state and eventually solid state, we developed the following procedure for the seven diamondoids and derivatives:
  • A MD simulation box with 5 × 5 × 5 = 125 molecules was built by the intrinsic tools in the Gromacs software, as shown in Figure 2a.
  • Then we performed a short MD simulation by letting all the molecules relax in a vacuum at 100 K (as in Figure 1) to obtain the stable self-assembled structure, as shown in Figure 2b.
  • We boxed those stable structures and set the distances from molecules to the boundaries of the simulation boxes to about +3 nm. In this manner we could build simulation systems with reasonable densities which were in the range of their gas states and also made the systems isolated from the adjacent simulation boxes under the periodic boundary conditions, as shown in Figure 2c.
  • Longer (more than 1,000 picoseconds) equilibrating simulations were performed in the NVT ensembles, and we applied the PME method with the cut-off at 2.0 nm and high temperatures (in the range of 500–700 K) in order to make sure the entire system went to the gaseous state, i.e., all the molecules separated from each other and distributed randomly in the simulation box as shown in Figure 2d.
    By applying the above four-step procedure we could equilibrate the system in a gas state. As a result, we eliminated the problem due to directly equilibrating the initially prepared molecular system, Figure 2a, at high temperatures. The equilibrated and stable-structure gaseous state, Figure 2d, was then ready for cooling down towards the self-assembly.
  • The next step of the simulation was the cooling down the equilibrated and stable-structure gaseous state, Figure 2d, from the gas to the liquid state as shown in Figure 2e.
  • Further cooling down of the system resulted in the complete self-assembly of all the molecules (to solid state) as shown in Figure 2f.
The snapshots of the gas-liquid-solid MD simulations for 125 molecules (stages d, e and f of Figure 2) of each of the seven diamondoids and derivatives (Adamantane, Diamantane, Amantadine, Rimantadine, Memantine, ADM•Na, DIM•Na) are reported in Figure 3.
From these snapshots we can directly observe clear phase transitions for each kind of molecules, from the gaseous state to the liquid state, and then their aggregation into a highly condensed (self-assembled) state.
In order to find the equilibrium configuration of the collection of the 125 molecules at every given temperature we used the simulation annealing procedure [1,19]. Every change of 1 K occurred within 10 picoseconds (5,000 time-steps and each time-step was 0.002 ps). With these settings we could observe the self-assembly behavior in the cooling step.
It should be mentioned that we applied the cut-off of 2.0 nm for the electrostatics instead of PME, which was used in step d (Figure 2), in order to: (i). make sure all the interactions among molecules are considered, even beyond 1.4 nm, which is the custom cut-off setting, since after 1.4 nm there is no significant van der waals force; (ii). thus save the computation time at the mean time, since PME requires more computation than the cut-off method.
According to Figure 3 adamantane, diamantane, ADM•Na and DIM•Na form ordered (crystalline) condensed states. However, amantadine, rimantadine and memantine, while they self-assemble, they do not seem to form clear ordered (crystalline) condensed states. We further produced the hydrogen bonds locations of the same self-assembled snapshots of amantadine, memantine and rimantadine as shown in Figure 4.
According to Figure 4 we do not observe any ordered format for the location of hydrogen-bonds, which is an indication of the non-crystalline self-assembled states of amantadine, rimantadine and memantine. Obviously the existence of partial hydrogen-bonds between these molecules is the reason for the lack of a clear crystalline state in their self-assemblies.
We also produced the radial distribution functions (RDFs) and structure factors (SFs) of the seven different diamondoids and derivatives as presented and discussed below. Study of the RDFs and SFs also reveal further about these phase transition features as discussed below. The RDFs we studied and report here are for the centers of the geometry/mass of molecules.

3.1. Adamantane and ADM•Na

The main feature of adamantane is that at low temperatures we can observe ordered crystal structures (the top-right image in Figures 3), which matches the previous experimental and theoretical studies [22]. From the RDF and SF graphs of adamantane at different temperatures (Figure 5) we can observe gas, liquid and solid characteristics.
In the RDF figures of adamantane, higher and sharper peaks can be observed as the temperature decreases and for the SF graphs, the intensity increases as phase transits from gas to liquid and to solid state. We observed similar features in the RDF and SF figures of of ADM•Na as shown in Figure 6.
The phase transition temperatures for ADM•Na are higher than those of adamantane. However, this is attributed to the presence of -Na ion affecting the molecular interactions among adamantanes, i.e., -Na ion in the ADM•Na structure causes stronger bonding than that of the respective -H ion in adamantane. We may also conclude that while the higher phase transition temperatures in ADM•Na, compared to those of adamantane, are due to –Na ion, the geometric structure is determined by the structure of adamantane.

3.2. Diamantane and DIM•Na

Diamantane molecules can self-assemble to a certain type of solid state structure which can be observed from simulation snapshots (Figure 3), however, its self-assembled structure is not as neat as the crystal structure of adamantane. In Figure 7 we report the RDF and structure factor of diamantane in various phases.
From the RDF and structure factor of diamantane, Figure 7, we may also conclude the following: At low temperatures diamantane does not have features of neat crystal structure. DIM•Na shows similar relationship to diamantane as ADM•Na to adamantane, i.e., higher phase transition temperatures due to sodium, while similar self-assembled crystal structures as diamantane molecules (See DIM•Na snapshots in Figure 3 and DIM•Na RDFs in and SFs in Figure 8).

3.3. Amantadine, Rimantadine and Memantine

These three adamantane derivatives self-assemble at higher temperatures compared with adamantane, but they do not seem to have well-organized self-assembled structures as shown in their snapshots in Figure 4. The reasoning for the special self-assembled structures of amantadine, rimantadine and memantine is as follows:
  • . Since there are nitrogen ions in the structure of amantadine, rimantadine and memantine, which makes their attractive intermolecular forces much larger than that of adamantane, higher temperatures should be applied to these systems in order to obtain initial gas state structures in the step d of Figure 2.
  • . During the cooling down process of amantadine, rimantadine and memantine, hydrogen-bonds are formed which, as expected, increases their phase transition temperatures compared to that of adamantane.
As it is observed from their self-assembled structures (snapshots in Figure 4) amantadine, rimantadine and memantine form certain types of self-assembled structure. However, those structures do not seem to be ordered as was in the case of adamantane, i.e., they have no apparent ordered crystalline structure. There may be two main reasons for this: (a). the -NH2 and -CH3 groups present in these molecules break the geometrical symmetry of adamantane; (b). we observed that the hydrogen bonds, due to -NH2 groups are randomly distributed in the bulk structures which makes the entire structures far from an ordered one. In Figures 811 we report the RDFs and SFs of amantadine, rimantadine and memantine for their gas, liquid and self-assembled states.
According to Figures 911 the self-assembled RDFs and SFs analysis of amantadine, rimantadine and memantine show that for most part there is less obvious ordered self-assembled features, i.e., due to lack of sharp peaks as was the case in adamantane. The liquid-state RDF graphs for amantadine, rimantadine and memantine also, for most part, have less obvious liquid-state features either.
We also performed calculations of the hydrogen-bonds saturation as a function of temperature (as shown in Figure 12) as well as hydrogen-bonds distance distribution (as shown in Figure 13-left) and hydrogen-bonds angle distribution (as shown in Figure 13-right) for amantadine, rimantadine and memantine.
According to Figure 12 in the process of phase transition and self-assembly of amantadine, rimantadine and memantine the number of their hydrogen-bonds will reach saturation limit at low temperatures.
From the hydrogen-bond distance-distribution graphs of amantadine, rimantadine and memantine at 60 K (Figure 13-left), we observe several peaks with their highest peaks at about 0.3 nm.
From the hydrogen-bond angle distributions graphs of amantadine, rimantadine and memantine at 60 K (Figure 13-right), we observe that the hydrogen-bond angles are quite randomly distributed for all the three molecules. We already demonstrated the hydrogen-bonds locations and orientations of this group at 50 K (Figure 4) which also showed no orderly features either. Both Figures 4 and 13 indicate that these structures are not orderly self-assembled structures contrary to the adamantane orderly self-assembled state.
In conclusion these three hydrogen-bonded adamantane derivatives just self-assemble at particular temperatures, but may not form well-organized self-assembled crystal structures. The reason is that the geometry structures of these molecules are not symmetric as was the case for adamantane; and their -NH2 and -CH3 segments, make them unable to pack orderly as adamantane does.
In Figures 14 and 15 we report the self-assembled (solid-state) radial distribution functions and structure factors, respectively, of all the seven diamondoids and derivatives with uniform coordinates scales for the purpose of their collective comparison.
According to Figures 14 and 15 we can observe that adamantane and ADM•Na have more sharp peaks than the other five, which is the indication of their orderly self-assembled structures. While other five molecules also show self-assembled characteristic peaks they have less number of such sharp peaks than adamantane and ADM•Na. These results match the image observation of simulations (Figure 3), i.e., adamantane and ADM•Na have more orderly structures in their self-assembled states due to their molecular symmetry, which thus proves that self-assembly of those molecules are structure-dependent.

4. Conclusion

We have performed a detailed molecular dynamics study of the self-assembly process of seven different diamondoids and derivatives due to temperature variations. From the MD simulation study of the seven diamondoids and derivatives, we may conclude the following: (1) The nature of self-assembly in these molecules is a structure-dependent phenomenon. (2) Final self-assembly structures depend on the different bonding types present in the molecular and intramolecular interactions of these various molecules. (3) The artificial molecules (ADM•Na., DIM•Na) still hold neat crystal structures. Although -Na ion increases the phase transition temperature, as does -NH2 in amantadine, rimantadine and memantine. (4) To a large extent the structural features of diamondoids are retained in ADM•Na and DIM•Na. The reasons for the latter might be that: A. The -Na ion has less topology effect than does the -NH2. B. There is no hydrogen-bonding in the structures of ADM•Na and DIM•Na; therefore they can aggregate to ordered self-assembled structures. This feature is very promising, since it allows us to build orderly-shaped self-assemblies suitable for NEMS and MEMS.

Acknowledgments

We are grateful to the computation facility of the University of Illinois at Chicago. All the MD simulation results are visualized by the VMD package [24].

References

  1. Mansoori, GA. Principles of Nanotechnology—Molecular-Based Study of Condensed Matter in Small Systems; World Scientific Pub. Co: Hackensack, NJ, USA, 2005. [Google Scholar]
  2. Mansoori, GA. Diamondoid molecules. Adv. Chem. Phys 2007, 136, 207–258. [Google Scholar]
  3. Zhang, GP; George, TF; Assoufid, L; Mansoori, GA. First-principles simulation of the interaction between adamantane and an atomic-force microscope tip. Phys Rev B 2007, 75. [Google Scholar]
  4. Ramezani, H; Mansoori, GA; Saberi, MR. Diamondoids-DNA nanoarchitecture: From nanomodules design to self-assembly. J. Comput. Theory Nanosci 2007, 4, 96–10. [Google Scholar]
  5. Mansoori, GA; George, TF; Zhang, GP; Assoufid, L. Molecular Building Blocks for Nanotechnology: From Diamondoids to Nanoscale Materials and Applications (Topics in Applied Physics); Springer: New York, NY, USA, 2007; Volume 109, pp. 44–71. [Google Scholar]
  6. Ramezani, H; Mansoori, GA. Molecular Building Blocks for Nanotechnology: Diamondoids as Molecular Building Blocks for Nanotechnology (Topics in Applied Physics); Springer: New York, NY, USA, 2007; Volume 109, pp. 44–71. [Google Scholar]
  7. Ramezani, H; Saberi, MR; Mansoori, GA. Diamondoids and DNA nanotechnologies. Int. J. Nanosci. Nanotechnol 2007, 3, 21–35. [Google Scholar]
  8. Xue, Y; Mansoori, GA. Quantum Conductance and Electronic Properties of Lower Diamondoids and their Derivatives. Int. J. Nanosci 2008, 7, 63–72. [Google Scholar]
  9. Marsusi, F; Mirabbaszadeh, K; Mansoori, GA. Opto-electronic properties of adamantane and hydrogen-terminated sila- and germa-adamantane: A comparative study. Phys. E 2009, 41, 1151–1156. [Google Scholar]
  10. Mansoori, GA; George, TF; Zhang, GP; Assoufid, L. Structure and opto-electronic behavior of diamondoids, with applications as mems and at the nanoscale level. In Progress Nanotechnology Research; Springer: New York, NY, USA, 2009; ; Chapter 1, pp. 1–19. [Google Scholar]
  11. Kohn, W; Sham, L. Self-consistent equations including exchange and correlation effects. Phys. Rev 1965, 140, A1133–A1138. [Google Scholar]
  12. Lee, C; Yang, W; Parr, RG. Development of the Colle-Salvetti correlation-energy formula into a functional of the electron density. Phys. Rev. B 1988, B37, 785–789. [Google Scholar]
  13. Woon, DE; Dunning, TH, Jr. Gaussian basis sets for use in correlated molecular calculations. III. The atoms aluminum through argon. J. Chem. Phys 1993, 98, 1358–1371. [Google Scholar]
  14. Curtiss, LA; McGrath, MP; Blaudeau, J-P; Davis, NE; Binning, RC, Jr; Radom, L. Extension of Gaussian-2 theory to molecules containing third-row atoms Ga-Kr. J. Chem. Phys 1995, 103, 6104–6113. [Google Scholar]
  15. Frisch, MJ; Trucks, GW; Schlegel, HB; Scuseria, GE; Robb, MA; Cheeseman, JR; Montgomery, JA, Jr; Vreven, T; Kudin, KN; Burant, M; Gill, PMW; Johnson, B; Chen, W; Wong, MW; Gonzalez, C; Pople, JA. Gaussian 03, Revision D.01; Gaussian Inc.: Wallingford, CT, USA, 2004. [Google Scholar]
  16. Carpenter, JE; Weinhold, FJ. Analysis of the geometry of the hydroxyl radical by the ldquo. Mol. Struct. (Theochem.) 1988, 169, 41–62. [Google Scholar]
  17. Jorgensen, WL; Maxwell, DS; Tirado-Rives, J. Development and testing of the opls all-atom force field on conformational energetics and properties of organic liquids. J. Am. Chem. Soc 1996, 118, 11225–11236. [Google Scholar]
  18. Damm, W; Frontera, A; Tirado-Rives, J; Jorgensen, WL. OPLS all-atom force field for carbohydrates. J. Comput. Chem 1997, 18, 1955–1970. [Google Scholar]
  19. Lindahl, E; Hess, B; van der Spoel, D. GROMACS 3.0: A package for molecular simulation and trajectory analysis. J. Mol. Mod 2001, 7, 306–317. [Google Scholar]
  20. van der Spoel, D; Lindahl, E; Hess, B; Groenhof, G; Mark, AE; Berendsen, HJC. GROMACS: Fast, flexible and free. J. Comp. Chem 2005, 26, 1701–1718. [Google Scholar]
  21. Hess, B; Kutzner, C; van der Spoel, D; Lindahl, E. GROMACS 4: Algorithms for highly efficient, load-balanced, and scalable molecular simulation. J. Chem. Theory Comput 2008, 4, 435–447. [Google Scholar]
  22. Greig, DW; Pawley, GS. Molecular dynamics simulations of the order-disorder phase transition in adamantane. Mol. Phys 1996, 89, 677–689. [Google Scholar]
  23. Ciccotti, G; Ferrario, M; Memeo, E; Meyer, M. Structural transition on cooling of plastic adamantane: A molecular-dynamics study. Phys. Rev. Lett 1987, 59, 2574–2577. [Google Scholar]
  24. Humphrey, W; Dalke, A; Schulten, K. VMD-visual molecular dynamics. J. Mol. Graph 1996, 14, 33–38. [Google Scholar]
Figure 1. Snapshots from MD simulations using various numbers of adamantine molecules (23, 33, 43, 53, 63, 73, 83, 93) in vacuum with the temperature set at 100 K.
Figure 1. Snapshots from MD simulations using various numbers of adamantine molecules (23, 33, 43, 53, 63, 73, 83, 93) in vacuum with the temperature set at 100 K.
Ijms 11 00288f1
Figure 2. Stages of simulation procedure: (a) Initial 5 × 5 × 5 MD simulation box; (b) Molecules relaxed in vacuum at 100 K after a short MD simulation; (c) Boxed simulation system with overall gas density; (d) Gas phase as a result of equilibrating simulation in the NVT ensemble; (e) The liquid state; (f) Final self-assembly of molecules to solid state.
Figure 2. Stages of simulation procedure: (a) Initial 5 × 5 × 5 MD simulation box; (b) Molecules relaxed in vacuum at 100 K after a short MD simulation; (c) Boxed simulation system with overall gas density; (d) Gas phase as a result of equilibrating simulation in the NVT ensemble; (e) The liquid state; (f) Final self-assembly of molecules to solid state.
Ijms 11 00288f2
Figure 3. Self-assembly snapshots (left to right) of 125 molecules of the seven diamondoids and derivatives (from top: Adamantane, Diamantane, Amantadine, Rimantadine, Memantine, ADM•Na., DIM•Na) as the temperature [in K shown on every snapshot] is decreased.
Figure 3. Self-assembly snapshots (left to right) of 125 molecules of the seven diamondoids and derivatives (from top: Adamantane, Diamantane, Amantadine, Rimantadine, Memantine, ADM•Na., DIM•Na) as the temperature [in K shown on every snapshot] is decreased.
Ijms 11 00288f3aIjms 11 00288f3b
Figure 4. MD snapshots of (from top to bottom) Amantadine, Rimantadine and Memantine at 50 K and their hydrogen bonds locations and orientations.
Figure 4. MD snapshots of (from top to bottom) Amantadine, Rimantadine and Memantine at 50 K and their hydrogen bonds locations and orientations.
Ijms 11 00288f4
Figure 5. Radial distribution functions (left) and structure factors (right) of adamantane at 450 K, 300 K and 150 K for its gas, liquid and solid states, respectively.
Figure 5. Radial distribution functions (left) and structure factors (right) of adamantane at 450 K, 300 K and 150 K for its gas, liquid and solid states, respectively.
Ijms 11 00288f5
Figure 6. Radial distribution functions (left) and structure factors (right) of ADM•Na at 600 K, 450 K and 200 K for its gas, liquid and solid states, respectively.
Figure 6. Radial distribution functions (left) and structure factors (right) of ADM•Na at 600 K, 450 K and 200 K for its gas, liquid and solid states, respectively.
Ijms 11 00288f6
Figure 7. Radial distribution functions (left) and structure factors (right) of diamantane at 500 K, 400 K and 300 K for its gas, liquid and solid states, respectively.
Figure 7. Radial distribution functions (left) and structure factors (right) of diamantane at 500 K, 400 K and 300 K for its gas, liquid and solid states, respectively.
Ijms 11 00288f7
Figure 8. Radial distribution functions (left) and structure factors (right) of DIM•Na at 600 K, 500 K and 400 K for its gas, liquid and solid states, respectively.
Figure 8. Radial distribution functions (left) and structure factors (right) of DIM•Na at 600 K, 500 K and 400 K for its gas, liquid and solid states, respectively.
Ijms 11 00288f8
Figure 9. Radial distribution functions (left) and structure factors (right) of amantadine at 450 K, 370 K and 300 K for its gas, liquid and solid states, respectively.
Figure 9. Radial distribution functions (left) and structure factors (right) of amantadine at 450 K, 370 K and 300 K for its gas, liquid and solid states, respectively.
Ijms 11 00288f9
Figure 10. Radial distribution functions (left) and structure factors (right) of rimantadine at 500 K, 430 K and 250 K for its gas, liquid and solid states, respectively.
Figure 10. Radial distribution functions (left) and structure factors (right) of rimantadine at 500 K, 430 K and 250 K for its gas, liquid and solid states, respectively.
Ijms 11 00288f10
Figure 11. Radial distribution functions (left) and structure factors (right) of Memantine at 500 K, 430 K and 250 K for its gas, liquid and solid states, respectively.
Figure 11. Radial distribution functions (left) and structure factors (right) of Memantine at 500 K, 430 K and 250 K for its gas, liquid and solid states, respectively.
Ijms 11 00288f11
Figure 12. The number of hydrogen bonds for amantadine, memantine and rimantadine vs. temperature. As temperature decreases the numbers of hydrogen-bonds increase, and tend to maximum numbers for the three derivatives.
Figure 12. The number of hydrogen bonds for amantadine, memantine and rimantadine vs. temperature. As temperature decreases the numbers of hydrogen-bonds increase, and tend to maximum numbers for the three derivatives.
Ijms 11 00288f12
Figure 13. Hydrogen-bond distance distribution (left) and hydrogen-bond angle distribution (right) at 60 K for amantadine, rimantadine and memantine. The most possible hydrogen-bond lengths are around 0.3 nm for the three molecules. The hydrogen-bond angles seem randomly distributed.
Figure 13. Hydrogen-bond distance distribution (left) and hydrogen-bond angle distribution (right) at 60 K for amantadine, rimantadine and memantine. The most possible hydrogen-bond lengths are around 0.3 nm for the three molecules. The hydrogen-bond angles seem randomly distributed.
Ijms 11 00288f13
Figure 14. Radial distribution functions of the seven diamondoids and derivatives (from left: Adamantane, Diamantane, Amantadine, Rimantadine, Memantine, ADM•Na, DIM•Na) in the self-assembled (solid) state.
Figure 14. Radial distribution functions of the seven diamondoids and derivatives (from left: Adamantane, Diamantane, Amantadine, Rimantadine, Memantine, ADM•Na, DIM•Na) in the self-assembled (solid) state.
Ijms 11 00288f14
Figure 15. Structure factors of the seven diamondoids and derivatives (from left: Adamantane, Diamantane, Amantadine, Rimantadine, Memantine, ADM•Na., DIM•Na) in the self-assembled (solid) state.
Figure 15. Structure factors of the seven diamondoids and derivatives (from left: Adamantane, Diamantane, Amantadine, Rimantadine, Memantine, ADM•Na., DIM•Na) in the self-assembled (solid) state.
Ijms 11 00288f15
Table 1. Molecular formulas and 3-D structures of Adamantane, Diamantane, Memantine, Rimantadine, Amantadine, Optimized ADM•Na and Optimized DIM•Na molecules. In these figures black spheres represent “–C”, whites represent “–H”, Blues represent “–N” and purples represent “–Na”.
Table 1. Molecular formulas and 3-D structures of Adamantane, Diamantane, Memantine, Rimantadine, Amantadine, Optimized ADM•Na and Optimized DIM•Na molecules. In these figures black spheres represent “–C”, whites represent “–H”, Blues represent “–N” and purples represent “–Na”.
Group 1Group 2Group 3
AdamantaneDiamantaneAmantadineRimantadineMemantineOptimized ADM•NaOptimized DIM•Na
C10H16C14H20C10H17NC11H20NC12H21NC10H15NaC14H19Na
Ijms 11 00288i1 Ijms 11 00288i2 Ijms 11 00288i3 Ijms 11 00288i4 Ijms 11 00288i5 Ijms 11 00288i6 Ijms 11 00288i7
Ijms 11 00288i8 Ijms 11 00288i9 Ijms 11 00288i10 Ijms 11 00288i11 Ijms 11 00288i12 Ijms 11 00288i13 Ijms 11 00288i14

Share and Cite

MDPI and ACS Style

Xue, Y.; Mansoori, G.A. Self-Assembly of Diamondoid Molecules and Derivatives (MD Simulations and DFT Calculations). Int. J. Mol. Sci. 2010, 11, 288-303. https://0-doi-org.brum.beds.ac.uk/10.3390/ijms11010288

AMA Style

Xue Y, Mansoori GA. Self-Assembly of Diamondoid Molecules and Derivatives (MD Simulations and DFT Calculations). International Journal of Molecular Sciences. 2010; 11(1):288-303. https://0-doi-org.brum.beds.ac.uk/10.3390/ijms11010288

Chicago/Turabian Style

Xue, Yong, and G. Ali Mansoori. 2010. "Self-Assembly of Diamondoid Molecules and Derivatives (MD Simulations and DFT Calculations)" International Journal of Molecular Sciences 11, no. 1: 288-303. https://0-doi-org.brum.beds.ac.uk/10.3390/ijms11010288

Article Metrics

Back to TopTop