Electron Density Functional Theory Research Articles

Designing a 3D framework NaVOPO4 as a high-power, low-strain and long-life positive electrode material for Na-ion batteries

Advances in sodium-ion batteries hugely rely on perfecting the performance of active electrode materials. In this paper, we offer a new NaVOPO4 polymorph adopting a KTiOPO4-type framework as a promising high-rate, low-strain and long-life positive electrode material for sodium-ion batteries. NaVOPO4 is prepared via a facile hydrothermally-assisted solid-state ion exchange. The crystal structure is refined based on synchrotron X-ray powder diffraction (XRD) and validated by X-ray absorption spectroscopy (XAS), transmission electron microscopy analysis and density functional theory (DFT) calculations. The electrochemical performance of NaVOPO4 is evaluated through galvanostatic charge/discharge tests, cyclic voltammetry, potentiostatic intermittent titration and electrochemical impedance spectroscopy. A carbon-coated NaVOPO4 demonstrates a specific capacity of ∼110 mAh g–1 at a C/10 rate at an average potential of ∼3.93 V vs. Na+/Na. The material exhibits decent capacity retention and rate capability, maintaining over 74 % of the initial capacity after 1000 cycles at a C/2 rate and around 87% at a 2C charge/discharge rate. Diffusion coefficients of 10–11 cm2 s–1 and DFT-calculated energy barriers of ∼0.15–0.35 eV indicate fast Na+ ion diffusion within the NaVOPO4 framework. The vanadium oxidation state and charge compensation mechanism in NaVOPO4 are studied by XAS and DFT. Moreover, the operando XRD analysis coupled with DFT elucidates several phase transitions occurring in NaVOPO4 with an exceptionally low volume variation of 2.4% in total, highlighting the reversibility and structural stability during Na+ de/insertion. Overall, NaVOPO4 exhibits attractive electrochemical performance and stability, making it a potential candidate for sodium-ion battery cathode materials.

Quantum symmetrization transition in superconducting sulfur hydride from quantum Monte Carlo and path integral molecular dynamics

AbstractWe study the structural phase transition, originally associated with the highest superconducting critical temperature Tc measured in high-pressure sulfur hydride. A quantitative description of its pressure dependence has been elusive for any ab initio theory attempted so far, raising questions on the actual mechanism leading to the maximum of Tc. Here, we estimate the critical pressure of the hydrogen bond symmetrization in the Im$$\bar{3}$$ 3 ¯ m structure, by combining density functional theory and quantum Monte Carlo simulations for electrons with path integral molecular dynamics for quantum nuclei. We find that the Tc maximum corresponds to pressures where local dipole moments dynamically form on the hydrogen sites, as precursors of the ferroelectric Im$$\bar{3}$$ 3 ¯ m-R3m transition, happening at lower pressures. For comparison, we also apply the self-consistent harmonic approximation, whose ferroelectric critical pressure lies in between the ferroelectric transition estimated by path integral molecular dynamics and the local dipole formation. Nuclear quantum effects play a major role in a significant reduction (≈50 GPa) of the classical ferroelectric transition pressure at 200 K and in a large isotope shift (≈25 GPa) upon hydrogen-to-deuterium substitution of the local dipole formation pressure, in agreement with the corresponding change in the Tc maximum location.

Nanoscale Localized Phonons at Al2O3 Grain Boundaries.

Nanoscale defects like grain boundaries (GBs) would introduce local phonon modes and affect the bulk materials' thermal, electrical, optical, and mechanical properties. It is highly desirable to correlate the phonon modes and atomic arrangements for individual defects to precisely understand the structure-property relation. Here we investigated the localized phonon modes of Al2O3 GBs by combination of the vibrational electron energy loss spectroscopy (EELS) in scanning transmission electron microscope and density functional perturbation theory (DFPT). The differences between GB and bulk obtained from the vibrational EELS show that the GB exhibited more active vibration at the energy range of <50 meV and >80 meV, and further DFPT results proved the wide distribution of bond lengths at GB are the main factor for the emergence of local phonon modes. This research provides insights into the phonon-defect relation and would be of importance in the design and application of polycrystalline materials.

Atomic structure of clusters and GP-zones in an Al-Mg-Si alloy

The structures of atomic clusters and GP-zones preceding precipitation in an Al-0.70Mg-0.85Si-0.15Fe-0.25Mn (wt. %) alloy have been investigated by annular dark field scanning transmission electron microscopy imaging and density functional theory calculations. One analysed condition consisted of one-month natural aging, and another was in a pre-baked state after 24 h aging at 90 °C. To quantify the effect of clusters and GP-zones on microstructure development during a subsequent isothermal artificial aging at 185 °C, hardness evolution was measured, and precipitate statistics was performed at the peak hardness corresponding to 5 h aging. Several types of clusters and GP-zones have been observed and based on structure and visual appearance named 'Disordered Frank-Kasper', 'Square', 'Binocular', and 1β", 2β", 3β", 4β". The first four are most common in the natural aged condition and are associated with a well-known influence on the initial hardness evolution. The last three are fractions of the hardening β" phase and refer to the number of structural units they contain. These GP-zones have an increased density in the pre-baked condition, leading to a rapid increase in hardness during artificial aging and to a finer microstructure of high-density small needle precipitates at peak hardness.

New Azo Derivative of β-Diketones and Its Cu(II), Co(II) Complexes: Synthesis, Theoretical Study and Biological Activity.

The paper is focused on biological activity and theoretical study of the structure and properties of a new azo derivative of β-diketones and its complexes with some metals. The aim of our work was to study the structure and properties of the newly synthesized compound as well as to theoretically determine the possibility of complex formation with the Cu(II) or Co(II) ions. A compound with the same substituents R1=R2=CH3 was chosen for the study. A synthesized azo compound based on 4-amino antipyrine and its complexes with Cu(II), Co(II) in solution and solid phase is reported. The structures of these compounds have been testified by X-ray, IR and NMR spectroscopy. The combined experimental and theoretical approach was used. To study the structure and properties of the synthesized compound, as well as its possible complex formation with the Cu(II), quantum-chemical calculations were carried out the 6-31G basis set and the electron density functional theory (DFT) method. These 3-(1-phenyl-2,3-dimethyl-pyrazolone-5) azopentadione-2,4 (PDPA) with Cu(II) and Co(II) complexes had effective inhibition against butyrylcholinesterase and acetylcholinesterase. IC50 values were found as 19.03, 3.64 μM for AChE and 28.47, 8.01 μM for BChE, respectively. Cholinesterase inhibitors work to slow down the acetylcholine's deterioration.

Structural study of nickelate based heterostructures

Heterostructures consisting of SmNiO3 and NdNiO3 alternating layers with additional LaAlO3 spacer layers were grown and fully characterized by means of x-ray diffraction, atomic force microscopy, and scanning transmission electron microscopy. A change in the orientation of the orthorhombic long-axis of the nickelate layers is observed when a single unit cell of LaAlO3 is inserted between SmNiO3 and NdNiO3, in agreement with density functional theory calculations. At the same time, the structure of the ultra-thin rhombohedral LaAlO3 layers is affected by their proximity to orthorhombic nickelate layers, with both scanning transmission electron microscopy studies and density functional theory calculations revealing a weak antipolar motion of the La-cation in the LaAlO3 layers that is not present in the bulk rhombohedral structure of this compound.

In Situ Modulation of Oxygen Vacancies on 2D Metal Hydroxide Organic Frameworks for High-Efficiency Oxygen Evolution Reaction.

The discovery of non-precious catalysts for replacing the precious metal of ruthenium in the oxygen evolution reaction (OER) represents a key step in reducing the cost of green hydrogen production. The 2D d-MHOFs, a new 2D materials with controllable oxygen vacancies formed by controlling the degree of coordination bridging between metal hydroxyl oxide and BDC ligands are synthesized at room temperature, exhibit excellent OER properties with low overpotentials of 207 mV at 10mA cm-2. High-resolution transmission electron microscopy images and density functional theory calculations demonstrate that the introduction of oxygen vacancy sites leads to a lattice distortion and charge redistribution in the catalysts, enhancing the OER activity of 2D d-MHOFs comprehensively. Synchrotron radiation and in situ Raman/Fourier transform infrared spectroscopy indicate that part of oxygen defect sites on the surface of 2D d-MHOFs are prone to transition to highly active metal hydroxyl oxides during the OER process. This work provides a mild strategy for scalable preparation of 2D d-MHOFs nanosheets with controllable oxygen defects, reveals the relationship between oxygen vacancies and OER performance, and offers a profound insight into the basic process of structural transformation in the OER process.

Identification and quantification of electron-generated atomic hydrogen through in-situ electron spin resonance and density functional theory

Electro-generated atomic hydrogen (H*) plays a crucial role in the electrochemical reduction process, serving as the key species that mediates the electrocatalytic hydrogenation reduction of stubborn contaminants in wastewater treatment. However, precisely identifying and quantifying transient atomic H* presents a significant challenge due to its limited lifespan and its existence solely within the boundary layer at the electrode/solution interface. Herein, we developed the electrodeposition of palladium nanoparticles onto carbon cloth and assessed its effectiveness as a cathode for generating and stabilizing atomic H*. The environmental application of atomic H* was validated through the dechlorination of 2, 4-dichlorophenol wastewater and the reduction of antimonite wastewater. Additionally, the identification of atomic H* was verified by electrochemical measurements, high-resolution mass spectra, and density functional theory. Moreover, introducing an excess of a spin trapping agent (5,5-dimethyl-1-pyrroline-N-oxide) and fast in-situ spin trapping facilitated creating favorable conditions for efficient trapping of atomic H* and subsequent electron spin resonance (ESR) spectroscopy quantification analysis. Subsequently, the quantification of atomic H* was achieved by double integration of the ESR signal of spin adduct and comparison with the external standard agent (4-hydroxy-2,2,6,6-tetramethyl-1-piperidine 1-oxyl). This study introduces a novel method for in-situ spin trapping and quantification of atomic H*, facilitating the advancement of electrochemical reduction technology and its application in wastewater treatment.

Novel complex compounds of nickel with 3-(1-phenyl-2,3-dimethyl-pyrazolone-5)azopentadione-2,4: synthesis, NBO analysis, reactivity descriptors and in silico and in vitro anti-cancer and bioactivity studies

A synthesized azo compound based on 4-amino antipyrine and its complexes with Ni(II) in solution and solid phase is reported. The structures of these compounds have been testified by IR and NMR spectroscopy. The combined experimental and theoretical approach was used. To study the structure and properties of the synthesized compound, as well as its possible complex formation with the Ni(II), ab initio quantum-chemical calculations were carried out using the Hartree-Fock (HF) method with the 6-31 G basis set and the electron density functional theory (DFT) method with hybrid three-parameter potential B3LYP and extended basis set 6-311++G(d,p) taking into account polarization and diffuse functions for all atoms. The geometric, energy, and electronic parameters were calculated and analyzed. The HOMO-LUMO energy gap has been calculated to determine chemical activity. Both complexes had effective inhibition against butyrylcholinesterase and acetylcholinesterase. IC50 values were found as 19.43 and 27.08 µM for AChE, 2.37 and 7.40 µM for BChE, respectively. For the anticancer outcome, high doses of compound E1 inhibited viability by about 40–45%, while this rate was around 65–70% for compound E2 at the same doses. Anticholinesterase and anticancer potential of compounds E1 and E2 also evaluated by in silico techniques. Both compounds show strong binding to VEGFR1, with E2 exhibiting superior inhibitory activity in hAChE and hBChE through shorter and stronger interactions. MD simulations suggest that E2 forms more stable complexes with hAChE and hBChE compared to E1, making it a promising candidate for further exploration in anticancer and anticholinesterase therapies. Communicated by Ramaswamy H. Sarma

Segregation-induced strength anomalies in complex single-crystalline superalloys

Pushing the maximum service temperature of aircraft engines and industrial gas turbines is the major pathway to improve their energy efficiency and reduce CO2 emissions. This maximum is mostly limited by the temperature capability of key-component materials, including superalloys. In this alloy class, segregation of elements facilitates plastic deformation and is generally considered to cause softening during high-temperature deformation. Here, we show that segregation-assisted processes can also lead to strengthening and induce an anomalous increase of the yield strength. Atomic-resolution transmission electron microscopy and density functional theory calculations reveal a segregation-assisted dissociation process of dislocations at precipitate-matrix interfaces in combination with atomic-scale reordering processes. These processes lead to an inhibition of athermal deformation mechanisms and a transition to stacking fault shearing, which causes the strengthening effect. Unraveling these elementary mechanisms might guide a mechanism-based alloy design of future superalloys with enhanced high-temperature capabilities.

Continuous ammonia synthesis from water and nitrogen via contact electrification

We synthesized ammonia (NH3) by bubbling nitrogen (N2) gas into bulk liquid water (200 mL) containing 50 mg polytetrafluoroethylene (PTFE) particles (~5 µm in diameter) suspended with the help of a surfactant (Tween 20, ~0.05 vol.%) at room temperature (25 °C). Electron spin resonance spectroscopy and density functional theory calculations reveal that water acts as the proton donor for the reduction of N2. Moreover, isotopic labeling of the N2 gas shows that it is the source of nitrogen in the ammonia. We propose a mechanism for ammonia generation based on the activation of N2 caused by electron transfer and reduction processes driven by contact electrification. We optimized the pH of the PTFE suspension at 6.5 to 7.0 and employed ultrasonic mixing. We found an ammonia production rate of ~420 μmol L-1 h-1 per gram of PTFE particles for the conditions described above. This rate did not change more than 10% over an 8-h period of sustained reaction.

Spectroscopic and computational investigations of Cobalt(II) binding to the innate immune protein human calprotectin.

Human calprotectin (CP) is an innate immune protein that participates in the metal-withholding response to infection by sequestering essential metal nutrients from invading microbial pathogens. CP is comprised of S100A8 (α subunit, 10.8kDa) and S100A9 (β subunit, 13.2kDa). Two transition-metal binding sites of CP form at the S100A8/S100A9 dimer interface. Site 1 is a His3Asp motif comprised of His83 and His87 from the S100A8 subunit and His20 and Asp30 from the S100A9 subunit. Site 2 is an unusual hexahistidine motif composed of S100A8 residues His17 and His27 and S100A9 residues His91, His95, His103, and His105. In the present study, the His3Asp and His6 sites of CP were further characterized by utilizing Co2+ as a spectroscopic probe. Magnetic circular dichroism spectroscopy was employed in conjunction with electron paramagnetic resonance spectroscopy and density functional theory computations to characterize the Co2+-bound S100A8(C42S)/S100A9(C3S) CP-Ser variant and six site variants that allowed the His3Asp and His6 sites to be further probed. Our results provide new insight into the metal-binding sites of CP-Ser and the effect of amino acid substitutions on the structure of site 2.

Size and temperature dependent shapes of copper nanocrystals using parallel tempering molecular dynamics.

We performed parallel-tempering molecular dynamics simulations to predict the temperature- and size-dependent equilibrium shapes of a series of Cu nanocrystals in the 100- to 200-atom size range. Our study indicates that temperature-dependent, solid-solid shape transitions occur frequently for Cu nanocrystals in this size range. Complementary calculations with electronic density functional theory indicate that vibrational entropy favors nanocrystals with a shape intermediate between a decahedron and an icosahedron. Overall, we find that entropy plays a significant role in determining the shapes Cu nanocrystals, so studies aimed at determining minimum-energy shapes may fail to correctly predict shapes observed at experimental temperatures. We also observe significant shape changes with nanocrystal size - sometimes with changes in a single atom. The information from this study could be useful in efforts to devise processing routes to achieve selective nanocrystal shapes.

Spin states of metal centers in electrocatalysis.

Understanding the electronic structure of active sites is crucial in efficient catalyst design. The spin state, spin configurations of d-electrons, has been frequently discussed recently. However, its systematic depiction in electrocatalysis is lacking. In this tutorial review, a comprehensive interpretation of the spin state of metal centers in electrocatalysts and its role in electrocatalysis is provided. This review starts with the basics of spin states, including molecular field theory, crystal field theory, and ligand field theory. It further introduces the differences in low spin, intermediate spin, and high spin, and intrinsic factors affecting the spin state. Popular characterization techniques and modeling approaches that can reveal the spin state, such as X-ray absorption microscopy, electron spin resonance spectroscopy, Mössbauer spectroscopy, and density functional theory (DFT) calculations, are introduced as well with examples from the literature. The examples include the most recent progress in tuning the spin state of metal centers for various reactions, e.g., the oxygen evolution reaction, oxygen reduction reaction, hydrogen evolution reaction, carbon dioxide reduction reaction, nitrogen reduction reaction, nitrate reduction reaction, and urea oxidation reaction. Challenges and potential implications for future research related to the spin state are discussed at the end.

Orbital electron delocalization of axial-coordinated modified FeN4 and structurally ordered PtFe intermetallic synergistically for efficient oxygen reduction reaction catalysis.

Regulating the chemical environment of materials to optimize their electronic structure, leading to the optimal adsorption energies of intermediates, is of paramount importance to improving the performance of electrocatalysts, yet remains an immense challenge. Herein, we design a harmonious axial-coordination Pt x Fe/FeN4CCl catalyst that integrates a structurally ordered PtFe intermetallic with an orbital electron-delocalization FeN4CCl support for synergistically efficient oxygen reduction catalysis. The obtained Pt2Fe/FeN4CCl with a favorable atomic arrangement and surface composition exhibits enhanced oxygen reduction reaction (ORR) intrinsic activity and durability, achieving a mass activity (MA) and specific activity (SA) of 1.637 A mgPt -1 and 2.270 mA cm-2, respectively. Detailed X-ray absorption fine spectroscopy (XAFS) further confirms the axial-coupling effect of the FeN4CCl substrate by configuring the Fe-N bond to ∼1.92 Å and the Fe-Cl bond to ∼2.06 Å. Additionally, Fourier transforms of the extended X-ray absorption fine structure (FT-EXAFS) demonstrate relatively prominent peaks at ∼1.5 Å, ascribed to the contribution of the Fe-N/Fe-Cl, further indicating the construction of the FeN4CCl moiety structure. More importantly, the electron localization function (ELF) and density functional theory (DFT) further determine an orbital electron delocalization effect due to the strong axial traction between the Cl atoms and FeN4, resulting in electron redistribution and modification of the coordination surroundings, thus optimizing the adsorption free energy of OHabs intermediates and effectively accelerating the ORR catalytic kinetic process.

THERMODYNAMICS OF COMPLEXATION WITH 18-CROWN-6 WITH L-CARNOSINE AND ITS CONFORMATIONAL ANALYSIS

The study of model biological "guest-host" processes between dipeptides and crown ethers allows us to understand the nature of molecular recognition in living systems. The reaction of the formation of the molecular complex of the crown ether 18-crown-6 (18K6) with L-carnosine (L-Carn) in water at T = 298.15 K was studied using isothermal titration calorimetry. The mixing heats of aqueous solutions of L-carnosine with aqueous solutions of 18-Carn-6 were determined on a TAM III titration calorimeter system (TA Instruments) equipped with a 20 mL titration cell. In the main experiment, the heat of mixing of aqueous 18K6 solution with aqueous L-Carn solution was determined. In separate experiments the heat effect of dilution of L-Carn solution in water was measured, which was taken into account in determining the thermal effect of complex formation. From the thermochemical data, the thermodynamic parameters of the molecular complex formation reaction were calculated. It was found that the presence of a histidine fragment in the L-Carn structure leads to an increase in the stability of the molecular complex L-Carn with 18K6 and to a decrease in the exothermicity of complex formation as compared to the complex formation of alanine with 18K6. Thermodynamic parameters of [L-Carn 18K6] formation reaction are similar to thermodynamics of 18K6 complex formation with histidine. Quantum-chemical calculations were performed within the framework of the electron density functional theory using the GAUSSSIAN 03 software package. The results were visualized using the ChemCraft software. We performed the conformational analysis of L-Carn within the framework of the density functional theory (version B3LYP/6-31G**) to analyze the reactivity of L-Carn structural fragments. For three stable conformers in the M062X/ aug-cc-pVTZ approximation the geometric parameters were obtained and the electron density distribution was analyzed within the framework of the QTAIM method as well as the NBO scheme. For citation: Usacheva T.R., Belova N.V., Saturina E.V., Krutova O.N., Lukanov M.M., Pavlova E.A. Thermodynamics of complexation with 18-crown-6 with l-carnosine and its conformational analysis. ChemChemTech [Izv. Vyssh. Uchebn. Zaved. Khim. Khim. Tekhnol.]. 2023. V. 66. N 5. P. 21-31. DOI: 10.6060/ivkkt.20236605.6781.

In Ukrainian

The results of theoretical modeling of the spatial structure and electronic absorption spectra of isolated small particles of polymorphic modifications of lead monoxide - litharge and massicot, which are characterized by the presence of hypercoordinated atoms, have been analyzed. The optical properties of many oligomers of lead oxide are systematically considered in the work. The results of quantum chemical calculations by the methods of electron density functional theory (DFT and TDDFT) on the spatial structure and energy characteristics, as well as the electronic absorption spectra of the considered models, are discussed. Theoretical results are compared with available experimental data. The calculation method used and the proposed cluster models with or without attached (chemisorbed) water molecules allow us to reproduce the spatial structure and energy characteristics of polymorphic modifications of lead oxide at a semi-quantitative level. It has been found that the addition of water molecules to models for α-lead oxide nanoparticles contributes to their stabilization. The calculated values of the cohesion energy for litharge nanoparticles are greater than those for massicot ones, which corresponds to the experimental data for lead monoxide crystals. The calculated electronic spectra of the litharge nanoparticle models are characterized by absorption bands which lie in the region that corresponds to the red color, and those of massicot relate to yellow; the corresponding values of the width of the energy gap for the litharge models are smaller than for the massicot ones, which is consistent with the experimental data for crystals. The results of calculations show that the presence of hydroxyl groups in nanoparticle models leads to a hypsochromic shift of absorption maxima, so that with a sufficient number of such groups, these species may lose their color.

Polar Polymorphism: A New Intermediate Structure toward the Thin-Film Phase in Asymmetric Benzothieno[3,2-b][1]-benzothiophene Derivatives.

Understanding structure and polymorphism is relevant for any organic device optimization, and it is of particular relevance in 7-decyl-2-phenyl[1]benzothieno[3,2-b][1]benzothiophene (Ph-BTBT-10) since high carrier mobility in Ph-BTBT-10 thin films has been linked to the structural transformation from the metastable thin-film phase to the thermodynamically stable bilayer structure via thermal annealing. We combine here a systematic nanoscale morphological analysis with local Kelvin probe force microcopy (KPFM) that demonstrates the formation of a polar polymorph in thin films as an intermediate structure for thicknesses lower than 20 nm. The polar structure develops with thickness a variable amount of structural defects in the form of individual flipped molecules (point defects) or sizable polar domains, and evolves toward the reported nonpolar thin-film phase. The direct experimental evidence is supported by electronic structure density functional theory calculations. The structure of the film has dramatic effects on the electronic properties, leading to a decrease in the film work function (by up to 1 eV) and a considerable broadening of the occupied molecular orbitals, attributed to electrostatic disorder. From an advanced characterization point of view, KPFM stands out as a valuable tool for evaluating electrostatic disorder and the conceivable emergence of polar polymorphs in organic thin films. The emergence of polar assemblies introduces a critical consideration for other asymmetric BTBT derivatives, which may be pivotal to understanding the structure-property relationships in organic field-effect transistors (OFETs). A precise determination of any polar assemblies close to the dielectric interface is critical for the judicious design and upgrading of high-performance OFETs.

A Density Potential Functional Approach for Mesoscopic Metal-Solution Interfaces

Mesoscopic metal-solution interfaces (MSI) are found at nanostructured electrodes in myriad electrochemical devices. Resolving the local reaction environment for electrochemical reactions occurring in mesoscopic MSIs requires a constant-potential, computationally economic method that properly treat both quantum mechanical metal electrons and almost classical electrolyte components. We have been developing a density potential functional theory (DPFT) to meet these harsh requirements [1-4]. Our DPFT method combines an orbital-free density functional theory of quantum mechanical metal electrons, a statistical field theoretic description of classical electrolyte ions and solvent, and Morse potentials describing short-range interactions between the two parts that are parameterized using first-principles based computations. In this talk, I will introduce the theoretical framework, its relationship to other similar methods including the jellium model and implicit solvation model, its parameterization and benchmark with experimental and ab initio simulation data, and recent applications to understand solvent effects on potential of zero charge, mechanoelectrochemistry of MSIs, and supported nanoparticle catalysts. The emerging opportunities, future development needs, and important caveats of the DPFT method will be discussed. Reference Huang J, Zhang Y, Li M, Groß A, Sakong S. Comparing Ab Initio Molecular Dynamics and a Semiclassical Grand Canonical Scheme for the Electric Double Layer of the Pt (111)/Water Interface. The Journal of Physical Chemistry Letters. 2023 Feb 27;14:2354-63.Huang J. Density-Potential Functional Theory of Electrochemical Double Layers: Calibration on the Ag (111)-KPF6 System and Parametric Analysis. Journal of Chemical Theory and Computation. 2023 Jan 18.Huang J, Chen S, Eikerling M. Grand-canonical model of electrochemical double layers from a hybrid density–potential functional. Journal of chemical theory and computation. 2021 Mar 31;17(4):2417-30.Huang J. Hybrid density-potential functional theory of electric double layers. Electrochimica Acta. 2021 Sep 1;389:138720. J.H. is supported by the Helmholtz Young Investigators Group Leader Program.

Atomically Unraveling Highly Crystalline Self-Intercalated Tantalum Sulfide with Correlated Stacking Registry-Dependent Magnetism.

In self-intercalated two-dimensional (ic-2D) materials, understanding the local chemical environment and the topology of the filling site remains elusive, and the subsequent correlation with the macroscopically manifested physical properties has rarely been investigated. Herein, highly crystalline gram-scale ic-2D Ta1.33S2 crystals were successfully grown by the high-pressure high-temperature method. Employing combined atomic-resolution scanning transmission electron microscopy annular dark field imaging and density functional theory calculations, we systematically unveiled the atomic structures of an atlas of stacking registries in a well-defined √3(a) × √3(a) Ta1.33S2 superlattice. Ferromagnetic order was observed in the AC' stacking registry, and it evolves into an antiferromagnetic state in AA/AB/AB' stacking registries; the AA' stacking registry shows ferrimagnetic ordering. Therefore, we present a novel approach for fabricating large-scale highly crystalline ic-2D crystals and shed light on a powerful means of modulating the magnetic order of ic-2D systems via stacking engineering, i.e., stackingtronics.