Predicting the Mechanical Properties of Carbon-Based Materials Using Molecular Dynamics

The development of innovative carbon-based materials can be greatly facilitated by molecular modeling techniques. Although molecular modeling has been used extensively to predict elastic properties of materials, modeling of more complex phenomenon such as fracture has only recently been possible with the development of new force fields such as ReaxFF, which is used in this work. It is not fully understood what molecular modeling parameters such as thermostat type, thermostat coupling, time step, system size, and strain rate are required for accurate modeling of fracture. Selection of modeling parameters to model fracture can be difficult and non-intuitive compared to modeling elastic properties using traditional force fields, and the errors generated by incorrect parameters may be non-obvious. These molecular modeling parameters are systematically investigated and their effects on the fracture of well-known carbon materials are analyzed. It is determined that for coupling coefficients of 250 fs and greater do not result in substantial differences in the stress-strain response of the materials using any thermostat type. A time step of 0.5 fs of smaller is required for accurate results. Strain rates greater than 2.2 ns-1 are sufficient to obtain repeatable results with slower strain rates for the materials studied. The results of this study indicate that further refinement of the Chenoweth parameter set is required to accurately predict the mechanical response of carbon-based systems. The ReaxFF has been used extensively to model systems in which bond breaking and formation occur. In particular ReaxFF has been used to model reactions of small molecules. Some elastic and fracture properties have been successfully modeled using ReaxFF in materials such as silicon and some metals. However, it is not clear if current parameterizations for ReaxFF are able to accurately reproduce the elastic and fracture properties of carbon materials. The stress-strain response of a new ReaxFF parameterization is compared to the previous parameterization and density functional theory results for well-known carbon materials. The new ReaxFF parameterization makes xv substantial improvements to the predicted mechanical response of carbon materials, and is found to be suitable for modeling the mechanical response of carbon materials. Finally, a new material composed of carbon nanotubes within an amorphous carbon (AC) matrix is modeled using the ReaxFF. Various parameters that may be experimentally controlled are investigated such as nanotube bundling, comparing multi-walled nanotube with single-walled nanotubes, and degree of functionalization of the nanotubes. Elastic and fracture properties are investigated for the composite systems and compared to results of pure-nanotube and pure-AC models. It is found that the arrangement of the nanotubes and degree of crosslinking may substantially affect the properties of the systems, particularly in the transverse directions.

Unraveling the Conformational Determinants of Peptide Dendrimers Using Molecular Dynamics Simulations

Peptide dendrimers are synthetic tree-like molecules composed of amino acids. There are at least two kinds of preferential structural behaviors exhibited by these molecules, which acquire either compact or noncompact shapes. However, the key structural determinants of such behaviors remained, until now, unstudied. Herein, we conduct a comprehensive investigation of the structural determinants of peptide dendrimers by employing long molecular dynamics simulations to characterize an extended set of third generation dendrimers. Our results clearly show that a trade-off between electrostatic effects and hydrogen bond formation controls structure acquisition in these systems. Moreover, by selectively changing the dendrimers charge we are able to manipulate the exhibited compactness. In contrast, the length of branching residues does not seem to be a major structural determinant. Our results are in accordance with the most recent experimental evidence and shed some light on the key molecular level interactions controlling structure acquisition in these systems. Thus, the results presented constitute valuable insights that can contribute to the development of truly tailor-made dendritic systems.

Accurate rotational constant and bond lengths of hexafluorobenzene by femtosecond rotational Raman coherence spectroscopy and ab initio calculations

The gas-phase rotational motion of hexafluorobenzene has been measured in real time using femtosecond (fs) time-resolved rotational Raman coherence spectroscopy (RR-RCS) at T = 100 and 295 K. This four-wave mixing method allows to probe the rotation of non-polar gas-phase molecules with fs time resolution over times up to ∼5 ns. The ground state rotational constant of hexafluorobenzene is determined as B 0 = 1029.740(28) MHz (2σ uncertainty) from RR-RCS transients measured in a pulsed seeded supersonic jet, where essentially only the v = 0 state is populated. Using this B 0 value, RR-RCS measurements in a room temperature gas cell give the rotational constants B v of the five lowest-lying thermally populated vibrationally excited states ν7/8, ν9, ν11/12, ν13, and ν14/15. Their B v constants differ from B 0 by between −1.02 MHz and +2.23 MHz. Combining the B 0 with the results of all-electron coupled-cluster CCSD(T) calculations of Demaison et al. [Mol. Phys.111, 1539 (2013)] and of our own allow to determine the C-C and C-F semi-experimental equilibrium bond lengths r e(C-C) = 1.3866(3) Å and r e(C-F) = 1.3244(4) Å. These agree with the CCSD(T)/wCVQZ r e bond lengths calculated by Demaison et al. within ±0.0005 Å. We also calculate the semi-experimental thermally averaged bond lengths r g(C-C)=1.3907(3) Å and r g(C-F)=1.3250(4) Å. These are at least ten times more accurate than two sets of experimental gas-phase electron diffraction r g bond lengths measured in the 1960s.

Mechanisms of recognition and binding of α-TTP to the plasma membrane by multi-scale molecular dynamics simulations

We used multiple sets of simulations both at the atomistic and coarse-grained level of resolution to investigate interaction and binding of α-tochoperol transfer protein (α-TTP) to phosphatidylinositol phosphate lipids (PIPs). Our calculations indicate that enrichment of membranes with such lipids facilitate membrane anchoring. Atomistic models suggest that PIP can be incorporated into the binding cavity of α-TTP and therefore confirm that such protein can work as lipid exchanger between the endosome and the plasma membrane. Comparison of the atomistic models of the α-TTP-PIPs complex with membrane-bound α-TTP revealed different roles for the various basic residues composing the basic patch that is key for the protein/ligand interaction. Such residues are of critical importance as several point mutations at their position lead to severe forms of ataxia with vitamin E deficiency (AVED) phenotypes. Specifically, R221 is main residue responsible for the stabilization of the complex. R68 and R192 exchange strong interactions in the protein or in the membrane complex only, suggesting that the two residues alternate contact formation, thus facilitating lipid flipping from the membrane into the protein cavity during the lipid exchange process. Finally, R59 shows weaker interactions with PIPs anyway with a clear preference for specific phosphorylation positions, hinting a role in early membrane selectivity for the protein. Altogether, our simulations reveal significant aspects at the atomistic scale of interactions of α-TTP with the plasma membrane and with PIP, providing clarifications on the mechanism of intracellular vitamin E trafficking and helping establishing the role of key residue for the functionality of α-TTP.

Parameterizing a Simple Potential Energy Function Appropriate for Vibrational Frequencies Using Molecular Dynamics

A Maple scheme for quickly parameterizing vibrational potential energy functions is presented. As an example, the potential energy function's parameters for the vibrational motions in H_2O_2 are obtained assuming the simplest potential energy function. This paper was originally written as a research paper, but rejected by the referees. It is therefore being edited into an ``educational'' paper for student usage.

Molecular Dynamics and the Diatomic Molecule as Harmonic Oscillator and Rigid Rotor

The HCl molecule is simulated (using Maple) in its dynamics, for both vibrational (and implied) rotational motions. A discussion of the center of mass transformations involved is part of the total presentation.

Molecular Dynamics Simulations of Lead and Lithium in Liquid Phase

Pb17Li is today a reference breeder material in diverse fusion R&D programs worldwide. Extracting dynamic and structural properties of liquid LiPb mixtures via molecular dynamics simulations, represent a crucial step for multiscale modeling efforts in order to understand the suitability of this compound for future Nuclear Fusion technologies. At present a Li-Pb cross potential is not available in the literature. Here we present our first results on the validation of two semi-empirical potentials for Li and Pb in liquid phase. Our results represent the establishment of a solid base as a previous but crucial step to implement a LiPb cross potential. Structural and thermodynamical analyses confirm that the implemented potentials for Li and Pb are realistic to simulate both elements in the liquid phase.

Liquid lead and Lithium. A Comprehensive molecular dynamics study

Pb17Li is today a reference breeder material in diverse fusion R&D programs worldwide. One of the main issues in these programs is the problem of liquid metals breeder blanket behavior. Structural material of the blanket should meet high requirements because of extreme operating conditions. Therefore the knowledge of eutectic properties like optimal composition, physical and thermodynamic behavior or diffusion coefficients of Tritium are extremely necessary for current designs. In particular, the knowledge of the function linking the tritium concentration dissolved in liquid materials with the tritium partial pressure at a liquid/gas interface in equilibrium, CT=f(PT), is of basic importance because it directly impacts all functional properties of a blanket determining: tritium inventory, tritium permeation rate and tritium extraction efficiency. Nowadays, understanding the structure and behavior of this compound is a real goal in fusion engineering and materials science. Simulations of liquids can provide much information to the community; not only supplementing experimental data, but providing new tests of theories and ideas, making specific predictions that require experimental tests, and ultimately helping to lead to the deeper understanding and better predictive behavior.