First-principles Calculations Research Articles

Acetonitrile real gas phase behavior from quasi-ideal gas to nanodroplets: A microscopical view.

We have applied a recently developed general purpose acetonitrile force field based on first-principles calculations to simulate acetonitrile in the gas phase at different temperatures and densities. These conditions range from nearly ideal to real gas phase behavior and condensation. The molecular dynamics simulation results agree fairly well with the experimental studies available in the literature on the gas samples. The structural analysis of aggregates and their associated interaction energies is examined and related to the early model proposed on molecular association and equilibrium determining the non-ideal behavior. The formation of dimers is mainly responsible for the non-ideal behavior of the gas at very low density, confirming suggested models based on previous experimental studies. However, when the density of the sample rises, the level of aggregation increases and the simple concept of dimerization does not hold anymore. The real behavior adopted by the gas is related to the distribution of molecular structures observed. The macroscopical view of a real gas as a generic interparticle interaction system without a defined form may then be rationalized on the basis of a defined molecular association originated by a distribution of aggregates at the low density regime. The sample with the highest density (∼1.4 × 103 molm-3) at the lowest temperature exhibits a massive aggregation where most of the acetonitrile (ACN) molecules in the simulation box form a big cluster. Its radial distribution function is similar to that of the liquid ACN. This strongly inhomogeneous distribution in the box can be considered a condensation in the gas phase under specific density-T conditions. This formation opens the door to the potential tuning of its solvent properties as a function of its size in these nanodroplets that in turn are controlled by the density-T conditions.

Evaluation of the Exchange Stiffness Constants of Itinerant Magnets at Finite Temperatures from the First-Principles Calculations

Evaluation of the Exchange Stiffness Constants of Itinerant Magnets at Finite Temperatures from the First-Principles Calculations

Vertical Dipole Dominates Charge Carrier Lifetime in Monolayer Janus MoSSe.

Two-dimensional semiconductor materials with vertical dipoles are promising photocatalysts as vertical dipoles not only promote the electron-hole separation but also enhance the carrier redox ability. However, the influence of vertical dipoles on carrier recombination in such materials, especially the competing relationship between vertical dipoles and band gaps, is not yet clear. Herein, first-principles calculations and nonadiabatic molecular dynamics simulations were combined to clarify the influence of band gap and vertical dipole on the carrier lifetime in Janus MoSSe monolayer. By comparing with the results of MoS2 and MoSe2 as well as exploring the carrier lifetime of MoSSe under strain regulation, it has been demonstrated that the vertical dipole, rather than the band gap, is the dominant factor affecting the carrier lifetime. Strikingly, a linear relationship between the carrier lifetime and vertical dipole is revealed. These findings have important implications for the design of high-performance photocatalysts and optoelectronic devices.

Realistic Simulation of Valence-Band Photoemission Spectra by Using First-Principles Calculations on the Basis of FLAPW Method

Realistic Simulation of Valence-Band Photoemission Spectra by Using First-Principles Calculations on the Basis of FLAPW Method

Effect of non-metal doping on the optoelectronic properties of ZrS2/ZrSe2 heterostructure under strain: a first-principles study.

In this paper, we systematically studied the effects of non-metallic element (B, C, N, O, F) doping and biaxial stretching on the photoelectric properties of ZrS2/ZrSe2 heterostructures by using the first-principles calculation method based on density functional theory. The results show that the p-type doping is realized by B, C, and N atom doping, and the n-type doping is realized by O and F atom doping. The doping of B and C atoms produces impurity energy levels in the band gap, which affects the conductivity of the heterostructure. The band gap of N and O atom-doped heterostructures increases under tensile strain, but it is still a direct band gap. The analysis of the optical properties of the heterostructures shows that the doping of non-metallic atoms can adjust the optical absorption rate and reflectivity of the heterostructures. Under the action of tensile strain, the optical properties of the doped heterostructures have changed significantly in the low-energy region. This article provides a theoretical basis for the future application of ZrS2/ZrSe2 heterostructures. This paper uses the first-principles calculation method based on density functional theory. The PBE exchange-correlation functional based on generalized gradient approximation (GGA) is selected for the specific calculation, and the crystal structure is geometrically optimized by the ultrasoft pseudopotential method. It is verified that when the cutoff energy of the ZrS2/ZrSe2 heterostructure is 500 eV, the K-point grid is selected to be 10 × 10 × 2 with the lowest energy, so the cutoff energy is selected to be 500 eV. The K-point grid is selected to be 10 × 10 × 2. The convergence limits for structural optimization are as follows: the maximum force between atoms is 0.01 eV/Å, the convergence threshold of the maximum energy change is set to 10-9 eV/atom, and the convergence threshold of the maximum displacement is 0.001 Å. In order to avoid the influence of atomic periodic motion between different atomic layers, a vacuum layer of 20 Å is added in the vertical direction. Considering the interaction of vdW between the interfaces, the DFT-D2 method is used to verify. The optical properties were calculated by the random phase approximation method, and the K-point grid was selected as 12 × 12 × 2.

Exploring the Molecular-Scale Structures at Solid/Liquid Interfaces of Li-Ion Battery Materials: A Force Spectroscopy Analysis with Sparse Modeling.

In this study, we clarify the liquid structure formed at the interface between LiCoO2 (LCO), the cathode material of Li-ion batteries, and propylene carbonate (PC), which is used as a solvent in the electrolyte, on a molecular scale. We apply sparse modeling-based modal analysis to force spectroscopy data measured by frequency modulation atomic force microscopy (FM-AFM) and show that each component in the FM-AFM force curve, such as oscillatory solvation force, background, and noise, can be automatically decomposed. Moreover, by combining detailed force curve analysis with solid/liquid interface simulations based on first-principles calculation, we have identified that there are distinct damped vibrational modes in the force curves at the LCO/PC interface with a period of about 0.57 nm and those with shorter periods, which likely correspond to the solvation forces associated with bulk-state PC molecules and those with PC molecules in "lying down" orientations.

Lithium Polysulfide Catalytic Mechanism of AlN/InN Heterojunction by First-Principles Calculation

Lithium Polysulfide Catalytic Mechanism of AlN/InN Heterojunction by First-Principles Calculation

Strain Engineering of Ion-Coordinated Nanochannels in Nanocellulose.

Expanding the interlayer spacing plays a significant role in improving the conductivity of a cellulose-based conductor. However, it remains a challenge to regulate the cellulose nanochannel expanded by ion coordination. Herein, starting from multiscale mechanics, we proposed a strain engineering method to regulate the interlayer spacing of the cellulose nanochannels. First-principles calculations were conducted to select the most suitable ions for coordination. Large-scale molecular dynamics simulations were performed to reveal the mechanism of interlayer spacing expansion by the ion cross-linking. Combining the shear-lag model, we established the relationship between interfacial cross-link density and interlayer spacing of an ion-coordinated cellulose nanochannel. Consequently, fast ion transport and current regulation were realized via the strain engineering of nanochannels, which provides a promising strategy for the current regulation of a cellulose-based conductor.

Investigating the structural, electronic, and magnetic properties of Cd1-xVxTe: Insights from First-Principles calculations

In our study, we aim to investigate the structural, electronic, and magnetic properties of CdVTe, a diluted magnetic semiconductor with a Cadmium telluride (CdTe) structure. To achieve this, we employ the full potential linear augmented plane wave method implemented in the CASTEP code, which allows us to accurately simulate the behavior of the system. To describe the electronic exchange and correlation effects, we adopt the generalized gradient approximation (GGA) within the density functional theory (DFT) framework. This choice of methodology ensures reliable and accurate calculations of the electronic and magnetic properties of Cd1-xVxTe. Our results reveal that the substitution of Vanadium (V) into the Cadmium (Cd) site of the CdTe lattice does not alter the zinc blende crystal structure, demonstrating the stability and structural integrity of the system. Furthermore, our calculations demonstrate that the introduction of Vanadium impurities leads to spin polarization within the material, resulting in the emergence of a magnetic moment. Notably, we find that a Vanadium concentration of approximately x≈0.12 exhibits the strongest magnetic properties, characterized by a significant magnetic moment. These findings provide valuable insights into the behavior and potential applications of CdVTe as a diluted magnetic semiconductor, shedding light on the interplay between structural, electronic, and magnetic properties in this material.

Signature of excitonic insulators in phosphorene nanoribbons.

Phosphorene is a recently developed two-dimensional (2D) material that has attracted tremendous attention because of its unique anisotropic optical properties and quasi-one-dimensional (1D) excitons. We use first-principles calculations combined with the maximally localized Wannier function tight binding Hamiltonian (MLWF-TB) and Bethe-Salpeter equation (BSE) formalism to investigate quasiparticle effects of 2D and quasi-1D blue and black phosphorene nanoribbons. Our electronic structure calculations shows that both blue and black monolayered phases are semiconductors. On the other hand black phosphorene zigzag nanoribbons are metallic. Similar behavior is found for very thin blue phosphorene zig-zag and armchair nanoribbon. As a general behavior, the exciton binding energy decreases as the ribbon width increases, which highlights the importance of quantum confinement effects. The solution of the BSE shows that the blue phosphorene monolayer has an exciton binding energy four times higher than that of the black phosphorene counterpart. Furthermore, both monolayers show a different linear optical response with respect to light polarization, as black phosphorene is highly anisotropic. We find a similar, but less pronounced, optical anisotropy for blue phosphorene monolayer, caused exclusively by the quasi-particle effects. Finally, we show that some of the investigated nanoribbons show a spin-triplet excitonic insulator behavior, thus revealing exciting features of these nanoribbons and therefore provides important advances in the understanding of quasi-one dimensional phosphorus-based materials.&#xD.

Tuning Electronic and Structural Properties of Lead-Free Metal Halide Perovskites: A Comparative Study of 2D Ruddlesden-Popper and 3D Compositions.

In recent decades, two-dimensional (\num{2}D) perovskites have emerged as promising semiconductors for next-generation photovoltaics, showing notable advancements in solar energy conversion. Herein, we explore the impact of alternative inorganic lattice \ce{BX}-based compositions ($\ce{B} = \ce{Ge}$ or \ce{Sn}, $\ce{X} = \ce{Br}$ or \ce{I}) on the energy gap and stability. Our investigation encompasses \ce{BA$_2$MA$_{n-1}$B$_n$X$_{3n+1}$} \num{2}D Ruddlesden-Popper perovskites (for $n = \num{1} - \num{5}$ layers) and \num{3}D bulk \ce{(MA)BX$_3$} systems, employing first-principles calculations with spin-orbit coupling (SOC), DFT-1/2 quasiparticle, and D3 dispersion corrections. The study unveils how atoms with smaller ionic radii induce anisotropic internal and external distortions within the inorganic and organic lattices. Introducing the \ce{BA} spacers in the low-layer regime reduces local distortions but widens band gaps. Our calculation protocol provides deeper insights into the physics and chemistry underlying \num{2}D perovskite materials, paving the way for optimizing environmentally friendly alternatives that can efficiently replace \ce{Pb} with sustainable materials.

Anisotropic optical property of ferroelectric Bi2O2X (X=S,Se,Te) monolayer and strain engineering.

Based on the first-principles calculations, ferroelectric Bi2O2X (X=S, Se, Te) monolayers with unequivalent in-plane lattice constants are confirmed to be the ground state, which is consistent with the experiment result (Nano Lett. 19, 5703(2019)), and the anisotropic optical property is firstly investigated. We find that the polarizations of Bi2O2X monolayers points along the direction of a-axis, and Bi2O2T e monolayer process the largest polarization. Furthermore, both the biaxial and uniaxial strains are favor for the enhancement of polarization of Bi2O2X monolayers. It should be mentioned that the type of band gap will convert from indirect to direct for Bi2O2T e monolayer when the a-axial tensile strain is larger than 2%. At last, the optical absorption coefficient for Bi2O2X monolayers are calculated, and we obtain that Bi2O2Te monolayer has the strongest optical absorption within the range of visible light, the anisotropy and possible strain engineering to improve the optical absorption are discussed in detail. Our findings are significant in fields of optoelectronics and photovoltaics.&#xD.

Anisotropic optoelectronic properties of lead-free Cs3Bi2I9 single-crystal photodetectors

Abstract Lead-free halide perovskites have attracted widespread research interest due to their excellent optoelectronic properties and environmental friendliness. As one of the bismuth iodide compounds, Cs3Bi2I9 perovskite has been extensively explored in the field of photovoltaic devices and radiation detectors due to their non-toxic lead-free components and excellent stability. In this work, we successfully grow large-sized Cs3Bi2I9 single crystals (SC) with (l00) and (00l) crystal exposure facets by inverse temperature crystallization method. Under 525 nm light illumination with the intensity 15 mW/cm2 and 7 V bias, the (l00) Cs3Bi2I9 SC shows 2 times higher photocurrent, 3.6 times higher responsivity, and 2.8 times higher detectivity than the (00l) Cs3Bi2I9 SC, respectively. Superior response time in the scale of millisecond is obtained in both (l00) and (00l) Cs3Bi2I9 SCs. Based on the first-principle calculation, the (l00) SC possesses a higher charge distribution density and a wider dispersion distribution than (00l) SC, suggesting that more electrons in the (l00) SC can be excited in a wider range. The tighter arrangement of Cs+ and [BiI6]− octahedra in the (l00) SC than (00l) SC leads to the anisotropic photoelectric performance in (l00) and (00l) Cs3Bi2I9 SCs. Our results provide a strategy for the oriental growth of Cs3Bi2I9 SCs and the design of anisotropic optoelectronic devices with excellent performance.

Anatomy of Hidden Dzyaloshinskii-Moriya Interactions and Topological Spin Textures in Centrosymmetric Crystals.

The Dzyaloshinskii-Moriya interaction (DMI) is understood to be forbidden by the symmetry of centrosymmetric systems, thus restricting the candidate types for investigating many correlated physical phenomena. Here, we report the hidden DMI existing in centrosymmetric magnets driven by the local inversion symmetry breaking of specific spin sublattices. The opposite DMI spatially localized on the inverse spin sublattice favors the separated spin spiral with opposite chirality. Furthermore, we elucidate that hidden DMI widely exists in many potential candidates, from the first-principles calculations on the mature crystal database. Interestingly, novel topological spin configurations, such as the anti-chirality-locked merons and antiferromagnetic-ferromagnetic meron chains, are stabilized as a consequence of hidden DMI. Our understanding enables the effective control of DMI by symmetry operations at the atomic level and enlarges the range of currently useful magnets for topological magnetism.

Rational design of hollow C@SnS2/SnS@C cubes with N-doping to realize stepwise sodium storage induced structural stability for highly stable cycling performance

Current application of SnS2-based anodes in sodium-ion batteries fails to meet the demand for higher Na+ storage capacity due to the structural pulverization induced by significant volume changes during cycling. Herein, a novel composite of hollow C@SnS2/SnS@C cubes with excellent diffusion kinetics and structural stability from stepwise Na+ storage is designed. The SnS with gentle pseudocapacitive reaction acts as the self-supporter to prevent the collapse of SnS2/SnS heterostructure during the Na+ insertion/extraction. Together with the clamping of double carbon layers and hollow construction, the volume change caused structural instability in heterostructure is effectively alleviated. The rapid transport of electrons and ions within the built-in electric field of heterostructure endows the composite an excellent diffusion kinetics, which is generated from the narrow bandgap and outstanding Na+ absorption of heterostructure according to the first-principle calculation. Thus, the as-prepared composite maintains a remarkable capacity of 547.9 mAh g−1 after 200 cycles at current density of 0.5 A g−1. Even a capacity of 502.5 mAh g−1 with a low attenuation rate of 0.035 % per cycle is achieved at 2 A g−1 after 1000 cycles. The work proposes an effective approach to achieve effective Na + storage of SnS2-based anodes with long-lasting structural stability.

(C8H7N2O2)2[Bi2Br8]·2H2O and (C8H7N2O2)6[Bi2Cl10]Cl2·2H2O: Exploring Birefringent Crystals in Hybrid Halide Systems.

In this study, new hybrid birefringent crystals of (C8H7N2O2)2[Bi2Br8]·2H2O and (C8H7N2O2)6[Bi2Cl10]Cl2·2H2O were successfully synthesized by introducing a new birefringent group [C8H7N2O2]+ by a simple aqueous solution evaporation method. They crystallize in the P21/n space group, and their structure consists mainly of the π-conjugated group [C8H7N2O2]+ and the octahedron centered on Bi3+. By first-principles calculations, the birefringence response comes from the [C8H7N2O2]+ group with a planar π-conjugated structure. Meanwhile, the synthesis, structure, first-principles calculations, and optical properties are reported in this paper.

Ion solvation free energy calculations based on first-principles molecular dynamics thermodynamic integration.

Numerous electrochemistry reactions require the precise calculation of the ion solvation energy. Despite the significant progress in the first-principles calculations for crystals and defect formation energies for solids, the liquid system free energy calculations still face many challenges. Ion solvation free energies can be calculated via different semiempirical ways, e.g., using implicit solvent models or cluster of explicit molecule models; however, systematically improving these models is difficult due to their lack of a solid theoretical base. A theoretically sound approach for calculating the free energy is to use thermodynamic integration. Nevertheless, owing to the difficulties of self-consistent convergence in first-principles calculations for unphysical atomic configurations, the computational alchemy approach has not been widely used for first-principles calculations. This study proposes a general approach to use first-principles computational alchemy for calculating the ion solvation energy. This approach is also applicable for other small molecules. The calculated ion solvation free energies for Li+, Na+, K+, Be2+, Mg2+, and Ca2+ are close to the experimental results, and the standard deviation due to molecular dynamics fluctuations is within 0.06eV.

Long-Range Influence of Cr on the Stacking Fault Energy of Cr-Containing Concentrated Solid-Solution Alloys

Single-phase concentrated solid-solution alloys (CSAs) have exhibited excellent mechanical and radiation tolerance properties, making them potential candidate materials for nuclear applications. These excellent properties are closely related to dislocation movements, which depend on the stacking fault energies (SFEs). In CSAs, SFEs show large fluctuations due to variations in the local atomic environments in the vicinity of the stacking faults. In this work, first-principle calculations were performed to investigate the origin of the fluctuations in the SFEs of the widely studied CSA, NiCoCr, which show a very wide distribution from about −200 mJ/m2 to 60 mJ/m2. Compared to the common understanding that only atoms in close proximity to the stacking fault influence the SFEs in pure metals and dilute alloys, charge redistribution can be observed in several nearby planes of the stacking fault in NiCoCr, indicating that atoms several atomic layers away from stacking fault also contribute to the SFEs. Our analysis shows that Cr plays a major role in the large fluctuation in the SFEs of NiCoCr based on both electronic and magnetic responses. The flexible electronic structure of Cr facilitates easier charge transfer with Cr in several nearby atomic planes near the stacking fault, leading to significant changes in the d-electron number, orbital occupation number, and magnetic moments of Cr.

Anisotropic materials with abnormal Poisson's ratios and acoustic velocities.

Isotropic materials are required to adhere to various mechanical principles due to their limited thermal stability. For instance, it is essential for Poisson's ratio to be within the range of -1 to 0.5, and the longitudinal wave velocity must exceed the transverse wave velocity. Nevertheless, perfect crystals, as anisotropic materials, have the ability to defy conventional rules. Through the integration of high-throughput processes and first-principles calculations, a comprehensive exploration of known materials was conducted, resulting in the establishment of a database featuring an extreme anisotropic mechanism. This included the identification of abnormal Poisson's ratios (with the directional Poisson's ratio ranging from -3.00 to 3.67), the discovery of extreme negative linear compressibility, the determination of the upper and lower limits of the sound velocity, and other associated properties. Several materials with abnormal Poisson's ratios (<-1 or >0.5) were listed, and their peculiar mechanical behavior, wherein the volume decreased counterintuitively with uniaxial tension, was discussed. Finally, this study focused on the velocities of longitudinal and transverse waves, with specific emphasis on materials exhibiting transverse wave velocities that exceeded the longitudinal wave velocities.

Tunable electronic, transport, and optical properties of fluorine- and hydrogen-passivated two-dimensional Ga2O3 by uniaxial strain

Two-dimensional (2D) semiconductors have shown great prospects for future-oriented optoelectronic applications, whereas the applications of conventional 2D materials are significantly impeded by their low electron mobility (⩽200 cm2 V−1 s−1). In this work, strain-mediated fluorine- and hydrogen-passivated 2D Ga2O3 systems (FGa2O3H) have been explored via using first-principles calculations with the Heyd–Scuseria–Ernzerh and Perdew–Burke–Ernzerhof functionals. Our results reveal a considerable high electron mobility of FGa2O3H up to 4863.05 cm2 V−1 s−1 as the uniaxial tensile strain reaches 6%, which can be attributed to the enhanced overlapping of wave functions and bonding features. Overall, when applying uniaxial strain monotonously along the a(b) direction from compressive to tensile cases, the bandgaps of 2D FGa2O3H increase initially and then decrease, which originates from the changes of σ* anti-bonding in the conduction band minimum and π bonding states in the valence band maximum accompanying the lengthening Ga–O bonds. Additionally, when the tensile strain is larger than 8%, the stronger π bonding at the G point leads to an indirect-to-direct transition. Besides the highest electron mobility observed in n-type doped 2D FGa2O3H with 6% tensile strain, the electrical conductivity is enhanced and further elevated as the temperature increases from 300 K to 800 K. The variations of the absorption coefficient in the ultraviolet region are negligible with increasing tensile strain from 0% to 6%, which sheds light on its applications in high-power optoelectronic devices.