Source author record

Adri C. T. van Duin

Adri C. T. van Duin appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

cond-mat.mtrl-sci physics.chem-ph cond-mat.mes-hall physics.app-ph physics.atm-clus

Catalog footprint

What is connected

9works

5topics

4close collaborators

Actions

Connect this record

Open graph Browse works

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

preprint2022arXiv

On the origin of non-classical ripples in draped graphene sheets

Ever since the discovery of graphene and subsequent explosion of interest in single atom thick materials, studying their mechanical properties has been an active area of research. New length scales often necessitate a rethinking of physical laws, making such studies crucial for understanding and ultimately utilizing novel material properties. Here we report on the investigation of nanoscale periodic ripples in suspended, single layer graphene sheets by scanning tunneling microscopy and atomistic scale simulations. Unlike the sinusoidal ripples found in classical fabrics, we find that graphene forms triangular ripples, where bending is limited to a narrow region on the order of a few unit cell dimensions at the apex of each ripple. This non-classical bending profile results in graphene behaving like a bizarre fabric, which regardless of how it is draped, always buckles at the same angle. Investigating the origin of such non-classical mechanical properties, we find that unlike a thin classical fabric, both in-plane and out-of-plane deformations occur in a graphene sheet. These two modes of deformation compete with each other, resulting in a strain-locked optimal buckling configuration when draped. Electronically, we see that this in-plane deformation generates pseudo electric fields creating a ~3 nm wide pnp heterojunction purely by strain modulation.

preprint2021arXiv

Atomistic insights into the degradation of halide perovskites: a reactive force field molecular dynamics study

Halide perovskites make efficient solar cells due to their exceptional optoelectronic properties, but suffer from several stability issues. The characterization of the degradation processes is challenging because of the limitations in the spatio-temporal resolution in experiments and the absence of efficient computational methods to study the reactive processes. Here, we present the first effort in developing reactive force fields for large scale molecular dynamics simulations of the phase instability and the defect-induced degradation reactions in inorganic CsPbI$_{3}$. We find that the phase transitions are driven by a combination of the anharmonicity of the perovskite lattice with the thermal entropy. At relatively low temperatures, the Cs cations tend to move away from the preferential positions with good contacts with the surrounding metal halide framework, potentially causing its conversion to a non-perovskite phase. Our simulations of defective structures reveal that, although both iodine vacancies and interstitials are very mobile in the perovskite lattice, the vacancies have a detrimental effect on the stability, initiating the decomposition reactions of perovskites to PbI$_{2}$. Our work puts ReaxFF forward as an effective computational framework to study reactive processes in halide perovskites.

preprint2021arXiv

ReaxFF Reactive Force Field Development for Cu/Si Systems and application to Copper Cluster Formation During Cu Diffusion Inside Silicon

Transition metal impurities such as nickel, copper, and iron, in solid-state materials like silicon have a significant impact on the electrical performance of integrated circuits and solar cells. To study the impact of copper impurities inside bulk silicon on the electrical properties of the material, one needs to understand the configurational space of copper atoms incorporated inside the silicon lattice. In this work, we performed ReaxFF reactive force field based molecular dynamics simulations, studying different configurations of individual and crystalline copper atoms inside bulk silicon by looking at the diffusional behavior of copper in silicon. The ReaxFF Cu/Si parameter set was developed by training against DFT data, including the energy barrier for an individual Cu-atom inside a silicon lattice. We found that the diffusion of copper atoms has a direct relationship with the temperature. Moreover, it is also shown that individual copper atoms start to clusterize inside bulk silicon at elevated temperatures. Our simulation results provide a comprehensive picture of the effects of temperature and copper concentration on the crystallization of individual copper inside silicon lattice. Finally, the stress-strain relationship of Cu/Si compounds under uniaxial tensile loading have been obtained. Our results indicate a decrease in the elastic modulus with increasing level of Cu-impurity concentration. We observe spontaneous microcracking of the Si during the stress-strain tests as a consequence of the formation of a small Cu clusters adjacent to the Si surface.

preprint2020arXiv

Polyacrylonitrile/Graphene Nanocomposite: Towards the Next Generation of Carbon Fibers

Carbon Fibers (CFs) are the key solution for the future lightweight vehicle with enhanced fuel efficiency and reduced emissions owing to their ultrahigh strength to weight ratio. However, the high cost of the current dominant PAN-based CFs hinders their application. The use of low-cost alternative precursors may overcome this issue. Unfortunately, low-cost CFs derived from cheaper single component precursors suffer from poor mechanical properties. Developing composite CFs by adding nanoadditives is very promising for low-cost CFs. Therefore, a fundamental understanding of carbonization condition impacts and polymer/additives conversion mechanisms during whole CF production are essential to develop low-cost CFs. In this work, we have demonstrated how the carbonization temperature affects the PAN/graphene CFs properties by performing a series of ReaxFF based molecular dynamics simulations. We found that graphene edges along with the nitrogen and oxygen functional groups have a catalytic role and act as seeds for the graphitic structure growth. Our MD simulations unveil that the addition of the graphene to PAN precursor modifies all-carbon membered rings in CFs and enhances the alignments of 6-member carbon rings in carbonization which leads to superior mechanical properties compare to PAN-based CFs. These ReaxFF simulation results are validates by experimental structural and mechanical characterizations. Interestingly, mechanical characterizations indicate that PAN/graphene CFs carbonized at 1250 C demonstrate 90.9% increase in strength and 101.9% enhancement in Young's modulus compare to the PAN-based CFs carbonized at 1500 C. The superior mechanical properties of PAN/graphene CFs at lower carbonization temperatures offers a path to both energy savings and cost reduction by decreasing the carbonization temperature and could provide key insights for the development of low-cost CFs.

preprint2020arXiv

ReaxFF reactive force field study of polymerization of polymer matrix in carbon nanotube-composite system

Human transport to Mars and deep space explorations demand the development of new materials with extraordinary high performance-to-mass ratios. Promising candidates to fulfill these requirements are ultrahigh strength lightweight (UHSL) materials, which consist of polymer matrices fortified by pristine carbon nanotubes (CNTs). Previous investigations have showed that with an increase in CNT diameter, its preferred configuration changes from a circular form to a flattened shape that can be obtained in high pressure or tension conditions. The ReaxFF reactive force field can reveal detailed chemical interactions at the atomistic scale. To enable ReaxFF simulations on CNT/polymer interfaces, we trained force field parameters to capture the proper structure of flattened carbon nanotubes (flCNTs), i.e. dumbbell-like shape CNTs, against available polymer consistent force field -- interface force field (PCFF-IFF) data which had good proximity to density functional theory (DFT) data. In this study we used accelerated ReaxFF molecular dynamics simulation using the optimized force field to study the polymerization of diglycidyl ether of bisphenol F (Bis F) and diethyltoluenediamine (DEDTA) molecules in vicinity of circular and flattened CNTs. Our results indicate that the flat regions of flCNT are more favorable spots for the polymers to settle compared to curved regions due to higher binding energies. Moreover, higher dimer generation around flCNT results in more effective coating of the CNT which leads to higher load transfer in compared to circular CNT. According to our results there is a high alignment between polymers and CNT surface which is due to strong pi-pi interactions of aromatic carbon rings in the polymers and CNTs. These ReaxFF simulations indicate the capability of this method to simultaneously observe the polymerization of monomers along with their interactions with CNTs.

preprint2015arXiv

Lithium Ion Solvation and Diffusion in Bulk Organic Electrolytes from First Principles and Classical Reactive Molecular Dynamics

Lithium-ion battery performance is strongly influenced by the ionic conductivity of the electrolyte, which depends on the speed at which Li ions migrate across the cell and relates to their solvation structure. The choice of solvent can greatly impact both solvation and diffusivity of Li ions. We use first principles molecular dynamics to examine the solvation and diffusion of Li ions in the bulk organic solvents ethylene carbonate (EC), ethyl methyl carbonate (EMC), and a mixture of EC/EMC. We find that Li ions are solvated by either carbonyl or ether oxygen atoms of the solvents and sometimes by the PF$_6^-$ anion. Li$^+$ prefers a tetrahedrally-coordinated first solvation shell regardless of which species are involved, with the specific preferred solvation structure dependent on the organic solvent. In addition, we calculate Li diffusion coefficients in each electrolyte, finding slightly larger diffusivities in the linear carbonate EMC compared to the cyclic carbonate EC. The magnitude of the diffusion coefficient correlates with the strength of Li$^+$ solvation. Corresponding analysis for the PF$_6^-$ anion shows greater diffusivity associated with a weakly-bound, poorly defined first solvation shell. These results may be used to aid in the design of new electrolytes to improve Li-ion battery performance.

preprint2014arXiv

Stabilized silicene within bilayer graphene: A proposal based on molecular dynamics and density-functional tight-binding calculations

Free standing silicene is predicted to display comparable electronic properties as graphene. However, the yet synthesized silicene-like structures have been only realized on different substrates which turned out to exhibit versatile crystallographic structures that are very different from the theoretically predicted buckled phase of freestanding silicene. This calls for a different approach where silicene is stabilized using very weakly interacting surfaces. We propose here a novel route by using graphene bilayer as a scaffold. The confinement between the flat graphene layers results in a planar clustering of Si atoms with small buckling, which is energetically unfavorable in vacuum. Buckled hexagonal arrangement of Si atoms similar to free-standing silicene is observed for large clusters, which, in contrast to Si atoms on metallic surfaces, is only very weakly van der Waals coupled to the graphene layers. These clusters are found to be stable well above room temperature. Our findings, which are supported by density functional tight-binding calculations, show that intercalating bilayer graphene with Si is a favorable route to realize silicene.

preprint2013arXiv

Thermal properties of fluorinated graphene

Large scale atomistic simulations using the reactive force field approach (ReaxFF) are implemented to investigate the thermomechanical properties of fluorinated graphene (FG). A new set of parameters for the reactive force field potential (ReaxFF) optimized to reproduce key quantum mechanical properties of relevant carbon-fluor cluster systems are presented. Molecular dynamics (MD) simulations are used to investigate the thermal rippling behavior of FG and its mechanical properties and compare them with graphene (GE), graphane (GA) and a sheet of BN. The mean square value of the height fluctuations $< h^2>$ and the height-height correlation function $H(q)$ for different system sizes and temperatures show that FG is an un-rippled system in contrast to the thermal rippling behavior of graphene (GE). The effective Young's modulus of a flake of fluorinated graphene is obtained to be 273 N/m and 250 N/m for a flake of FG under uniaxial strain along arm-chair and zig-zag direction, respectively.

preprint2011arXiv

Reactive force field for lithium-aluminum silicates with applications to eucryptite phases

We have parameterized a reactive force field (ReaxFF) for lithium aluminum silicates using density functional theory (DFT) calculations of structural properties of a number of bulk phase oxides, silicates, and aluminates, as well as of several representative clusters. The force field parameters optimized in this study were found to predict lattice parameters and heats of formation of selected condensed phases in excellent agreement with previous DFT calculations and with experiments. We have used the newly developed force-field to study the eucryptite phases in terms of their thermodynamic stability and their elastic properties. We have found that (a) these ReaxFF parameters predict the correct order of stability of the three crystalline polymorphs of eucryptite, α, β, and γ, and (b) that upon indentation, a new phase appears at applied pressures \geq 7 GPa. The high pressure phase obtained upon indentation is amorphous, as illustrated by the radial distribution functions calculated for different pairs of elements. In terms of elastic properties analysis, we have determined the elements of the stiffness tensor for α- and β- eucryptite at the level of ReaxFF, and discussed the elastic anisotropy of these two polymorphs. Polycrystalline average properties of these eucryptite phases are also reported to serve as ReaxFF predictions of their elastic moduli (in the case of α-eucryptite), or as tests against values known from experiments or DFT calculations (β- eucrypite). The ReaxFF potential reported here can also describe well single-species systems (e.g., Li-metal, Al-metal, and condensed phases of silicon), which makes it suitable for investigating struc- ture and properties of suboxides, atomic-scale mechanisms responsible for phase transformations, as well as oxidation-reduction reactions.