Kohei Shimamura

A divide-conquer-recombine algorithmic paradigm for large spatiotemporal quantum molecular dynamics simulations.

We introduce an extension of the divide-and-conquer (DC) algorithmic paradigm called divide-conquer-recombine (DCR) to perform large quantum molecular dynamics (QMD) simulations on massively parallel supercomputers, in which interatomic forces are computed quantum mechanically in the framework of density functional theory (DFT). In DCR, the DC phase constructs globally informed, overlapping local-domain solutions, which in the recombine phase are synthesized into a global solution encompassing large spatiotemporal scales. For the DC phase, we design a lean divide-and-conquer (LDC) DFT algorithm, which significantly reduces the prefactor of the O(N) computational cost for N electrons by applying a density-adaptive boundary condition at the peripheries of the DC domains. Our globally scalable and locally efficient solver is based on a hybrid real-reciprocal space approach that combines: (1) a highly scalable real-space multigrid to represent the global charge density; and (2) a numerically efficient plane-wave basis for local electronic wave functions and charge density within each domain. Hybrid space-band decomposition is used to implement the LDC-DFT algorithm on parallel computers. A benchmark test on an IBM Blue Gene/Q computer exhibits an isogranular parallel efficiency of 0.984 on 786 432 cores for a 50.3 × 10(6)-atom SiC system. As a test of production runs, LDC-DFT-based QMD simulation involving 16 661 atoms is performed on the Blue Gene/Q to study on-demand production of hydrogen gas from water using LiAl alloy particles. As an example of the recombine phase, LDC-DFT electronic structures are used as a basis set to describe global photoexcitation dynamics with nonadiabatic QMD (NAQMD) and kinetic Monte Carlo (KMC) methods. The NAQMD simulations are based on the linear response time-dependent density functional theory to describe electronic excited states and a surface-hopping approach to describe transitions between the excited states. A series of techniques are employed for efficiently calculating the long-range exact exchange correction and excited-state forces. The NAQMD trajectories are analyzed to extract the rates of various excitonic processes, which are then used in KMC simulation to study the dynamics of the global exciton flow network. This has allowed the study of large-scale photoexcitation dynamics in 6400-atom amorphous molecular solid, reaching the experimental time scales.

Effects of twins on the electronic properties of GaAs

Generation of twin defects during the growth of semiconductor nanowires is a major concern, but their effects on electronic properties are not well understood. Here, combined quantum-mechanical and molecular-dynamics simulations reveal that the radiative decay time of an exciton increases due to twin. Furthermore, the twin-scattering contribution to electron mobility is found to be significant, in conformity with photoconductivity measurements. In addition to acting as a carrier-scattering source, twins in nanowires are found to modify the mobility by changing strain and thereby the effective mass. These effects should have profound impacts on the efficiency of nanowire-based devices.

The nature of free-carrier transport in organometal halide perovskites

Organometal halide perovskites are attracting great attention as promising material for solar cells because of their high power conversion efficiency. The high performance has been attributed to the existence of free charge carriers and their large diffusion lengths, but the nature of carrier transport at the atomistic level remains elusive. Here, nonadiabatic quantum molecular dynamics simulations elucidate the mechanisms underlying the excellent free-carrier transport in CH3NH3PbI3. Pb and I sublattices act as disjunct pathways for rapid and balanced transport of photoexcited electrons and holes, respectively, while minimizing efficiency-degrading charge recombination. On the other hand, CH3NH3 sublattice quickly screens out electrostatic electron-hole attraction to generate free carriers within 1 ps. Together this nano-architecture lets photoexcited electrons and holes dissociate instantaneously and travel far away to be harvested before dissipated as heat. This work provides much needed structure-property relationships and time-resolved information that potentially lead to rational design of efficient solar cells.

Enhanced charge recombination due to surfaces and twin defects in GaAs nanostructures

Power conversion efficiency of gallium arsenide (GaAs) nanowire (NW) solar cells is severely limited by enhanced charge recombination (CR) at sidewall surfaces, but its atomistic mechanisms are not well understood. In addition, GaAs NWs usually contain a high density of twin defects that form a twin superlattice, but its effects on CR dynamics are largely unknown. Here, quantum molecular dynamics (QMD) simulations reveal the existence of an intrinsic type-II heterostructure at the (110) GaAs surface. Nonadiabatic quantum molecular dynamics (NAQMD) simulations show that the resulting staggered band alignment causes a photoexcited electron in the bulk to rapidly transfer to the surface. We have found orders-of-magnitude enhancement of the CR rate at the surface compared with the bulk value. Furthermore, QMD and NAQMD simulations show unique surface electronic states at alternating (111)A and (111)B sidewall surfaces of a twinned [111]-oriented GaAs NW, which act as effective CR centers. The calculated large surface recombination velocity quantitatively explains recent experimental observations and provides microscopic understanding of the underlying CR processes.

Dissociation dynamics of ethylene molecules on a Ni cluster using ab initio molecular dynamics simulations.

The atomistic mechanism of dissociative adsorption of ethylene molecules on a Ni cluster is investigated by ab initio molecular-dynamics simulations. The activation free energy to dehydrogenate an ethylene molecule on the Ni cluster and the corresponding reaction rate is estimated. A remarkable finding is that the adsorption energy of ethylene molecules on the Ni cluster is considerably larger than the activation free energy, which explains why the actual reaction rate is faster than the value estimated based on only the activation free energy. It is also found from the dynamic simulations that hydrogen molecules and an ethane molecule are formed from the dissociated hydrogen atoms, whereas some exist as single atoms on the surface or in the interior of the Ni cluster. On the other hand, the dissociation of the C-C bonds of ethylene molecules is not observed. On the basis of these simulation results, the nature of the initial stage of carbon nanotube growth is discussed.

Meteorite impact-induced rapid NH3 production on early earth: Ab initio molecular dynamics simulation

NH3 is an essential molecule as a nitrogen source for prebiotic amino acid syntheses such as the Strecker reaction. Previous shock experiments demonstrated that meteorite impacts on ancient oceans would have provided a considerable amount of NH3 from atmospheric N2 and oceanic H2O through reduction by meteoritic iron. However, specific production mechanisms remain unclear, and impact velocities employed in the experiments were substantially lower than typical impact velocities of meteorites on the early Earth. Here, to investigate the issues from the atomistic viewpoint, we performed multi-scale shock technique-based ab initio molecular dynamics simulations. The results revealed a rapid production of NH3 within several picoseconds after the shock, indicating that shocks with greater impact velocities would provide further increase in the yield of NH3. Meanwhile, the picosecond-order production makes one expect that the important nitrogen source precursors of amino acids were obtained immediately after the impact. It was also observed that the reduction of N2 proceeded according to an associative mechanism, rather than a dissociative mechanism as in the Haber-Bosch process.

Nanocarbon synthesis by high-temperature oxidation of nanoparticles

High-temperature oxidation of silicon-carbide nanoparticles (nSiC) underlies a wide range of technologies from high-power electronic switches for efficient electrical grid and thermal protection of space vehicles to self-healing ceramic nanocomposites. Here, multimillion-atom reactive molecular dynamics simulations validated by ab initio quantum molecular dynamics simulations predict unexpected condensation of large graphene flakes during high-temperature oxidation of nSiC. Initial oxidation produces a molten silica shell that acts as an autocatalytic ‘nanoreactor’ by actively transporting oxygen reactants while protecting the nanocarbon product from harsh oxidizing environment. Percolation transition produces porous nanocarbon with fractal geometry, which consists of mostly sp2 carbons with pentagonal and heptagonal defects. This work suggests a simple synthetic pathway to high surface-area, low-density nanocarbon with numerous energy, biomedical and mechanical-metamaterial applications, including the reinforcement of self-healing composites.

Crystalline anisotropy of shock-induced phenomena: Omni-directional multiscale shock technique

We propose an omni-directional multiscale shock technique (OD-MSST) to study the shock waves in an arbitrary direction of crystalline materials, atomistically based on the molecular dynamics simulation method. Using OD-MSST, we found transitions from elastic to shear-banding to plastic behaviors for a model covalent crystal. In addition to such a shock “phase diagram,” a transition from inter-molecular to intra-molecular mechanochemical reaction pathways was found as a function of crystallographic orientation in an energetic van der Waals crystal.

A crossover in anisotropic nanomechanochemistry of van der Waals crystals

In nanoscale mechanochemistry, mechanical forces selectively break covalent bonds to essentially control chemical reactions. An archetype is anisotropic detonation of layered energetic molecular crystals bonded by van der Waals (vdW) interactions. Here, quantum molecular dynamics simulations reveal a crossover of anisotropic nanomechanochemistry of vdW crystal. Within 10−13 s from the passage of shock front, lateral collision produces NO2 via twisting and bending of nitro-groups and the resulting inverse Jahn-Teller effect, which is mediated by strong intra-layer hydrogen bonds. Subsequently, as we transition from heterogeneous to homogeneous mechanochemical regimes around 10−12 s, shock normal to multilayers becomes more reactive, producing H2O assisted by inter-layer N-N bond formation. These time-resolved results provide much needed atomistic understanding of nanomechanochemistry that underlies a wider range of technologies.

Ab initio study of dissociation reaction of ethylene molecules on Ni cluster

The dissociation reaction of ethylene molecules on the Ni cluster surface is investigated by ab initio molecular dynamics simulations. We observe that hydrogen atoms are generated from ethylene molecules at a rate of about 20 ps−1. The activation energy for the dissociation of a hydrogen atom is estimated to be about 0.52 eV, which corresponds to a rate of only about 0.1 ps−1. We find that the adsorption energy of an ethylene molecule on the Ni cluster is more than 1.5 eV, which is three times greater than the activation energy for the hydrogen dissociation. It is, therefore, suggested that the adsorption energy is responsible for the increase of the rate of the dissociation reaction. Based on these results, we discuss the microscopic process of the reaction of ethylene molecules on the Ni cluster in detail.