2026-09-29 | | Total: 23
The objective of this study is to simulate the propagation of the shock wave in water due to an explosion. The study is part of a global research program on the development of an alternative stimulation technique to conventional hydraulic fracturing in tight gas reservoirs aimed at inducing a distributed state of microcracking of rocks instead of localized fracture. We consider the possibility of increasing the permeability of rocks with dynamic blasts. The blast is a shock wave generated in water by pulsed arc electrohydraulic discharges. The amplitude of these shock waves is prescribed by the electrohydraulic discharges which generate high pressures of several kilobars within microseconds. A simplified method has been used to simulate the injected electrical energy as augmentation of enthalpy in water locally. The finite element code EUROPLEXUS is used to perform fluid fast dynamic computation. The predicted pressure is consistent with the experimental results. In addition, shock wave propagation characteristics predicted with simulation can be valuable reference for design of underwater structural elements and engineering of underwater explosion.
Strong light-matter coupling in optical cavities offers a non-invasive route to modify molecular properties, and it is usually reached by collectively coupling many molecules to the same mode. Here we consider the photonic counterpart of the collective regime, in which a single molecule interacts with a cavity mode holding many photons, a steady state of the pumped cavity. Since the eigenstates of the Pauli-Fierz Hamiltonian carry no transverse electric field, such a state cannot be described as an excited state of the light-matter system. We therefore constrain the photon number of the strong coupling QED Hartree-Fock wave function with a Lagrange multiplier and self consistently minimize the energy with respect to all parameters. The constrained reference carries the field of a coherent state with an average number of photons, generated by a cavity mode whose effective free-field frequency is lowered by the multiplier. We apply this framework to benzene and a benzene-water complex, where the pumped cavity polarizes the molecules and modifies their interaction as with a static classical field. For hydrogen peroxide, the field can reshape the torsional potential and the far-infrared torsional spectrum. The classical driven field limit thus emerges from a quantized mean field description, and the photon density becomes a control parameter for cavity-modified chemistry.
Molecular crystal packing shapes properties from drug bioavailability to charge transport, yet generating realistic structures requires coordinating molecular flexibility, intermolecular interactions and symmetry. Here we introduce SALA, a flow-matching model built on state--context separation: it evolves only asymmetric-unit (ASU) atomic coordinates and symmetry-compatible lattice parameters while predicting their updates from the full periodic environment. Conditioned on molecular graphs and a space group, SALA generates molecular conformation, lattice geometry and packing end to end while preserving the specified symmetry. Across 799 held-out crystals spanning four chemical classes and all seven crystal systems, SALA achieves a 72.0\% hit rate from 50 candidates per target, compared with 8.1\% for an atomistic full-cell generative baseline. Mean best-candidate packing root-mean-square deviation decreases from 4.48 to 1.70~Å, with the advantage widening as system size and symmetry multiplicity increase. SALA thus enables end-to-end generation of complex, symmetry-constrained molecular crystals, with potential applications in CSP, structure-based pharmaceutical polymorph design and organic molecular materials discovery.
Generative models produce three-dimensional molecular structures with high geometric accuracy, yet they are trained on distributions of atomic coordinates without explicit access to the underlying physics. The potential energy surface (PES) describes how the energy changes around a molecular geometry, and the RMSD between a generated structure and its reference captures this local landscape only indirectly. Supervising the PES during training could close this gap, but conformers and transition states (TSs) are stationary points at which the forces vanish, so first-order quantities carry little information about them. What remains is the local curvature, encoded in the Hessian of the PES. Quantum-chemical Hessians are too expensive to be available at the scale required for training, but recently proposed machine-learned predictors make curvature inexpensive enough to supervise directly. We introduce CurvFM, a flow matching (FM) formulation whose training objective constrains the curvature at the generated structure. A frozen Hessian prediction model is evaluated at the generated endpoint and at the reference structure, which removes the need for quantum-chemical Hessians and keeps the residual cheap enough to impose at every training step. Across TS and conformer generation, CurvFM substantially improves physical fidelity at unchanged global structural accuracy, reducing the median deviation in maximum force by 63% on GEOM-QM9 and raising the fraction of valid TSs by up to 10% across three reaction datasets.
We present NNLap, a machine-learned Laplacian-level exchange-correlation (XC) functional that augments PBE with a neural-network correction depending on the electron density, its gradient, and its Laplacian. The model is trained on exact XC potentials and energies, obtained through inverse density-functional theory (DFT) calculations on configuration-interaction densities. Despite training on only a few systems -- five atoms and three molecules -- the model achieves remarkable accuracy on thermochemistry benchmarks, competing with the meta-GGA functionals SCAN and r2SCAN. It also attains accurate total energies, comparable to SCAN and better than r2SCAN and B3LYP. This shows that a Laplacian-level model, trained on exact XC potentials and energies, can reach the accuracy of meta-GGAs without their orbital dependence.
Floquet theory provides insight into the inversion of time-dependent Kohn-Sham density functional theory. Specifically, mathematical derivations show that the fundamental frequencies of a time-dependent wavefunction solution are the leading-order harmonics for the TD-KS state. Numerical tests of the resulting ansatz in 1D and 3D for atomic and molecular cases demonstrate its utility. In particular, low $L_{2}$ errors in the time-dependent density and longitudinal current were found, even though currents were not an explicit optimization objective. In all, the proposed inversion ansatz provides exchange-correlation potentials from time-dependent wavefunctions, is highly interpretable, and may significantly help in the development of nonadiabatic density functionals.
Quantum computing offers a promising route for electronic-structure calculations, but practical condensed-matter applications often require solving large families of related Hamiltonians arising from $\mathbf{k}$-point sampling, compositional variation, defects, and changes in system size. Here, we develop strategies for transferring variational information between related Hamiltonians to accelerate quantum simulations. For Hamiltonians of the same embedded size, optimized variational parameters are transferred between related calculations. For different sizes, we introduce Augmented Ansatz Reuse for Target Initialization (AARTI), which transfers compatible Pauli-generator structure and optimized source parameters, then augments the ansatz with target-specific generators. Our workflow combines density functional theory, Wannier downfolding, embedded tight-binding Hamiltonians, variational quantum algorithms, and neural-network quantum states, with quantum-circuit simulations implemented using NVIDIA's CUDA-Q platform. Using Li$_x$CoO$_2$ as a representative battery cathode material, we study Hamiltonian families spanning multiple delithiation levels, dense $\mathbf{k}$-point meshes, and increasing sizes. Parameter and ansatz transfer substantially reduce the optimization steps required to solve these related Hamiltonians while maintaining the common target accuracy, with lower terminal errors in several cases. The calculated electronic-structure evolution with delithiation exhibits pronounced non-rigid-band behavior consistent with experiment. These results establish a framework for exploiting similarity among related Hamiltonians to improve the efficiency and scalability of quantum workflows for realistic condensed-matter systems.
An improved algorithm to implement our previously proposed static quantum embedding method MPCC (J. Chem. Phys. 161, 164107 (2024)) is presented. It uses the density-fitting (DF) approximation for the electron-repulsion integrals (ERIs) across all layers of the embedding method: the low-level method, the screened interaction, and the high-level method. Additionally, an approximate factorized solution of the Sylvester equation is employed based on the Laplace transformation for the low-level method. To build the screened interaction, we have explored an approximate construction of it based on the distinguishable cluster approximation (DCA) by Kats and Manby (J. Chem. Phys. 141, 061101 (2014)), which amounts to neglecting certain higher-order exchange and ladder diagrams. Combining all these techniques, we have obtained an $\mathcal{O}(N^4)$ algorithm for the low-level method; a non-iterative $N_{\rm frag} \mathcal{O}(N^4)$ algorithm for the screened interaction; and an iterative $\mathcal{O}(N_{\rm frag}^6)$ algorithm for the fragment solver, where $N_{\rm frag}$ is the number of fragment orbitals and $N$ is the total number of orbitals. The resulting low-scaling implementation of the MPCC method is applied to trans polyacetylenes, that is, C$_{2n}$H$_{2n+2}$ molecules for chain lengths $n=1-8$, and the timing and accuracy of the method are compared to the original unfactorized implementation. Moreover, the approximations are further validated by applying the method to the potential energy surface of the \ce{O3} molecule. Finally, we apply the method to the selected molecules in the S22 dataset for which the performance of the MP2 method is known to be poor. The results demonstrate that the DF-based implementation of the MPCC method provides a significant reduction in computational cost while maintaining accuracy comparable to that of the original implementation.
Chemistry's phenomenological theories multiply with its observables: the same physical reality is repeatedly named in different measurement contexts, and the relations among theories have never had an auditable carrier. Here we present PhotoLean, a machine-checked genealogy of seventeen photochemical and photophysical theories formalized in Lean 4 over a shared kernel. Every pair of theories carries either a verified relation (equivalence, implication, composition with registered premises, definitional certificate) or a registered reason for its absence, covering all 136 pairs. Three adjudications resolve inherited conflations with if-and-only-if boundaries and witnesses: the symmetry factor $β$ equals one-half exactly when reactant and product force constants coincide; Kasha's rule and the Kasha-Vavilov rule agree only under closed quantification with a loss channel; and intensity-only Stern-Volmer data cannot distinguish static from dynamic quenching, whereas the lifetime channel discriminates exactly. As proposition generation shifts to machines, such machine-checked relational ledgers offer a replicable audit layer for theory itself.
Transition-state searches remain a major bottleneck in reaction discovery, as identifying valid saddle-point structures requires numerous expensive quantum-chemical calculations. Generative models can reduce this burden by proposing candidates from reactant and product geometries, but supervised training on geometries alone does not enforce the physical conditions required of a transition state, while enforcing them during sampling is costly. We introduce Trèfle, a one-step flow model that amortizes physics-guided generation into training: a mean-flow objective replaces multi-step sampling, a low-cost surrogate replaces density functional theory (DFT) reward evaluations, and reinforcement-learning post-training with rewards that enforce those conditions replaces inference-time guidance. Trèfle generates accurate transition states at $30\times$ the throughput of multi-step baselines. It transfers few-shot to ten unseen transition metals and generalizes zero-shot to molecules far larger than any in training. With DFT refinement and reaction-path validation, Trèfle-seeded searches recover all 13 unseen $γ$-ketohydroperoxide channels, while post-training raises per-attempt recovery from 28% to 50% across 85 unseen bimolecular channels. Trèfle thus provides fast, reliable seeds for automated reaction discovery.
Molecular dynamics simulations provide a computational microscope for atomic scale processes but remain restricted to relatively small spatial and temporal scales. Coarse-grained models extend their reach by representing groups of atoms as effective interaction sites. The number of particles is systematically reduced by embedding the collective behavior of groups of atoms into coarse-grained sites governed by effective potentials. However, constructing accurate, highly coarsened potentials, even with the assistance of modern machine-learning methods, is a grand challenge because coarse grained forces depend on a statistical distribution of composite atomic configurations making direct calculation computationally prohibitive. Here, we introduce an alternative framework based on implicit score matching that directly avoids the need for explicit coarse-grained forces. This approach both substantially reduces the training data required relative to conventional force-matching methods and enables high-fidelity models in which each coarse-grained site represents hundreds to thousands of atoms. We demonstrate the method's versatility and effectiveness across a diverse range of chemical systems and phenomena, developing bottom-up, machine-learned models. These enable efficient modeling at micrometer and millimeter length scales, while retaining significant amounts of the atomic scale fidelity. These results establish implicit score matching as a practical route towards chemically accurate simulations that bridge molecular and mesoscopic scales.
The ultrafast photochemistry of heterocyclic compounds is central to photobiology and materials chemistry, yet direct experimental observation of their structural dynamics remains rare. Here we present the first real-time structural characterization of photoinduced ring opening and subsequent fragmentation in the isomeric pair oxazole and isoxazole using MeV ultrafast electron diffraction (UED), complemented by non-adiabatic molecular dynamics simulations. Upon photoexcitation at 200 nm, both isomers undergo ring opening followed by fragmentation into various products, but with strikingly different dynamics governed by heteroatom positioning. Trajectory Surface Hopping simulations reproduce the experimental diffraction signatures, which are surprisingly similar for both isomers, and reveal distinct mechanistic pathways for the two isomers. For isoxazole, all trajectories exclusively undergo N-O bond cleavage within 40 fs, followed by sequential fragmentation into HCN + ketene and HCO + vinyl nitrene channels on the hundreds-of-femtoseconds timescale. Oxazole, by contrast, shows significantly slower ring opening (290 fs) with only 85% efficiency, proceeding primarily through O-C cleavage and accessing a richer landscape of intermediates including nitrile ylide and O-pyramidalized structures. Simulated diffraction patterns derived from trajectory ensembles, convolved with the experimental instrument response function, are in agreement with the branching ratios observed by UED. This synergy between UED and trajectory surface hopping provides an atomistic picture of how the simple interchange of heteroatom connectivity in structural isomers fundamentally reshapes excited-state potential surfaces, conical intersection accessibility, and photochemical outcome.
AI agents are advancing the automation of molecular dynamics (MD) research, with the potential to accelerate discovery and design of lithium-battery electrolytes. However, whether they can autonomously complete scientifically reliable electrolyte MD studies remains insufficiently evaluated. We introduce ElectrolyteMD-Bench, which unifies dilute, high-concentration, and localized high-concentration electrolytes in a comparative study across distinct solvation regimes. The benchmark connects local coordination, ionic association, single-particle dynamics, and collective transport within one continuous research process. Agents autonomously perform system construction, equilibration, production sampling, property analysis, and interpretation under a frozen scientific model, while deciding whether to correct errors, continue sampling, or stop. We audit each study along four dimensions: simulation validity, analysis validity, evidence integrity, and outcome calibration. Across eight model--harness configurations, all 24 system runs produced valid production trajectories, yet none of the eight studies satisfied all scientific requirements. Failures included missing analyses, incorrect physical definitions or numerical implementations, insufficient statistical support, and completion claims inconsistent with retained evidence. These results expose a gap between successful MD workflow execution and reliable scientific study completion. ElectrolyteMD-Bench therefore provides a foundational test and defines a key capability threshold that autonomous electrolyte simulation must cross to advance toward computation-first materials exploration.
Methane pyrolysis experiments were performed in a micro flow reactor to investigate the chemical mechanism sensitivities at elevated pressures and temperatures. The micro flow reactor used in this study provides independent control over reaction parameters. Experiments were performed at temperatures from 1400 to 1800 K and pressures from 1 to 10 bar at a fixed residence time of 60 ms. The yields of hydrogen, unreacted methane, and solid carbon were quantified for comparison to predictive models. The mass of heterogeneous depositions was also quantified. The presence of C2 and PAH species were measured by GC-MS. The results were compared to 1D simulations using a discrete sectional model for capturing the kinetics of particulate formation and surface deposition. The simulations replicated the pressure trend well at high temperatures (1800 K), but the sensitivity to temperature was in disagreement. Simulations generally found lower carbon and hydrogen yields and a stronger sensitivity to temperature than the experiment. An abundance of PAH intermediates was measured at lower temperatures and higher pressures. The simulated surface deposition rates were lower than the measured values but showed similar profiles. The results demonstrate the utility of micro flow reactors in recreating pyrolysis environments for hydrogen production from methane. The predictive capabilities of simulations could be improved by advancing the fuel decomposition and PAH formation and consumption kinetics at low temperatures and by detailed treatment of surface site density.
Several researchers have reported an atypical inversion of the lowest singlet ($S_{1}$) and triplet ($T_{1}$) excited states in nitrogen-substituted planar hydrocarbons. Contrary to Hund's rule of maximum multiplicity, poly-aza derivatives of phenalenyl, such as cyclazine and heptazine, have long been known to possess a negative singlet-triplet gap, $ΔE_{ST}=E(S_{1})-E(T_{1})$. Using the well-known Pariser-Parr-Pople (PPP) model, a semi-empirical, wave-function-based Hamiltonian for correlated $π$-electrons, we investigate the singlet-triplet gap in ten N-substituted phenalenyl motifs. Two schemes of long-range Coulomb interactions are considered, and key parameters, including the on-site orbital energy, Hubbard U, and hopping integrals, are systematically tuned to benchmark the model for reliably estimating negative singlet-triplet gaps. Full configuration-interaction (FCI) method is employed to incorporate electron-correlation effects. Ground-state and $T_{1}$ absorption spectra are subsequently calculated using the optimized PPP parameters that best reproduce the available experimental data. For validation, first-principles TDDFT calculations using various exchange--correlation functionals are also performed. Conventional adiabatic TDDFT fails to reproduce the inverted ($S_{1}-T_{1}$) ordering, indicating that the inclusion of the correlation driven multiconfigurational effects are essential for describing it. We find that the singlet TDDFT ground-state absorption spectra reproduce several qualitative features of the PPP-FCI spectra, although the TDDFT excitation energies are systematically higher. Detailed analysis of the ground-state and $T_{1}$ absorption spectra reveals significant configuration mixing, providing a comprehensive characterization of the optical response of the N-substituted phenalenyl derivatives studied.
Weak, non-covalent interactions such as electrostatics, induction, and dispersion are commonly treated as distinct physical mechanisms, often using separate models and approximations. Here, we present a unified quantum-electrodynamical framework in which all weak interactions between neutral or charged subsystems emerge from a single interaction Hamiltonian and are expressed entirely in terms of electromagnetic Green functions. Starting from macroscopic QED in arbitrary linear environments, we derive compact master formulas that generate permanent electrostatic, mixed induction, and fully induced (dispersion) interactions on equal footing. Monopoles, dipoles, and quadrupoles are treated consistently within the same multipole expansion, leading to a common structural form in which all interaction channels are obtained from contractions of multipole moments and response functions with Green tensors. Environmental effects enter universally through dressed Green functions. The standard non-retarded free-space limits, including Keesom, Debye, and London/van der Waals interactions, are recovered as a consistency check of the formalism. Beyond this validation, the formulation shows that once the general expressions are established, different physical scenarios, such as interfaces, confinement, or polarisable continua, can be treated by direct substitution of the corresponding Green function, without reformulating the interaction for each case. This elevates Green functions to a universal organising principle for weak interactions in complex environments.
The hierarchical equations of motion (HEOM) approach is one of the most powerful methods avail- able to simulate the dynamics of open quantum systems. Although the HEOM approach has been shown to be stable and efficient in many parameter regimes and physical scenarios, it is also well known that the unavoidable finite truncation of the hierarchy can lead to fundamental numeri- cal instabilities for bosonic environments. In this work, we apply recently developed techniques to analyze these issues. We use an auxiliary Fock-space picture of the hierarchy to first identify the unbalanced hierarchy-raising terms responsible for the strong transient amplification associated with the non-normal algebraic structure of the HEOM generator. We then apply a non-unitary similarity transformation that converts these terms into balanced combinations of hierarchy-raising and hierarchy-lowering terms, thereby suppressing their numerical impact. We first demonstrate the efficacy of this approach via a spin-boson model with a simple Brownian oscillator spectral den- sity and subsequently extend the analysis to a structured spectral density, showing for both cases that the transformed HEOM approach allows simulations of much larger system-bath coupling than would otherwise be possible.
We develop a rational-function approximation for the structural properties of multicomponent fluids with narrow square-well and square-shoulder interactions, together with weakly nonadditive hard-sphere mixtures. The theory is obtained by reformulating the analytical Percus--Yevick solution for additive sticky hard-sphere mixtures so as to account for finite interaction widths, leading to a unified semianalytical framework for these classes of fluids. The approximation satisfies the exact low-density limit and preserves the continuity of the cavity functions; a simple local correction restores the continuity of their first derivatives. Comparison with Monte Carlo simulations for binary mixtures shows that the theory accurately describes both the radial distribution functions and their Fourier transforms for narrow wells and shoulders and for weak nonadditivity. The proposed formulation extends previous sticky-hard-sphere-based approaches while retaining their analytical simplicity.
Helium dimer in its metastable triplet electronic state, He$_2(a~^3Σ_u^+)$, is a promising candidate for direct laser cooling and precision measurements. Interactions and collisions involving He$_2(a~^3Σ_u^+)$ remain unexplored in ultracold and cold energy domains. Here, we present the first theoretical investigation of intermolecular interactions and collision dynamics of He$_2(a~^3Σ_u^+)$+Ne van der Waals complex. For the ground electronic state, $X~^3\mathrm{A}'$, we employ a spin-restricted coupled cluster method with single, double, and noniterative triple excitations. The resulting interaction potential exhibits an extremely shallow well with a depth of 4.612 cm$^{-1}$. We construct a three-dimensional vibrationally averaged potential energy surface. Using this surface, we perform full close-coupling calculations including fine-structure levels of He$_2(a~^3Σ_u^+)$ to investigate collisional excitation and quenching induced by Ne. We analyze the low-energy scattering dynamics and model resonant features associated with quasi-bound states supported by interaction potential. Propensity rule for fine-structure conserving transitions is established. We additionally consider He$_2$+Ar collisions and compare their energy dependence and resonance structure with those of He$_2$+Ne. The He$_2$+Ar system exhibits a richer resonance structure, highlighting the sensitivity of low-energy scattering to interaction potential and reduced mass. Calculated PES and scattering properties provide a reliable benchmark for description of He$_2$+Ne van der Waals complex and may be useful for interpretation of future experimental spectra and collision measurements.
We implement the machine-learned Skala exchange-correlation (XC) functional natively in CP2K using its Gaussian and plane-wave (GPW) and Gaussian and augmented-plane-wave (GAPW) methods. A joint one-centre reconstruction of density, density gradients, and kinetic-energy density before functional evaluation preserves mixed gradient terms and nonlocal couplings between smooth and atom-local contributions in all-electron and pseudopotential calculations. An XC-specific GAPW representation resolves rapidly varying local contributions on atom-centred grids, reducing plane-wave requirements for pseudopotential calculations. The implementation includes Brillouin-zone sampling and point-group symmetry reduction. With sufficiently flexible orbital bases, it retains the molecular benchmark accuracy of our earlier GauXC formulation. All-electron Skala-D3(BJ) calculations yield a mean absolute error (MAE) of 1.54 kJ mol$^{-1}$ against diffusion Monte Carlo for crystalline CO$_2$, NH$_3$, and urea, spanning dispersion, quadrupolar electrostatics, and hydrogen bonding. For all thirteen DMC-ICE13 phases, the MAEs are 1.16 kJ mol$^{-1}$ for absolute lattice energies and 0.82 kJ mol$^{-1}$ for the twelve relative energies to ice Ih. At fixed experimental geometries, all-electron/mixed-core calculations give a band-gap MAE of 0.43 eV on a 15-material set, substantially below commonly used semilocal functionals and comparable to widely used hybrid functionals. The LC10 benchmark nevertheless reveals systematically underestimated equilibrium lattice constants, indicating structural overbinding. This framework connects molecular Skala to condensed-phase electronic structure and enables future development using periodic many-body reference data.
We present a first-quantization, all-electron, full-potential Kohn-Sham density functional theory method for quantum computers, based on a classically optimized quantum-circuit eigensolver for the self-consistent field (SCF) eigenproblems in any orthonormal basis. The plane-wave-orthonormalized-orbital (PWOO) basis combines plane waves with atomic orbitals (AOs) orthonormalized against them and one another. The eigensolver has a dual-basis form: after full-basis bootstrap steps, it solves each SCF step in a contracted basis of $K$ functions ($\lceil\log_2 K\rceil$ qubits), updated in the full basis when the residual outside it exceeds a threshold. For an H$_8$ chain (SVWN, $K=8$), the state-vector-emulated SCF agrees with classical calculations within 0.02 mHa for 32768 plane waves and PWOO bases up to cc-pVTZ; contraction changes these PWOO energies by at most 5.7 $μ$Ha. With emulated shot noise, the STO-3G PWOO SCF converges near the noise-free fixed point. On IonQ Forte-1, three-qubit-register readout reproduces the noise-free PWOO-AO energy difference within one standard deviation.
Structural observations in polymer chemistry are often rationalized through underlying mechanistic insights. The study of families of mechanisms - and the phenomenology and lawlike behavior that can accompany them - is less common. Here we construct a theory that links a mechanistic feature - the number of nearest neighbors $k$ that influence monomer incorporation - with the accessible copolymer sequence space. We find that the reaction mechanisms can be mapped 1-to-1 to a regular de Bruijn graph of degree $k$. We calculate the periodic copolymer sequences as function of neighbor index $k$, and number of monomers $s$. We show that the number of periodic sequences grows explosively with $s,k$, bounded from below by $s!^{s^{k-1}}/s^k$. We find that periodic k-polymerization is governed by several simple laws describing allowed sequences, number of specific reactions, repeating unit length and the allowed and disallowed sequence coexistence in block copolymers. We formulate Simple Assembly Rules (SAR) that capture commonly observed behavior and design of literature systems, and characterize the possible k-polymerizations they are allowed to access ($k\leq2$) as function of symmetries and illustrate our findings with reported examples from the coordination chemistry and supramolecular chemistry literature. We discuss chemical-structural contributions by which nontrivial ($k>1$) sequences are formed, highlighting the importance steric hindrance, long-range interactions and degrees of freedom.
In this work, we introduce a precursor-level strategy by substitutionally doping the Ti3AlC2 MAX precursor with a small amount of Yttrium during top-down MAX preparation. This allows permanent structural changes that are retained in the resulting Ti3C2Tx MXene. Density functional theory (DFT) calculations identify the outer Ti(4f) site, facing the Al layer, as the preferred incorporation site for Y under the Al-excess synthesis conditions used here, rather than the inner Ti(2a) or Al(2b) sites. This site-specific substitution is consistent with the experimentally observed lattice expansion and carbide-like Y bonding and provides a microscopic basis for the structural changes retained in the resulting Ti3C2Tx MXene. Y incorporation expands the MXene interlayer spacing due to an increase in the amount of hydroxyl surface terminations, creating more favorable pathways for lithium-ion transport. Comprehensive structural characterization (including XRD, XPS and TEM) confirms these modifications. As a result, the Y-doped Ti3C2Tx cathode shows diffusion-controlled charge storage with a Li diffusion coefficient of 10-9-10-11 cm2/s. It delivers a reversible capacity of about 130 mAh/g at 0.2 C and retains 74 mAh/g after 1000 cycles at 2 C with nearly 100% Coulombic efficiency. The electrode also achieves a high power density of up to 3970 W/kg, outperforming the undoped MXene and previously reported multilayered MXene cathodes. These results show that substitutional doping of the MAX precursor is a promising strategy for controlling the interlayer structure of MXenes and improving ion transport, providing a new route for the development of high-power energy-storage electrodes.