Crystal Orbitals Research Articles

Strain Engineering of ScYCCl2 MXene Monolayer and Intercalation of Metal-ions on MXene Surface: A DFT Study

The present work theoretically examines the influence of bi-axial strain on functionalized ScYCCl2 monolayer. The indirect to direct band gap transition on tensile conditions and the phase transition from semiconductor to metal on compressive strain have been studied. The analysis of extreme conditions of 10% compressive and 10% tensile strain through phonon and AIMD simulations underscore the kinetic and thermal stability, respectively of the monolayer under strain. These findings assure the possibility of experimental synthesis of the ScYCCl2 monolayer. After metallization, ScYCCl2 MXene is used as an anode of metal-ion (−Na, −K, −Li, −Mg) batteries as it has high theoretical storage capacity and low open circuit voltage. Work function engineering and the strain-dependent optical behavior of the ScYCCl2 monolayer have been examined. The work function of the ScYCCl2 monolayer has been raised under compressive strain and decreased under tensile strain. The Crystal Orbital Hamiltonian Population has been simulated under tensile and compressive strain to check the bond- strengths. Hence, the ScYCCl2 monolayer has the capability to alter its characteristics under strain. The improved optical characteristics recommend its applications in low-dimensional photonic devices and metallization after compressive strain recommends its energy storage applications in metal-ion batteries.

Optical Activities of Organic Crystals: Crystal Orbital Formulation of Anisotropic Gyration.

To describe the optical activities of crystals, gyration tensors are necessary for all wavelengths of incident light. To date, few studies on direct calculations of gyration tensors from Bloch states have been conducted. Herein, a practical procedure to obtain gyration tensors of organic crystals is presented using the crystal orbital method. The extended Hückel method was adopted to evaluate the gyration tensors, epitomizing the one-electron formulation. To confirm the validity of the formulation, the optical rotatory power of alanine and γ-glycine was examined. The reproduced profiles of the optical rotatory power were consistent with the results of recent experiments. This is a general formulation of the one-electron theory of optical activities for three-dimensional crystals. In principle, the optical rotatory strength tensor is not invariant with translation. For systems with small unit cells, however, the formalism is quasi-invariant with respect to translation. The quasi origin-independent formalism is sufficiently substantial to be applicable to modern crystal optics.

Coordination Engineering of Fe-Centered Catalysts for Superior Li-S Battery Performance.

Iron-nitrogen functionalized graphene has emerged as a promising cathode host for rechargeable lithium-sulfur batteries (RLSBs) due to its affordability and enhanced battery performance. To optimize its catalytical efficiency, we propose a novel approach involving coordination engineering. Our investigation spans a plethora of catalysts with varied coordination environments, focusing on elements B, C, N and O. We revealed that Fe-C4 and Fe-B2C2-h are particularly effective for promoting Li2S oxidation, whereas Fe-N4 excels in catalyzing the sulfur reduction reaction (SRR). Importantly, our study identified specific descriptors - namely, the Integrated Crystal Orbital Hamilton Population (ICOHP) and the bond length between Fe and S in Li2S adsorbed state - as the most effective predictive descriptors for Li2S oxidation barriers. Meanwhile, Li2S adsorption energy emerges as a reliable descriptor for assessing the SRR barrier. These identified descriptors are expected to be instrumental in rapidly identifying promising cathode hosts across various metal-centered systems with diverse coordination environments. Our findings not only offer valuable insights into the role of coordination environment, but also present an effective path for rapidly identifying high performance catalysts for RLSBs, enabling the acceleration of advanced RLSBs development.

Computational insights into the structure and decomposition behaviors of 2,4,6-triamino-5-nitropyrimidine-1,3-dioxide under high pressure up to 10 GPa.

Inspired by the recent successful synthesis of the energetic compound 2,4,6-triamino-5-nitropyrimidine-1,3-dioxide (ICM-102), which displayed a good balance between high energy and sensitivity, the response of the structure and decomposition behaviors of ICM-102 to high pressure was systematically investigated using first principle calculations. ICM-102 exhibited a graphite-like layer structure, with the c-axis and the a-axis mainly contributing to the distance between the molecular planes. As the pressure increased from 1atm to 10 GPa, this distance decreased from 3.166 to 2.689 Ǻ. The hydrogen bonds had the most contribution to the non-covalent interactions within the same molecular planes, resulting in the b-axis discontinuity. However, van der Waals interactions gradually appeared between molecular planes as the pressure increased to 2.5 GPa. Based on the analysis of crystal orbitals, the distribution of π bonds and the Laplacian bond order (LBO), it was determined that the generation mechanism of H2O molecules involved the cleavage of N-Oc (coordinated oxygen atoms), followed by intermolecular hydrogen transfer reactions, and ultimately the formation of H2O molecules through competition with H atoms in the amino groups within the same molecules. More importantly, the pressure dependence of LBO values for N-Oc revealed that high pressure could inhibit the ICM-102 decomposition process due to reinforcing hydrogen bonds and van der Waals interactions. This work will deepen our understanding of the stability of ICM-102 under high pressure and provide a helpful reference for its potential detonation applications. All simulations, including geometry optimization and vibration analysis under quasi-hydrostatic pressure, were conducted using the CP2K code. The PBE function and the Goldk-Teter-Hutter (GTH) pseudopotential with the double-ζ-with-polarization (DZVP) basis set were employed. Additionally, the semiempirical dispersion correction D3 (BJ) was used to account for the intermolecular dispersion force. The simulations were performed under periodic boundary conditions, with a finest grid level cutoff set to 400 Ry for the Γ point. The Broyden-Flecher-Goldfarb-Shanno (BFGS) optimization method was used, with tighter convergence criteria applied for the subsequent calculations of infrared spectra. Finally, the wave-function analysis, such as non-covalent interaction and LBO, was conducted using the Multiwfn and VMD packages.

Ab initio simulation of the dynamic shock response of single crystal and lightweight multicomponent alloy

The dynamic response of shock wave impact on single crystal aluminium and lightweight multicomponent alloy Al-Cu-Li-Mg is simulated by using the combination of Ab initio Molecular Dynamics (AIMD) and Multi-Scale Shock Technique (MSST), with the analysis carried out at the atomic/electronic levels. The simulation is verified by comparing the particle velocity of single crystal obtained in this work with the data in literature. The shock compression process not only involves the migration of atoms, but also is related to electronic transition. Two stages could be found in the shock compression process: oscillatory compression of the crystal cell and oscillatory migration of the atoms. The crystal structure of the multicomponent alloy could be disordered even at low shock speed, due to the difference in the ability to migrate between different kinds of atoms. As the sample is shock-compressed, the contribution proportion of crystal orbitals shows a sharp decrease for D orbital, while it increases significantly for S orbital and P orbital. The electron structure shows a quicker response to the shock wave compression process than the crystal structure. The orbital contribution from P orbital of the crystal is mainly due to the P orbital of Al atoms, while the orbital contribution from D orbital of the crystal is mainly due to the D orbital of Cu atoms. Total Density of States (TDOS) is mainly contributed by the Projected Density of State (PDOS) of Cu atoms in the occupied state of energy levels, while it is close to the PDOS of Al atoms in the non-occupied state of energy levels.

Exploring the key role of electronic effects in Pd catalysts supported on Ni-modified N-doped porous carbon for direct synthesis of H2O2 under atmospheric pressure

Developing high-performance catalysts for direct synthesis of hydrogen peroxide (H2O2; DSHP) remains challenging. The paper proposes a strategy to optimize Pd catalysts by enriching Ni-modified N-doped porous carbon. The theoretical studies indicate that Ni modulates the electronic structure of Pd through charge transfer and strain effects, effectively reducing the adsorption energy of reaction intermediates. Furthermore, the Crystal Orbital Hamilton Population (COHP) analysis demonstrated that the energy-integrated COHP (ICOHP) of the OO bonds in O2* (* denotes the adsorbed state), OOH*, and HOOH* were positively correlated with dissociation activation energies and negatively correlated with hydrogenation activation energies. The modulation of the electronic structure of Pd by Ni reduces the hydrogenation activation energies of O2* and OOH* and increases the dissociation activation energies of O2*, OOH*, and HOOH*, thereby enhancing the selectivity and productivity of H2O2. The experimental results aligned with the theoretical predictions, showcasing that the optimized Pd-Ni7.5/NPCs catalyst exhibited remarkable H2O2 selectivity and productivity of 70.5 % and 230.31 mol kgcat-1·h-1, respectively, surpassing the performance of Pd/NPCs (47.1 % and 52.3 mol kgcat-1·h-1). This study offers profound insights into the crucial role of electron effects in DSHP using nickel-modified palladium catalysts, providing strong evidence into the design and optimization of Pd catalysts.

Exploring the Nature of Ag–Ag Interactions in Different Tellurides by Means of the Crystal Orbital Bond Index (COBI)

Exploring the Nature of Ag–Ag Interactions in Different Tellurides by Means of the Crystal Orbital Bond Index (COBI)

A comparative study of structural, mechanical, electronic and optical properties of InTe monolayer & homo-bilayer

We have done a comparative study between InTe monolayer and van der Waals homo-bilayer. Four different stacking AA, AAʹ, AB, ABʹ are taken into consideration in which possible three stackings AAʹ, AB, ABʹ are considered for the calculations. The homo-bilayer is simulated for the 10ps via ab-intio molecular dynamics to ensure kinetic stability. The mechanical properties show that the InTe homo-bilayer is more elastic than its monolayer. The indirect to direct bandgap modulations are seen from the electronic properties. The Telluride atom has a major role in the electron transfer between the layer of the homo-bilayer. The Crystal Orbital Hamilton Populations calculation is done to observe the bonding nature. Enhancement is seen in the absorption coefficient of order 105 cm−1. The highest refractive index of 3.88 is observed in the visible region. Reflectivity increases in the visible and ultraviolet regions whereas the transmittivity decreases in homo-bilayer than monolayer.

Identification of Cu(111) as Superior Active Sites for Electrocatalytic NO Reduction to NH3 with High Single‐Pass Conversion Efficiency

AbstractOpting for NO as an N source in electrocatalytic NH3 synthesis presents an intriguing approach to tackle energy and environmental challenges. However, blindly pursuing high NH3 synthesis rates and Faradaic efficiency (FE) while ignoring the NO conversion ratio could result in environmental problems. Herein, Cu nanosheets with exposed (111) surface is fabricated and exhibit a NO‐to‐NH3 yield rate of 371.89 μmol cm−2 h−1 (flow cell) and the highest FE of 93.19±1.99 % (H‐type cell). The NO conversion ratio is increased to the current highest value of 63.74 % combined with the development of the flow cell. Additionally, Crystal Orbital Hamilton Population (COHP) clearly reveals that the “σ‐π* acceptance‐donation” is the essence of the interaction between the Cu and NO as also supported by operando attenuated total reflection infrared spectroscopy (ATR‐IRAS) in observing the key intermediate of NO−. This work not only achieves a milestone NO conversion ratio for electrocatalytic NO‐to‐NH3, but also proposes a new descriptor that utilizes orbital hybridization between molecules and metal centers to accurately identify the real active sites of catalysts.

Identification of Cu(111) as Superior Active Sites for Electrocatalytic NO Reduction to NH3 with High Single-Pass Conversion Efficiency.

Opting for NO as an N source in electrocatalytic NH3 synthesis presents an intriguing approach to tackle energy and environmental challenges. However, blindly pursuing high NH3 synthesis rates and Faradaic efficiency (FE) while ignoring the NO conversion ratio could result in environmental problems. Herein, Cu nanosheets with exposed (111) surface is fabricated and exhibit a NO-to-NH3 yield rate of 371.89 μmol cm-2 h-1 (flow cell) and the highest FE of 93.19±1.99 % (H-type cell). The NO conversion ratio is increased to the current highest value of 63.74 % combined with the development of the flow cell. Additionally, Crystal Orbital Hamilton Population (COHP) clearly reveals that the "σ-π* acceptance-donation" is the essence of the interaction between the Cu and NO as also supported by operando attenuated total reflection infrared spectroscopy (ATR-IRAS) in observing the key intermediate of NO- . This work not only achieves a milestone NO conversion ratio for electrocatalytic NO-to-NH3 , but also proposes a new descriptor that utilizes orbital hybridization between molecules and metal centers to accurately identify the real active sites of catalysts.

Positive and Negative Contribution from Lead–Oxygen Groups and Halogen Atoms to Birefringence: A First Principles Investigation

Oxyhalides, containing oxygen and halogen atoms and combining the advantages of oxides and halides in geometry and optical response, have great potential in optical materials. In this study, the electronic structures and the optical properties of the Pb3O2X2 (X = Cl, Br, I) compounds have been investigated using the first principles method. The results show that these compounds have birefringence at 0.076, 0.078, and 0.059 @ 1064 nm, respectively. And, the asymmetric stereochemical active lone pair electrons were found around lead atoms, which were confirmed by the projected density of states, the electronic localization functions, and the crystal orbitals. The contribution of atoms and polyhedra to birefringence was further evaluated using the Born effective charge. The results show that halogen atoms give negative contribution, and lead—oxygen polyhedra give positive contribution. The spin—orbit coupling effect is also investigated, and the downshift of the conduction band and variation in the valence band are found after relevant spin—orbit coupling (SOC), which leads to a reduction in the band gap and birefringence.

The Zintl Concept Applied to Intergrowth Structures: Electron‐Hole Matching, Stacking Preferences, and Chemical Pressures in Pd5InAs

AbstractEnumerating the potential stacking sequences of layers is a fundamental way to account for the structure diversity of solid state compounds. In many cases, these stacking variations represent polymorphs with only small energetic differences. Here, we examine a compound for which the preferred stacking pattern instead reveals key aspects about its chemical bonding: Pd5InAs. Its structure is based on the intergrowth of slabs of the AuCu3 and PtHg2 (or alternatively, fluorite) structure types. Two basic stacking arrangements are available to this compound, represented by the Pd5TlAs and HoCoGa5 structure types. DFT total energy calculations reveal that the former outcompetes the latter by a staggering 0.65 eV/formula unit. Through a combination of DFT‐reversed approximation Molecular Orbital (DFT‐raMO) and DFT‐Chemical Pressure (DFT‐CP) analysis we trace this preference to two factors. First, with DFT‐raMO analysis, we derive a Zintl‐like bonding scheme of Pd5InAs. This scheme, along with the inspection of selected crystal orbitals, is then connected to preferred stacking through the coordination environments of the Pd atoms at the interface between the Pd−In and Pd−As layers. In the hypothetical HoCoGa5‐type and observed Pd5InAs‐type structures, different Pd coordination environments arise at the interfaces. The hypothetical structure features square planar PdIn2As2 units, in each of which the same 4d‐orbital serves in the Pd sublattice's role as both Lewis acid (for interactions with the As) and Lewis base (for interactions with the In). In the observed structure, tetrahedral PdIn2As2 units occur instead, so that these contradictory roles are distributed to separate 4d‐orbitals, leading to more effective bonding. DFT‐CP analysis illustrates that this driving force for the Pd5TlAs‐type arrangement is supplemented by a favorable alignment of the packing tensions in the parent structures. Altogether, the resulting picture demonstrates how the reaction of simple intermetallic structures to form intergrowths can be guided by recognizable chemical interactions.

Phase transitions in transition-metal dichalcogenides with strain: insights from first-principles calculations

It is well known that monolayer transition-metal dichalcogenides (MX2, M = Mo, W and X = S, Se, Te) could exist in three common structures, i.e. 1T, , and 1H phases. In order to reveal their possible phase transitions driven by biaxial strain, we use first-principles calculations to determine the energy landscapes associated with these three phases. Due to its intrinsic dynamical instability, the centrosymmetric 1T phase is known to be metastable and will transform into the phase where pairs of metal atoms pull together toward each other. Moreover, controlling the metallic and semiconducting 1H phases is of particular interest as this can introduce novel opportunities in a series of applications. When a biaxial strain is simultaneously compressed along the zigzag direction and stretched along the armchair direction, phase transitions from 1H to have occurred in MSe2 and MoTe2, but for MSe2 the biaxial strain is much difficult to realize in experiments. For WTe2, the structure will transform into the 1H form when a biaxial strain is just compressed along the armchair direction. The transitions between 1H and phases could be attributed to the changes of metal-chalcogen bonds along the armchair direction by analyzing the Crystal Orbital Hamilton Population. Only half M-X bonds along the armchair direction is the main factor, because their lengths are robust in phase and decrease in 1H form with the tensile strain applied along the armchair direction. Our findings provide a guideline for the phase engineering of transition-metal dichalcogenides with biaxial strain.

Impact of d-states on transition metal impurity diffusion in TiN

In this work, we studied the energetics of diffusion-related quantities of transition-metal impurities in TiN, a prototype ceramic protective coating. We use ab-initio calculations to construct a database of impurity formation energies, vacancy-impurity binding energies, migration, and activation energies of 3d and selected 4d and 5d elements for the vacancy-mediated diffusion process. The obtained trends suggest that the trends in migration and activation energies are not fully anti-correlated with the size of the migration atom. We argue that this is caused by a strong impact of chemistry in terms of binding. We quantified this effect for selected cases using the density of electronic states, Crystal Orbital Hamiltonian Population analysis, and charge density analysis. Our results show that the bonding of impurities in the initial state of a diffusion jump (equilibrium lattice position), as well as the charge directionality at the transition state (energy maximum along the diffusion jump pathway), significantly impact the activation energies.

Computational Screening of Two-Dimensional Metal-Organic Frameworks as Efficient Single-Atom Catalysts for Oxygen Reduction Reaction.

The development of efficient and inexpensive oxygen reduction reaction (ORR) catalysts is crucial for renewable energy technologies. Herein, using density functional theory (DFT) methods and microkinetic simulations, we systematically investigated the ORR catalytic performance of a series of 2D metal-organic frameworks M3 (HADQ)2 (HADQ=2,3,6,7,10,11-hexaamine dipyrazino quinoxaline). It shows that all 2D M3 (HADQ)2 (M=Cr, Mn, Fe, Co, Ni, Cu, Ru, Rh and Pd) monolayers are metallic, due to π-conjugated crystal orbitals centered on the central metals and ligand N atoms. The catalytic activity of M3 (HADQ)2 depends on the binding strength between ORR intermediates and metal species, and can be tuned via changing the central metals. Among these candidates, Rh3 (HADQ)2 and Co3 (HADQ)2 show superior ORR performance to Pt (111) with high half-wave potentials of 0.99 and 0.93 V, respectively. Moreover, the screened two catalysts have excellent intermediate-tolerance ability for dynamic coverage of oxygenated species on the active sites. Our work provides a new path towards developing efficient ORR electrocatalysts.

Comprehensive Mechanism for CO Electroreduction on Dual-Atom-Catalyst-Anchored N-Doped Graphene.

Carbon neutrality has drawn increasing attention for realizing the carbon cyclization and reducing the greenhouse effect. Although the C1 products, such as CO, can be achieved with a high Faraday efficiency, the targeted production of C2 fuels as well as the mechanism have not been systematically investigated. In this work, we carry out a first-principles study to screen dual-atom catalysts (DACs) for producing C2 fuels through the electrocatalytic carbon monoxide reduction reaction (e-CORR). We find that methanol, ethanol and ethylene can be produced on both DAC-Co and DAC-Cu, while acetate can be achieved on DAC-Cu only. Importantly, methanol and ethylene are preferred on DAC-Co, while acetate and ethylene on DAC-Cu. Furthermore, we show that the explicit solvent can enhance the adsorption and influence the protonation steps, which subsequently affects the protonation and dimerization behavior as well as the performance and selectivity of e-CORR on DACs. We further demonstrate that the C-C coupling is easy to be formed and stabilized if the Integrated Crystal Orbital Hamilton Population (ICOHP) is low because of the low energy barrier. Our findings provide not only guidance on the design of novel catalysts for e-CORR, but an insightful understanding on the reduction mechanism.

Enhanced magnetism in Ru‐doped hybrid improper perovskite Ca3Mn2O7 via experimental and first‐principles study

AbstractThe quasi‐2D Ruddlesden–Popper layered perovskites with single‐phase multiferroicity have attracted much attention in recent magnetoelectric memory applications. We have systematically investigated the optical and magnetic properties of Ru‐doped hybrid improper perovskite Ca3Mn2O7 both experimentally and using first‐principles calculations. The Ca3RuxMn2−xO7 (x = 0, 0.02, 0.06, 0.10) powders were successfully synthesized using the sol–gel method. Ru doping enhanced the ferromagnetism of Ca3Mn2O7. The quasi‐2D antiferromagnetic (AFM) fluctuation effect was observed in Ca3Mn2O7 and in certain doped samples. We also found that the optical bandgaps of the doped samples were reduced after doping, agreeing with the results of the first‐principles calculations. Infrared absorption spectra indicated the distortion of the Mn–O bonds was possibly due to the Jahn–Teller effect. To analyze the electronic structures and to understand the detailed atomic contributions of magnetic moments, the Crystal Orbital Hamilton Population (COHP) and Integrated COHP of the (Mn,Ru)O6 octahedra were also calculated. The Ru‐doped materials with the coexistence of ferromagnetic and AFM orderings are expected to be excellent candidates for magnetoelectric devices.

Tuning Structural and Optical Properties of Porphyrin-based Hydrogen-Bonded Organic Frameworks by Metal Insertion.

Herein, a simple way of tuning the optical and structural properties of porphyrin-based hydrogen-bonded organic frameworks (HOFs) is reported. By inserting transition metal ions into the porphyrin cores of GTUB-5 (p-H8 -TPPA (5,10,15,20-Tetrakis[p-phenylphosphonic acid] HOF), the authors show that it is possible to generate HOFs with different band gaps, photoluminescence (PL) life times, and textural properties. The band gaps of the resulting HOFs (viz., Cu-, Ni-, Pd-, and Zn-GTUB-5) are measured by diffuse reflectance and PL spectroscopy, as well as calculated via DFT, and the PL lifetimes are measured. Across the series, the band gaps vary over a narrow range from 1.37 to 1.62 eV, while the PL lifetimes vary over a wide range from 2.3 to 83 ns. These differences ultimately arise from metal-induced structural changes, viz., changes in the metal-to-nitrogen distances, number of hydrogen bonds, and pore volumes. DFT reveals that the band gaps of Cu-, Zn-, and Pd- GTUB-5 are governed by highest occupied/lowest unoccupied crystal orbitals (HOCO/LUCO) composed of π- orbitals on the porphyrin linkers, while that of Ni-GTUB-5 is governed by a HOCO and LUCO composed of Ni dorbitals. Overall, our findings show that metal-insertion can be used to optimize HOFs for optoelectronics and small-molecule capture applications.

Quantum anomalous Hall effect in M2X3 honeycomb Kagome lattice

The quantum anomalous Hall (QAH) effect has recently drawn great attention in spintronics with extraordinary property of chiral edge states without dissipation in absence of magnetic field. In M2X3 honeycomb Kagome lattice, numerous two-dimensional materials are predicted to be QAH insulators including metal oxides/sulfides and metal organic lattice. In this work, we proposed a general model to explain the mechanism of Dirac half metal with absence of spin orbital coupling and the nontrivial topological property with spin orbital coupling, which could be induced by combination of electron counting rule, crystal field effect and orbitals hybridization. Based on the mechanism, we further predict that triphenyl-metal lattice M2(C6H4)3 (M = V, Nb, Ta) are all QAH insulators with high Curie temperature and large nontrivial band gap for triphenyl-Nb and triphenyl-Ta lattice.

Crystal structure and electronic properties of BrF under high-pressure

The bromine monofluoride (BrF)’s structural and electronic properties are widely investigated under high pressure through ab initio computations using the density functional theory. At first, five new high-pressure BrF structures are uncovered with space groups P21, Cmc21, Pmc21, P2/m, and R-3m symmetry under pressure. According to the elastic constants and phonon scattering curves, the mechanical and dynamical stability of Cmc21, Pmc21, P2/m, and R-3m structures under pressures can be demonstrated. The band structure and density of states show that the P21 and the Cmc21 structures are non-metallic. In contrast, the Pmc21, P2/m, and R-3m structures are metallic. Band structure computations indicate that the Cmc21 structure is a semiconductor containing a direct band gap (Eg) of 0.90 eV (PBE) or 2.19 eV (HSE06). The electron localization function (ELF), Crystal Orbital Hamilton Population (COHP), and charge density difference calculations indicate that the Cmc21 structure is covalent. Besides, we find that the BrF (at ambient pressure) transition temperature of 144 K (∼−129 °C) is compatible with the previous experiment that reported that the temperature when BrF precipitates into a fine yellow powder is below −130 °C. The results enrich the BrF's high-pressure phases and establish the basis for subsequent theoretical and experimental studies.