Search NASA⌕ Search

SEARCH · Search NASA

Results for “zero-point vibrational energy”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

17 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

The Remarkable I 2 O 3 Molecule: A New View from Theory

Atmospheric iodine chemistry has garnered increasing attention as a result of increased iodine emissions. A key subset of this chemistry involves iodine oxides (I 2 O 2–5 ), which serve as precursors to particle formation. Among these, I 2 O 3 is the simplest iodine oxide involved in particle formation, but it has remained undetected in the atmosphere. Previous theoretical studies have characterized this peculiar molecule, primarily using energies to refine geometries obtained at low levels of theory. Due to the reemerging interest in I 2 O 3 , this study presents geometries optimized at the CCSD(T)/aug-cc-pwCVTZ-PP level of theory─marking the first instance, to the best of our knowledge, where this system has been studied exclusively with CCSD(T). Harmonic vibrational frequencies were computed at the same level of theory. Final energetics were obtained using the very high level CCSDT(Q) method with basis sets up to quintuple-zeta cardinality (aug-cc-pwCV5Z-PP) and extrapolated to the CBS limit to yield CCSDT(Q)/CBS//CCSD(T)/aug-cc-pwCVTZ-PP energies. These energies include harmonic zero-point vibrational energy corrections and scalar relativistic energy corrections. Additionally, this study discovers new isomers along the I 2 O 3 potential energy surface, a novel contribution to the field. The performance of different computational methods and DFT functionals commonly used in atmospheric chemistry is also assessed relative to high-level theoretical methods.

basis sets↗
DOE OSTI2023

Analytical harmonic vibrational frequencies with VV10-containing density functionals: Theory, efficient implementation, and benchmark assessments

VV10 is a powerful nonlocal density functional for long-range correlation that is used to include dispersion effects in many modern density functionals, such as the meta-generalized gradient approximation (mGGA), B97M-V, the hybrid GGA, ωB97X-V, and the hybrid mGGA, ωB97M-V. While energies and analytical gradients for VV10 are already widely available, this study reports the first derivation and efficient implementation of the analytical second derivatives of the VV10 energy. The additional compute cost of the VV10 contributions to analytical frequencies is shown to be small in all but the smallest basis sets for recommended grid sizes. Here, this study also reports the assessment of VV10-containing functionals for predicting harmonic frequencies using the analytical second derivative code. The contribution of VV10 to simulating harmonic frequencies is shown to be small for small molecules but important for systems where weak interactions are important, such as water clusters. In the latter cases, B97M-V, ωB97M-V, and ωB97X-V perform very well. The convergence of frequencies with respect to the grid size and atomic orbital basis set size is studied, and recommendations are reported. Finally, scaling factors to allow comparison of scaled harmonic frequencies with experimental fundamental frequencies and to predict zero-point vibrational energy are presented for some recently developed functionals (including r2SCAN, B97M-V, ωB97X-V, M06-SX, and ωB97M-V).

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

A matrix completion algorithm for efficient calculation of quantum and variational effects in chemical reactions

This work examines the viability of matrix completion methods as cost-effective alternatives to full nuclear Hessians for calculating quantum and variational effects in chemical reactions. The harmonic variety-based matrix completion (HVMC) algorithm, developed in a previous study (https://doi.org/10.1063/5.0018326), exploits the low-rank character of the polynomial expansion of potential energy to recover, using a small sample, vibrational frequencies (square roots of nuclear Hessian eigenvalues) constituting the reaction path. Furthermore, these frequencies are essential for calculating rate coefficients using variational transition state theory with multidimensional tunneling (VTST-MT). HVMC performance is examined for four SN2 reactions and five hydrogen transfer reactions, with each H-transfer reaction consisting of at least one vibrational mode strongly coupled to the reaction coordinate. HVMC is robust and captures zero-point energies, vibrational free energies, zero-curvature tunneling, and adiabatic ground state and free energy barriers as well as their positions on the reaction coordinate. For medium to large reactions involving H-transfer, with the exception of the most complex Ir catalysis system, less than 35% of total eigenvalue information is necessary for accurate recovery of key VTST-MT observables.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

High-Accuracy Heats of Formation for Alkane Oxidation: From Small to Large via the Automated CBH-ANL Method

It is generally challenging to obtain high-accuracy predictions for the heat of formation for species with more than a handful of heavy atoms, such as those of importance in standard combustion mechanisms. To this end, we construct the CBH-ANL approach and illustrate that, for a set of 194 alkane oxidation species, it can be used to produce ΔH f (o K) values with 2σ uncertainties of 0.2-0.5 kcal mol -1 . This set includes the alkanes, hydroperoxides, and alkyl, peroxy, and hydroperoxyalkyl radicals for 17 representative hydrocarbon fuels containing up to 10 heavy atoms with various degrees of branching in the alkane backbone. The CBH-ANL approach, automated in the QTC and AutoMech software suites, builds balanced chemical equations for the calculation of ΔH f (o K), in which the reference species may be up to five heavy atoms. The high-level ANLo and ANL1 reference ΔH f (o K) values are further refined for even the largest of these reference species with a novel laddering approach. We perform a comprehensive quantification of the uncertainties for both the individual reference species (the largest of which is 0.15 kcal mol -1 ) and the propagation of those uncertainties when used in the calculation of ΔH f (o K) for the 194 target species. Here we examine the sensitivity of the predicted ΔH f (o K) values to (i) electronic energies from various methods, including ωB97X-D/cc-pVTZ, B2PLYP-D3/cc-pVTZ, CCSD(T)-F12b/cc-pVDZ-F12//B2PLYP-D3/cc-pVTZ, and CCSD(T)-F12b/cc-pVTZ-F12//B2PLYP-D3/cc-pVTZ; (ii) the zero-point vibrational energies (ZPVEs), where we consider harmonic ZPVEs as well as two scaling-based estimates of the anharmonic ZPVEs, all implemented for both ωB97X-D/cc-pVTZ and B2PLYP-D3/cc-pVTZ calculations; (iii) the particular CBH-ANL scheme employed; and (iv) the procedure for choosing the reference conformer for the analyses. The discussion concludes with a summary of the estimated overall uncertainty in the predictions and a validation of the predictions for the alkane subset.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

Convergent ab initio analysis of the multi-channel HOBr + H reaction

High-level potential energy surfaces for three reactions of hypobromous acid with atomic hydrogen were computed at the CCSDTQ/CBS//CCSDT(Q)/complete basis set level of theory. Focal point analysis was utilized to extrapolate energies and gradients for energetics and optimizations, respectively. The H attack at Br and subsequent Br–O cleavage were found to proceed barrierlessly. The slightly submerged transition state lies −0.2 kcal mol−1 lower in energy than the reactants and produces OH and HBr. The two other studied reaction paths are the radical substitution to produce H2O and Br with a 4.0 kcal mol−1 barrier and the abstraction at hydrogen to produce BrO and H2 with an 11.2 kcal mol−1 barrier. The final product energies lie −37.2, −67.9, and −7.3 kcal mol−1 lower in energy than reactants, HOBr + H, for the sets of products OH + HBr, H2O + Br, and H2 + BrO, respectively. Additive corrections computed for the final energetics, particularly the zero-point vibrational energies and spin–orbit corrections, significantly impacted the final stationary point energies, with corrections up to 6.2 kcal mol−1.

Chemistry↗
DOE OSTI2023

Isotope Effects in the Zundel–Eigen Isomerization of H + (H 2 O) 6

The isomerization pathway between the energetically low-lying Zundel and Eigen isomers of the protonated water hexamer was investigated using high-level ab initio calculations including a treatment of zero-point corrections. On the basis of these calculations, the Zundel–Eigen isomerization was found to proceed through a stable intermediate isomer, which consists of a four-membered ring with two single acceptor water molecules. Further, the inclusion of vibrational zero-point energy is shown to be important for accurately establishing the relative energies of the three relevant isomers involved in the Zundel–Eigen isomerization. Diffusion Monte Carlo calculations including anharmonic vibrational effects show that all three isomers of H + (H 2 O) 6 and D + (D 2 O) 6 have well-defined structures. The energetic ordering of the three isomers changes upon deuteration. The implications of these results for the vibrational spectra of H + (H 2 O) 6 and D + (D 2 O) 6 are also discussed.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

Transformation of an Embedded Five-Membered Ring in Polycyclic Aromatic Hydrocarbons via the Hydrogen-Abstraction–Acetylene-Addition Mechanism: A Theoretical Study

Five-membered rings are constituents of many polycyclic aromatic hydrocarbons (PAHs) and their presence on the edges of large PAHs has been repeatedly observed experimentally. However, modern kinetic combustion models often do not consider the growth of PAHs through the transformation of the five-membered rings. Here, in connection with the above, we carried out a theoretical study of the mechanism of Hydrogen-Abstraction-aCetylene-Addition (HACA) transformation of an embedded five-membered ring on the armchair PAH edge to a six-membered ring, considering cyclopenta[d,e,f]phenanthrene (4,5-methylenephenanthrene) as a prototype system for this process. The potential energy surface for the reactions of cyclopenta[d,e,f]phenanthrenyl radicals produced by direct H abstractions from cyclopenta[d,e,f]phenanthrene with acetylene has been compiled at the G3(MP2,CC)//B3LYP/6-311G(d,p) level of theory including zero-point vibrational energy corrections. The computed energies and molecular parameters were then used to solve the Rice-Ramsperger-Kassel-Marcus Master Equation (RRKM-ME) in order to calculate reaction rate at various pressures and temperatures, which were fitted to modified Arrhenius equation for further kinetic modeling. The results show that the HACA transformation of the embedded five-membered ring to a six-membered ring is possible, albeit slow. The most viable reaction mechanism involves the R2 + C 2 H 2 reaction, where the acetylene molecules adds to a σ-radical in the six-membered ring adjacent to the five-membered ring via a low entrance barrier. The predominant product of R2 + C 2 H 2 is predicted to be 3 ethynyl 4H cyclopenta[def]phenanthrene Pr5 via immediate H elimination from the initial addition complex. Next, Pr5 undergoes H-assisted isomerization to 4aH-pentaleno[4,3,2,1-cdef]phenanthrene Pr4, and the latter adds an H atom eventually forming the 1-pyrenylmethyl radical Pr3: R2 + C 2 H 2 ⇌ 3 ethynyl 4H cyclopenta[def]phenanthrene (Pr5) + H ⇌ 4aH-pentaleno[4,3,2,1-cdef]phenanthrene (Pr4) + H Pr5 + H ⇌ Pr4 + H Pr4 + H → 1-pyrenylmethyl (Pr3) This HACA sequence may be competitive with the methyl radical addition to the R1 radical formed by H abstraction from the CH 2 group in the five-membered ring of cyclopenta[d,e,f]phenanthrene, which provides a pathway to pyrene following two H atom losses. Relative contributions of the two mechanisms of the five- to six-membered ring transformation would strongly depend on the branching ratios of the R1 and R2 radicals produced by the H abstractions and the available concentration of C 2 H 2 vs. CH 3 and hence differ in different flames.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

The water trimer reaction OH + (H 2 O) 3 → (H 2 O) 2 OH + H 2 O

All important stationary points on the potential energy surface (PES) for the reaction OH + (H 2 O) 3 → (H 2 O) 2 OH + H 2 O have been fully optimized using the “gold standard” CCSD(T) method with the large Dunning correlation-consistent cc-pVQZ basis sets. Three types of pathways were found. For the pathway without hydrogen abstraction, the barrier height of the transition state (TS1) is predicted to lie 5.9 kcal mol -1 below the reactants. The two major complexes (H 2 O) 3 ∙∙∙OH (CP1 and CP2a) are found to lie 6.3 and 11.0 kcal mol -1 , respectively, below the reactants [OH + (H 2 O) 3 ]. For one of the H-abstraction pathways the lowest classical barrier height is predicted to be much higher, 6.1 kcal mol -1 (TS2a) above the reactants. For the other H-abstraction pathway the barrier height is even higher, 15.0 (TS3) kcal mol -1 . Vibrational frequencies and the zero-point vibrational energies connected to the PES are also reported. Fianlly, the energy barriers for the H-abstraction pathways are compared with those for the OH + (H 2 O) 2 and OH + H 2 O reactions, and the effects of the third water on the energetics are usually minor (0.2 kcal mol -1 ).

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

Energetic and Spectroscopic Properties of Astrophysically Relevant MgC 4 H Radicals Using High-Level Ab Initio Calculations

Considering the importance of magnesium-bearing hydrocarbon molecules (MgC n H; n = 2, 4, and 6) in the carbon-rich circumstellar envelopes (e.g., IRC+10216), a total of 28 constitutional isomers of MgC 4 H have been theoretically investigated using density functional theory (DFT) and coupled-cluster methods. The zero-point vibrational energy corrected relative energies at the ROCCSD(T)/cc-pCVTZ level of theory reveal that the linear isomer, 1-magnesapent-2,4-diyn-1-yl (1, 2 Σ + ), is the global minimum geometry on the MgC 4 H potential energy surface. The latter has been detected both in the laboratory and in the evolved carbon star, IRC+10216. The calculated spectroscopic data for 1 match well with the experimental observations (error ~ 0.78%) which validates our theoretical methodology. Plausible isomerization processes happening among different isomers are examined using DFT and coupled-cluster methods. CASPT2 calculations have been performed for a few isomers exhibiting multireference characteristics. The second most stable isomer, 1-ethynyl-1λ 3 -magnesacycloprop-2-ene-2,3-diyl (2, 2 A 1 , μ = 2.54 D), is 146 kJ mol –1 higher in energy than 1 and possibly the next promising candidate to be detected in the laboratory or in the interstellar medium in future.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

Artificial Neural Networks as Mappings between Proton Potentials, Wave Functions, Densities, and Energy Levels

We report artificial neural networks (ANNs) have become important in quantum chemistry. Herein, applications to nuclear quantum effects, such as zero-point energy, vibrationally excited states, and hydrogen tunneling, are explored. ANNs are used to solve the time-independent Schrödinger equation for single- and double-well potentials representing hydrogen-bonded molecular systems capable of proton transfer. ANN mappings are trained to predict the lowest five proton vibrational energies, wave functions, and densities from the proton potentials and to predict the excited state proton vibrational energies and densities from the proton ground state density. For the inverse problem, ANN mappings are trained to predict the proton potential from the proton vibrational energy levels or the proton ground state density. This latter mapping is theoretically justified by the first Hohenberg-Kohn theorem establishing a one-to-one correspondence between the external potential and the ground state density. ANNs for two- and three-dimensional systems are also presented to illustrate the straightforward extension to higher dimensions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

Eliminating Imaginary Vibrational Frequencies in Quantum-Chemical Cluster Models of Enzymatic Active Sites

In constructing finite models of enzyme active sites for quantum-chemical calculations, atoms at the periphery of the model must be constrained to stop unphysical rearrangements during geometry relaxation. A simple fixed-atom or “coordinate-lock” approach is commonly employed but leads to undesirable artifacts in the form of small imaginary frequencies. These preclude evaluation of finite-temperature free-energy corrections, limiting thermochemical calculations to enthalpies only. Full-dimensional vibrational frequency calculations are possible by replacing the fixed-atom constraints with harmonic confining potentials. Here, we compare that approach to an alternative strategy in which fixed-atom contributions to the Hessian are simply omitted. While the latter strategy does eliminate imaginary frequencies, it tends to underestimate both the zero-point energy and the vibrational entropy while introducing artificial rigidity. Harmonic confining potentials eliminate imaginary frequencies and provide a flexible means to construct active-site models that can be used in unconstrained geometry relaxations, affording better convergence of reaction energies and barrier heights with respect to the model size, as compared to models with fixed-atom constraints.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

Reliable p K a Prediction through Efficient Incorporation of Anharmonicity within the Nuclear–Electronic Orbital Framework

Accurate pK a prediction is critical for understanding chemical reactivity and molecular properties across a wide range of applications. Computational approaches usually invoke a harmonic treatment of the vibrational modes for zero-point energies, as well as thermal and entropic contributions. Herein, we present a general protocol for relative pK a prediction that incorporates the significant anharmonic effects using nuclear–electronic orbital (NEO) theory. This protocol is validated against experimental data for a range of molecules in acetonitrile, including protonated nitrogen bases, nitrophenols, anilines, and diamines, as well as cobalt electrocatalysts. For simple acids, the NEO approach offers only a slight improvement over conventional density functional theory with the standard harmonic vibrational treatment, whereas for hydrogen-bonded acids, the NEO approach offers more significantly improved performance at a comparable computational cost. This accessible methodology provides a practical route for accurate pKa prediction in challenging systems and is extendable to related thermodynamic properties such as hydricities and proton-coupled redox potentials.

Density functional theory↗
DOE OSTI2025

Anharmonic phonons, thermal expansion, and nuclear quantum effects in Zn

The anharmonic behavior of phonons and thermal expansion of hexagonal zinc were studied from 15 to 690 K by inelastic neutron scattering (INS) and ab initio simulations. Phonon spectra were measured for Q-points covering the full Brillouin zone, giving the phonon density of states (DOS) and dispersions along high-symmetry directions. The dispersions were sharp at 15 K, but diffuse intensity was observed at energies above them. The dispersions broadened with temperature T and the diffuse intensity became stronger. This diffuse intensity appeared in all INS measurements and simulations, except for classical molecular dynamics at 15 K. The temperature-dependent effective potential (TDEP) method, which included the nuclear quantum effect from zero-point vibrational dynamics, was used to calculate the free energy and thermal expansion. For T < 100 K nuclear quantum effects were important for obtaining the correct negative thermal expansion, and path integral molecular dynamics (PIMD) was particularly effective for obtaining the negative thermal expansion in the basal plane. A Heisenberg-Langevin model for interacting phonons coupled to a thermal bath was able to reproduce the shape and intensity of the diffuse spectral features.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

The Low-Lying Electronic States of NO 2 : Potential Energy and Dipole Surfaces, Bound States, and Electronic Absorption Spectrum

Nitrogen dioxide, NO 2 , is a free radical composed of the two most abundant elements in Earth’s atmosphere, nitrogen and oxygen, and is relevant to atmospheric and combustion chemistry. The electronic structure of even its lowest-lying states is remarkably complex, with various conical intersections and Renner–Teller pairings, giving rise to complex and perturbed vibronic states. Here we report some analysis of the 18 molecular states of doublet spin-multiplicity formed by combining ground-state N( 4 S u ) and O( 3 P g ) atoms. In this work, three-dimensional potential energy surfaces were fit at the MRCI(Q)-F12/VTZ-F12 level, describing the lowest four ($\tilde X$, $\tilde A$, $\tilde B$, and $\tilde C$) electronic states. A properties-based diabatization procedure was applied to accommodate the intersections, producing energies in a quasidiabatic representation and yielding couplings that were also fit into surfaces. The low-lying vibrational levels on the ground $\tilde X$ state were computed and compared with experimental measurements. Compared to experiment, the lowest 125 calculated vibrational levels (up to 8500 cm –1 above the zero-point energy) have a root-mean-squared error of 16.5 cm –1 . In addition, dipole moments for each of the lowest four electronic states—and the transition dipoles between them—were also computed and fit. With the coupled energy and dipole surfaces, the electronic spectrum was calculated in absolute intensity and compared with experimental measurements. Detailed structure in the experimental spectrum was successfully reproduced, and the total integrated intensity matches experiment to an accuracy of ~1.5% with no empirical adjustments.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

Nonlinear nanomechanical resonators approaching the quantum ground state

It is an open question whether mechanical resonators can be made nonlinear with vibrations approaching the quantum ground state. This requires the engineering of a mechanical nonlinearity far beyond what has been realized so far. Here we discover a mechanism to boost the Duffing nonlinearity by coupling the vibrations of a nanotube resonator to single-electron tunnelling and by operating the system in the ultrastrong-coupling regime. We find that thermal vibrations become highly nonlinear when lowering the temperature. The average vibration amplitude at the lowest temperature is 13 times the zero-point motion, with approximately 42% of the thermal energy stored in the anharmonic part of the potential. Our work may enable the realization of mechanical Schrödinger cat states, mechanical qubits and quantum simulators emulating the electron–phonon coupling.

36 MATERIALS SCIENCE↗
DOE OSTI2024

High-resolution infrared spectroscopy of jet cooled cyclobutyl in the α-CH stretch region: large-amplitude puckering dynamics in a 4-membered ring radical

We report gas-phase cyclobutyl radical (c-C 4 H 7 ) is generated at a rotational temperature of T rot = 26(1) K in a slit-jet discharge mixture of 70% Ne/30% He and 0.5–0.6% cyclobromobutane (c-C 4 H 7 Br). A fully rovibrationally resolved absorption spectrum of the α-CH stretch fundamental band between 3062.9 cm −1 to 3075.7 cm −1 is obtained and analyzed, yielding first precision structural and dynamical information for this novel radical species. The α-CH stretch band origin is determined to be 3068.7887(4) cm −1 , which implies only a modest (≈0.8 cm −1 ) blue shift from rotationally unresolved infrared spectroscopic studies of cyclobutyl radicals in liquid He droplets [ A. R. Brown, P. R. Franke and G. E. Douberly, J. Phys. Chem. A , 2017, 121 , 7576–7587]. Of particular dynamical interest, a one-dimensional potential energy surface with respect to the ring puckering coordinate is computed at CCSD(T)/ANO2 level of theory and reveals a double minimum C s puckered geometry, separated by an exceedingly shallow planar C 2v transition state barrier (E barr ≈ 1 cm −1 ). Numerical solutions on this double minimum potential yield a zero-point energy for the ground state (E zero-point ≈ 27 cm −1 ) greatly in excess of the interconversion barrier. This is indicative of highly delocalized, large amplitude motion of the four-membered ring structure, for which proper vibrationally averaging of the moment of inertia tensor reproduces the experimentally determined inertial defect remarkably well. Finally, intensity alternation in the experimental spectrum due to nuclear spin statistics upon exchange of three indistinguishable H atom pairs (IH = ½) matches K a + K c = even : odd = 36 : 28 predictions, implying that the unpaired electron in the radical center lies in an out-of-plane p π orbital. Thus, high-resolution infrared spectroscopy provides first experimental confirmation of a shallow double minimum ring puckering potential with a highly delocalized ground state wave function peaked at a planar C 2v transition state geometry consistent with a cyclobutyl π radical.

Chemistry↗
DOE OSTI2021

Spectroscopy and electronic structure of the low-energy states of ThN

The electronic spectrum of ThN over the range 19,600–21,200 cm -1 is remarkably congested, showing at least twenty vibronic bands originating from the ground state zero-point level. Rotationally resolved spectra for eleven of these bands have been examined to probe the underlying reasons for the dense manifold of states. Dispersed fluorescence spectra and fluorescence decay lifetimes were measured to provide addition insights. Our experimental measurements were complemented by electronic structure calculations for the low-energy doublet and quartet states. The ab initio calculations yielded a density of excited states that was consistent with the congested spectrum. These states were derived from the formal configurations Th 3+ (6d)N 3- , Th 2+ (7s 2 )N 2- (2p 5 ) and Th 2+ (7s6d)N 2 -(2p 5 ). The calculations indicated that the rotational and vibrational constants would be characteristic of the electronic configurations. This prediction was not borne out by the experimental data, suggesting that the full range of inter-configurational state interactions were not captured by the computational model.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗