Ab Initio Simulations Research Articles

New tetragonal phases of titanium under shock loading - Predicted by molecular dynamics and ab initio calculations

Based on molecular dynamics simulations of the long-term evolution of polycrystalline titanium structures under shock waves, a novel phase (ζ) has been identified to emerge at shock pressures of approximately 20 GPa. This phase originates from the initial HCP phase and the disordered structure induced by the shock. The energy change characteristics during this phase transition are explained based on the Helmholtz free energy under volume constraints. The ζ phase comprises framework atoms (F-atoms) with fixed positions and unfixed interstitial atoms (I-atoms), exhibiting incomplete symmetry. Validation of this phase structure using ab initio simulations and machine learning potentials supports its existence as a metastable structure. The transformation from the initial α phase to ζ phase is diffusional transformation, while the reverse transformation is multistage displacive transformation.

Crystal orbital overlap population based on all-electron ab initio simulation with numeric atom-centered orbitals and its application to chemical-bonding analysis in Li-intercalated layered materials

Crystal orbital overlap population (COOP) is one of the effective tools for chemical-bonding analysis, and thus it has been utilized in the materials development and characterization. In this study, we developed a code to perform the COOP-based chemical-bonding analysis based on the wave function obtained from a first principles all-electron calculation with numeric atom-centered orbitals. The chemical-bonding analysis using the developed code was demonstrated for F2, Si, CaC6, and metals including Ti and Nb. Furthermore, we applied the method to analyze the chemical-bonding changes associated with a Li intercalation in three representative layered materials: graphite, MoS2, and ZrNCl, because of their great industrial importance, particularly for the applications in battery and superconducting materials. The COOP analysis provided some insights for understanding the intercalation mechanism and the stability of the intercalated materials from a chemical-bonding viewpoint.

The anodic dissolution kinetics of Mg alloys in water based on ab initio molecular dynamics simulations

AbstractThe corrosion susceptibility of magnesium (Mg) alloys presents a significant challenge for their broad application. Although there have been extensive experimental and theoretical investigations, the corrosion mechanisms of Mg alloys are still unclear, especially the anodic dissolution process. Here, a thorough theoretical investigation based on ab initio molecular dynamics and metadynamics simulations has been conducted to clarify the underlying corrosion mechanism of Mg anode and propose effective strategies for enhancing corrosion resistance. Through comprehensive analyses of interfacial structures and equilibrium potentials for Mg(0001)/H2O interface models with different water thicknesses, the Mg(0001)/72 H2O model is identified to be reasonable with −2.17 V vs. standard hydrogen electrode equilibrium potential. In addition, utilizing metadynamics, the free energy barrier for Mg dissolution is calculated to be 0.835 eV, enabling the theoretical determination of anodic polarization curves for pure Mg that aligns well with experimental data. Based on the Mg(0001)/72 H2O model, we further explore the effects of various alloying elements on anodic corrosion resistance, among which Al and Mn alloying elements are found to enhance corrosion resistance of Mg. This study provides valuable atomic‐scale insights into the corrosion mechanism of magnesium alloys, offering theoretical guidance for developing novel corrosion‐resistant Mg alloys.

Coordination number impact on magnetic properties of Schiff base Co(II) complexes

Two mononuclear Co(II) complexes with a tetrahedral and octahedral coordination geometry with a Schiff base ligand L− (anion of 4-[(2-fluoroanilino)-phenyl-methylene]-5-methyl-2-phenyl-pyrazol-3-one) having the formulas [CoL2]·2MeCN (1) and [CoL2(H2O)2]·2MeCN (2) have been synthesized and structurally characterized. Notably, their coordination geometry and magnetic relaxation behaviors can be tuned by variation of reaction conditions. Static magnetic studies reveal that complex 1 exhibits considerable easy-axis magnetic anisotropy. Frequency-dependent susceptibility measurements revealed slow magnetic relaxation behavior with spin-reversal barriers of 79.2 K at an applied external field of 500 Oe. Clear magnetic hysteresis loops were observed for 1 at 2 K. Complex 2 show field-induced slow magnetic relaxation with Raman and direct relaxation processes. The observed magnetic behavior of 1 and 2 is confirmed by high-level ab initio calculations and simulation of χMT(T) curves.

Adsorption of Glycine on TiO2 in Water from On-the-fly Free-Energy Calculations and In Situ Electrochemical Impedance Spectroscopy.

We report here an experimental-computational study of hydrated TiO2 anatase nanoparticles interacting with glycine, where we obtain quantitative agreement of the measured adsorption free energies. Ab initio simulations are performed within the tight binding and density functional theory in combination with enhanced free-energy sampling techniques, which exploit the thermodynamic integration of the unbiased mean forces collected on-the-fly along the molecular dynamics trajectories. The experiments adopt a new and efficient setup for electrochemical impedance spectroscopy measurements based on portable screen-printed gold electrodes, which allows fast and in situ signal assessment. The measured adsorption free energy is -30 kJ/mol (both from experiment and calculation), with preferential interaction of the charged NH3+ group which strongly adsorbs on the TiO2 bridging oxygens. This highlights the importance of the terminal amino groups in the adsorption mechanism of amino acids on hydrated metal oxides. The excellent agreement between computation and experiment for this amino acid opens the doors to the exploration of the interaction free energies for other moderately complex bionano systems.

Pressure effects on the structural, electronic, elastic, optical, and vibrational properties of YMg intermetallic compounds: a first-principles study

The properties of YMg in B2 structure have been comprehensively analysed using the first-principles plane-wave pseudopotential method. Specifically, the structural, electronic, elastic, vibrational, and optical properties were investigated using the generalized gradient approximation (GGA) method in the context of density functional theory. The Vienna ab initio simulation package (VASP) was utilized for these calculations. The computed lattice parameter (3.803 Å) and bulk modulus (41.33 GPa) are consistent with the earlier data on ambient pressure. The electronic band structure and energy-dependent density of states reveal the metallic nature of the titled compounds. The Born stability requirements confirmed the mechanical stability. The analysis of Pugh’s and Poisson’s ratios and Cauchy’s pressure reveals that YMg is ductile under the pressures in consideration. According to several anisotropy indices, the compound is noticeably anisotropic both in ambient and under pressure. Our investigation includes an analysis of several fundamental mechanical parameters of the material, including the bulk modulus, the pressure derivative of the Zener anisotropy factor, Poisson’s ratio, isotropic shear modulus and Young’s modulus with a particular focus on their dependence on pressure. We have determined that the elastic constants obtained remain mechanically stable, satisfying the Born Stability conditions even at high pressures of up to 60 GPa. To explore the dynamic stability of YMg, we analysed the material’s phonon dispersion curves. The examined compound displays stability under dynamic conditions from 0 GPa to 30 GPa, as evidenced by its positive vibration frequencies. However, this stability is not sustained under higher pressure, as the compound becomes unstable after 30 GPa to 70 GPa. The electronic band structure and density of states diagrams demonstrate YMg’s metallic properties. At atmospheric pressure (0 GPa), the total density of states (TDOS) near the Fermi level is approximately 1.63 states/eV, with pressure application reducing DOS. The dielectric function, refractive index, and energy loss spectra are examined within the 0–20 eV energy range. YMg exhibits its highest absorption between 4 and 11 eV. The peak optical conductivity is observed around 0.78 eV (equivalent to 1589.5409 nm), while the most significant energy loss occurs at 11.90 eV, roughly corresponding to 2.8 Hz in the ultraviolet spectrum. Moreover, we extensively analyzed the material’s phonon thermodynamic and optical properties, providing insights into its behavior under various conditions. The outcomes acquired at zero pressure are generally coherent with the current theoretical values.

Reduction of precious metal ions in aqueous solutions by contact-electro-catalysis

Precious metals are core assets for the development of modern technologies in various fields. Their scarcity poses the question of their cost, life cycle and reuse. Recently, an emerging catalysis employing contact-electrification (CE) at water-solid interfaces to drive redox reaction, called contact-electro-catalysis (CEC), has been used to develop metal free mechano-catalytic methods to efficiently degrade refractory organic compounds, produce hydrogen peroxide, or leach metals from spent Li-Ion batteries. Here, we show ultrasonic CEC can successfully drive the reduction of Ag(ac), Rh3+, [PtCl4]2-, Ag+, Hg2+, Pd2+, [AuCl4]-, and Ir3+, in both anaerobic and aerobic conditions. The effect of oxygen on the reaction is studied by electron paramagnetic resonance (EPR) spectroscopy and ab-initio simulation. Combining measurements of charge transfers during water-solid CE, EPR spectroscopy and gold extraction experiments help show the link between CE and CEC. What’s more, this method based on water-solid CE is capable of extracting gold from synthetic solutions with concentrations ranging from as low as 0.196 ppm up to 196 ppm, reaching in 3 h extraction capacities ranging from 0.756 to 722.5 mg g−1 in 3 h. Finally, we showed CEC is employed to design a metal-free, selective, and recyclable catalytic gold extraction methods from e-waste aqueous leachates.

Unraveling the Phosphorus Adsorption Mechanisms in Three-dimensional Reduced Graphene Oxide Materials.

To prevent eutrophication, controlling the phosphate concentration levels is one of the most important issues in surface water management. One of the most utilized methods is phosphate adsorption. However, its application faces a bottleneck due to the unclear understanding of adsorption and interaction mechanisms. The present work unlocks the phosphorus adsorption mechanisms in three-dimensional reduced graphene oxide with different reduction levels and pore sizes to remove phosphate from water using experiments and multiscale simulations. Experiments were performed to evaluate the influence of pH, ionic strength, and temperature on the adsorption. Molecular Dynamics and Ab Initio simulations evaluated the influence of the pore size and oxidation degrees of the materials. We show that the adsorption capacity of the materials increases with increasing pH and ionic strength and decreasing temperature. It is observed that the more oxidized the material and the less compact the structure, the better the adsorption. These results are theoretically explained in terms of the interaction of functional groups and the clustering of phosphate ions, which results in better adsorption in materials with larger pores. The underlying mechanisms for the 3D-reduced graphene oxide performance were confirmed by spectroscopy analysis. All the results show that 3D-reduced graphene oxide can sorb phosphate in different complex water remediation systems with characteristics that can be modulated by changing the material synthesis method.

Surface-Enriched Room-Temperature Liquid Bismuth for Catalytic CO2 Reduction.

Bismuth-based electrocatalysts are effective for carbon dioxide (CO2) reduction to formate. However, at room temperature, these materials are only available in solid state, which inevitably suffers from surface deactivation, declining current densities, and Faradaic efficiencies. Here, the formation of a liquid bismuth catalyst on the liquid gallium surface at ambient conditions is shown as its exceptional performance in the electrochemical reduction of CO2 (i.e., CO2RR). By doping a trace amount of bismuth (740ppm atomic) in gallium liquid metal, a surface enrichment of bismuth by over 400 times (30 at%) in liquid state is obtained without atomic aggregation, achieving 98% Faradic efficiency for CO2 conversion to formate over 80 h. Ab initio molecular simulations and density functional theory calculations reveal that bismuth atoms in the liquid state are the most energetically favorable sites for the CO2RR intermediates, superior to solid Bi-sites, as well as joint GaBi-sites. This study opens an avenue for fabricating high-performing liquid-state metallic catalysts that cannot be reached by elementary metals under electrocatalytic conditions.

Unraveling the Hydration Shell Structure and Dynamics of Group 10 Aqua Ions.

We present ab initio simulations based on subsystem DFT of group 10 aqua ions accurately compared against experimental data on hydration structure. Our simulations provide insights into the molecular structures and dynamics of hydration shells, offering recalibrated interpretations of experimental results. We observe a soft, but distinct second hydration shell in Palladium (Pd) due to a balance between thermal fluctuations, metal-water interactions, and hydrogen bonding. Nickel (Ni) and platinum (Pt) exhibit more rigid hydration shells. Notably, our simulations align with experimental findings for Pd, showing axial hydration marked by a broad peak at about 3 Å in the Pd-O radial distribution function, revising the previously sharp "mesoshell" prediction. We introduce the "hydrogen bond dome" concept to describe a resilient network of hydrogen-bonded water molecules around the metal, which plays a critical role in the axial hydration dynamics.

Phase segregation and nanoconfined fluid O2 in a lithium-rich oxide cathode

Lithium-rich oxide cathodes lose energy density during cycling due to atomic disordering and nanoscale structural rearrangements, which are both challenging to characterize. Here we resolve the kinetics and thermodynamics of these processes in an exemplar layered Li-rich (Li1.2–xMn0.8O2) cathode using a combined approach of ab initio molecular dynamics and cluster expansion-based Monte Carlo simulations. We identify a kinetically accessible and thermodynamically favourable mechanism to form O2 molecules in the bulk, involving Mn migration and driven by interlayer oxygen dimerization. At the top of charge, the bulk structure locally phase segregates into MnO2-rich regions and Mn-deficient nanovoids, which contain O2 molecules as a nanoconfined fluid. These nanovoids are connected in a percolating network, potentially allowing long-range oxygen transport and linking bulk O2 formation to surface O2 loss. These insights highlight the importance of developing strategies to kinetically stabilize the bulk structure of Li-rich O-redox cathodes to maintain their high energy densities.

Water Solvent Reorganization upon Ultrafast Resonant Stimulated X-ray Raman Excitation of a Metalloporphyrin Dimer.

We propose an X-ray Raman pump-X-ray diffraction probe scheme to follow solvation dynamics upon charge migration in a solute molecule. The X-ray Raman pump selectively prepares a valence electronic wavepacket in the solute, while the probe provides information about the entire molecular ensemble. A combination of molecular dynamics and ab initio quantum chemistry simulations is applied to a Zn-Ni porphyrin dimer in water. Using time-resolved X-ray diffraction and pair distribution functions, we extracted solvation shell dynamics.

Neutron scattering and neural-network quantum molecular dynamics investigation of the vibrations of ammonia along the solid-to-liquid transition

Vibrational spectroscopy allows us to understand complex physical and chemical interactions of molecular crystals and liquids such as ammonia, which has recently emerged as a strong hydrogen fuel candidate to support a sustainable society. We report inelastic neutron scattering measurement of vibrational properties of ammonia along the solid-to-liquid phase transition with high enough resolution for direct comparisons to ab-initio simulations. Theoretical analysis reveals the essential role of nuclear quantum effects (NQEs) for correctly describing the intermolecular spectrum as well as high energy intramolecular N-H stretching modes. This is achieved by training neural network models using ab-initio path-integral molecular dynamics (PIMD) simulations, thereby encompassing large spatiotemporal trajectories required to resolve low energy dynamics while retaining NQEs. Our results not only establish the role of NQEs in ammonia but also provide general computational frameworks to study complex molecular systems with NQEs.

Uniform Tendency of Surface Dipoles Across Silicon Doping Levels and Types of H‐Terminated Surfaces

AbstractThe termination of surface‐dangling bonds on silicon through hydrogen atoms, also known as Si–H, can achieve chemical passivation and reduce surface states in the electronic bandgap, thus altering electronic properties. Through a comprehensive study of doping levels (1014–1020 cm−3) and types (n and p), a consistent surface dipole trend induced by Si–H termination is discovered. It is achieved by redistributing surface charges and establishing thermal equilibrium with the chemical bond. To resolve this, the surface work function, surface electron affinity, and the energy difference between the valence band and the Fermi level are measured by employing the Kelvin probe, X‐ray photoelectron spectroscopy, and photoelectron yield spectroscopy methods. These findings are further validated through ab initio simulations. This finding has immense implications not only for eliminating electronic defects at semiconductor interfaces, which is crucial in microelectronics but also for developing and engineering hybrid interfaces and heterojunctions with controlled electronic properties.

Performance assessment of monolayer Black Phosphorus DG-JLFET

Two-dimensional materials are very promising for ultra-short channel future devices. This paper investigates, for the first time, the viability of the junction less transistors based on 2-D materials for future state-of-the-art technology nodes. Specifically, we investigate the performance of a Junction less monolayer black phosphorus (BP) FET (JLFET) with 12nm gate length using ab initio quantum transport simulations. The electrostatic control mechanism of the device and various intrinsic static and dynamic characteristics of the device are studied. The results reveal that BP JLFET performance can fulfill the benchmark requirement of International Roadmap for Devices and Systems (IRDS 2021) for 2028 in terms of HP and HD applications. Therefore, JLFET based on 2D materials can be a promising alternative for nano scale future device applications.

Machine-Learning driven STM images prediction of doped/defective graphene: Towards optimized tools for 2D nanomaterials characterization

Scanning tunneling microscopy (STM) is a key characterization technique that allows for visualization of specific features at nanostructures, given its ability to reveal changes in local charge densities in doping or defective sites. Besides the experimental acquisition of STM data from real samples, there is the theoretical route, which allows for STM images simulation from ab-initio calculations on defined nanostructures. This work presents an alternative route for theoretical STM image prediction by means of Machine Learning (ML) methods, based on the training of deep convolutional neural networks using simulated STM data. The ML-based STM predictions for defective, B- and N-doped graphene quantum dots are compared to ab-initio STM simulations, in terms of accuracy of structural rearrangements, charge density localizations, and computational resources requirements. The results demonstrate that STM images predicted by ML-methods are accurate at doping sites, whereas they show good similarity to simulations of large structural vacancies. This work represents a novel alternative for the use of ML methods in image-based characterization routines of doped/defective graphene, and potentially, for other 2-dimensional nanomaterials.

High-performance Schottky-barrier field-effect transistors based on two-dimensional GaN with Ag or Au contacts

The performance of two-dimensional (2D) monolayer GaN Schottky barrier field effect transistors (SBFETs) with four different metals (Ag, Au, Al, and Pt) as electrodes were studied by using ab initio simulations. The N-type Schottky contact is formed on Ag-GaN, Au-GaN, and Pt-GaN, characterized by electron Schottky barrier heights (SBH) of 0.6, 0.5, and 0.38 eV, respectively. Whereas, Al-GaN contact formed P-type Schottky contact with hole SBH of 1.42 eV. The 5.1 nm GaN SBFETs with four metal electrodes all could overcome the short-channel effect. Additionally, GaN SBFETs with Ag and Au electrodes have excellent performance, whose ON-currents are 1151.1 μA/μm and 1258.9 μA/μm, respectively. They could satisfy the demands of International Technology Roadmap for Semiconductors for high-performance transistor. Furthermore, research indicates that the device's current increases with increasing temperature. Notably, under the constant bias and gate voltage, the current is unaffected by temperature variations between the left and right electrodes.

Astro-electrochemistry of NH3 clusters and ice: e− trapping, stability, and electron transfer

ABSTRACT Quantifying electron trapping and transfer to small molecules is crucial for interfacial chemistry. In an astrochemical context, we study how NH3 clusters in both crystalline and amorphous forms can capture low-energy electrons to form ammoniated electrons. Electron affinities, vertical detachment energies, and vertical attachment energies were computed via ab initio static and dynamics simulations, (DFT, DLPNO-CCSD(T);AIMD), for (NH3)n clusters (n = 4, 5, 6, 8, 14, 23, and 38). Our results indicate that the clusters could trap and stabilize the unpaired electron which is always externally localized on the clusters. Interactions of the ammoniated electron clusters with astrochemically relevant molecules indicate that electron transfer to water and methanol are feasible, forming the radical anions (H2O)−· and (CH3OH)−·. The trapping of electrons by both crystalline and amorphous NH3 ices, and subsequent transfer to small molecules, highlights ‘astro-electrochemical’ reactions, and has implications for both astrochemistry as well as terrestrial cluster science.

The basic intermediate-range order structures in the FLiNaK–LaF3–Li2O melt: identification and characterisation

ABSTRACT In this paper, we investigate the structure of FLiNaK–LaF3 melt with small addition of Li2O through the theoretical approach combining ab initio molecular dynamics and simulation of the isolated structural units. The structure of the melt was studied in detail by calculation of radial distribution functions along with coordination numbers and radii. We found that oxyfluoride intermediate-range structure is formed in the melt. Our observations support the suggestions of other authors on formation of intermediate-range structures in fluoride melts containing metal oxides. Given structure shows long lifetime, the consideration of the isolated specie is valuable. We performed a simulation of ‘gas-phase’ [O2La4Fn] structure in order to calculate its Raman spectrum. It was found that the spectrum of the isolated structure agree in semi-quantitative manner with the vibrational density of states of the oxygen ions in the melt.

Thermoelectric properties and thermal transport in two-dimensional GaInSe3 and GaInTe3 monolayers: A first-principles study

We here report the electronic structure calculation of GaInSe3 and GaInTe3 monolayers with the P3m1 (no. 156) space group. The electronic structure and thermoelectric properties of the monolayers are calculated through the Vienna Ab initio Simulation Package and BoltzTraP2 codes. The dynamic and thermodynamic stabilities were verified by calculating their phonon spectra and simulating ab initio molecular dynamics. The monolayers were found to have a direct bandgap, with both PBE + SOC and HSE06 + SOC potentials. The lattice thermal conductivity of GaInTe3 monolayer calculated using Phono3py code shows ultra-low values due to enhanced phonon–phonon scattering. Combining electrical and thermal transport, the values have been evaluated. Importantly, the p-type GaInTe3 has excellent thermoelectric properties at 700 K, with a zT value of 2, indicating that the p-type GaInTe3 has potential application in the field of thermoelectricity.