Molecular States Research Articles

Spectroscopic properties under vibrational strong coupling in disordered matter from path-integral Monte Carlo simulations.

Vibrational strong coupling (VSC), the strong coupling between a Fabry-Perrot cavity and molecular vibrations at mid-infrared frequencies, has received important attention in the last years due to its capacity of modifying both vibrational spectra and chemical reactivity. VSC is a collective effect, and in this work, we introduce Path Integral Monte Carlo (PIMC) simulations that not only take into account the quantum character of the molecular vibrations and of the optical resonance of the cavity but also reproduce this collective behavior by considering multiple replicas of the molecular system. Moreover, we show that it is possible to extract from the PIMC simulations the decomposition of the hybrid optical and molecular states in terms of the bare molecular modes. On a model system of an ensemble of disordered Morse oscillators coupled to a single cavity through the Pauli-Fierz Hamiltonian, PIMC can retrieve known features obtained from analytical modes such as the Tavis-Cummings model and obtain a very close agreement with exact diagonalization for a small number of Morse oscillators. We also find that notwithstanding the anhamonic character of the Morse oscillators, the collective mode coupled to the cavity behaves as a harmonic oscillator, following the quantum central limit theorem.

TMOD-13. IDENTIFYING MOLECULAR SIGNATURES OF PPM1D-MUTANT DIPGS

Abstract Diffuse Intrinsic Pontine Glioma (DIPGs) are universally fatal pediatric tumors with a median survival of only 9-12 months. PPM1D-truncating mutations (PPM1Dtr), which promote its phosphatase activity and suppress TP53 function, are found in approximately 15% of DIPG tumors. PPM1Dtr are mutually exclusive to TP53 mutations, and present unique therapeutic opportunities to reactivate wild-type TP53. Another distinct feature of PPM1Dtr DIPGs is the co-occurrence of PI3K/AKT/mTOR pathway alterations, such as PIK3CA. There is an unmet need to elucidate the underlying molecular mechanisms of DIPG and identify novel therapeutic vulnerabilities that will improve outcomes of this devastating tumor. We hypothesize that PPM1D-mutant DIPGs can be therapeutically targeted by reactivating the p53 signaling pathway and other complementary pathways which cooperate with PPM1Dtr. To investigate this, we have generated Trp53- and Ppm1d-mutant isogenic DIPG murine models to determine how specific oncogenic alterations function in DIPG formation, growth and treatment sensitivity. Current studies aim to examine differences in sensitivity to standard of care treatment modalities and leverage next generation sequencing techniques to uncover additional vulnerabilities. We found notable differences upon addition of PIK3CA mutations in the context of PPM1Dtr DIPG tumor formation, growth, survival, and lineage identity. Ongoing studies using whole genome transcriptomics to compare tumor models will provide insight into the molecular signals and cell state changes driven by individual genetic alterations. Furthermore, by comparing isogenic TP53 and PPM1Dtr mutant DIPG models, we aim to further delineate the molecular mechanisms that drive differential sensitivity of TP53 wild-type tumors. The immune-competent nature of our murine models also provides us with a unique opportunity to delineate tumor immune microenvironment changes. Leveraging this information and our novel murine models, our long-term goal is to develop combinatorial therapeutic strategies to improve outcomes for this deadly disease.

DDDR-33. TRANSCRIPTOMICS-BASED NOVEL DRUG DISCOVERY IN DIFFUSE MIDLINE GLIOMA

Abstract Introduction. Pediatric high-grade gliomas are the most common cause of cancer-related death in children, of which the most malignant and devastating tumors include diffuse midline glioma (DMG). Recent availability of tumor samples and advances in next-generation sequencing enable the profiling of thousands of molecular features in DMG. We employ a system-based approach that uses gene expression as a representation of the molecular state, providing an avenue for therapeutic discovery. Previously, we used this approach to identify repurposed candidates for DMG; however, the final candidate had poor brain penetration, prompting us to discover novel compounds with better brain penetration. To apply it to novel compound discovery, we need the expression profiles of all library compounds, which is impractical for millions of compounds. To overcome this, we have developed a platform, Gene expression profiles Predictor on chemical Structures (GPS), that incorporates machine learning algorithms which leverage existing drug-induced gene expression profiles to predict gene expression based on compound structure. methods. DMG RNA-Seq samples were acquired from St. Jude Children’s Research Hospital and Children’s Brain Tumor Network (CBTN). We developed deep learning autoencoder to search reference normal tissues from Genotype-Tissue Expression (GTEx). A DMG meta-signature was created and fed into GPS to predict drug candidates from the Enamine CNS library, seven of which were experimentally tested in vitro and the most promising one subsequently validated in vivo. Results. Cell proliferation analyzed by MTS assay showed three compounds with IC50< 50 µM in K27M-mutant DMG cell lines (SF8628, DIPG007, SU-DIPG36). Toxicity study showed no adverse effects in the mice treated with 100-200 mg/kg of the most promising compound (IC50 = 1.9 µM in DMG cells). Conclusion. Our novel computational framework leverages deep learning to fill the gap in drug discovery in DMG. Work is underway to determine antitumor activity, brain penetration (via HPLC) and survival benefit in the mice bearing DMG PDX treated with the lead compound.

Production mechanism of the hidden charm pentaquark states Pcc¯

We investigate hidden charm pentaquark states using an off-shell coupled-channel formalism involving heavy meson and singly heavy baryon scattering. Our approach utilizes an effective Lagrangian to construct the kernel amplitudes, which respect both heavy quark symmetry and hidden local symmetry. After solving the coupled integral equations, we obtain the transition amplitudes for J/ψN scattering and various heavy meson and singly heavy baryon scattering processes. We identify seven distinct peaks related to molecular states of heavy mesons D¯ (D¯*) and singly heavy baryons Σc (Σc*). Four of these peaks can be associated with the known Pcc¯ states: Pcc¯(4312), Pcc¯(4380), Pcc¯(4440), and Pcc¯(4457). We predict two additional resonances with masses around 4.5 GeV, which we interpret as D¯*Σc* molecular states, and identify one cusp structure. Additionally, we predict two P-wave pentaquark states with positive parity, which may be candidates for genuine pentaquark configurations. Notably, these pentaquark states undergo significant modifications in the J/ψN elastic channel, with some even disappearing due to interference from the positive parity channel. The present investigation may provide insight into the absence of pentaquark states in J/ψ photoproduction observed by the GlueX Collaboration. Published by the American Physical Society 2024

Cross Peaks on Two-Dimensional Optical Spectra Arising from Quantum Cross-Correlation Functions.

Cross peaks on 2D optical spectra are indicative of interactions between molecular excitonic states. Currently, the two conventional assignments of cross peaks are direct coupling and population transfer between excitonic states. Here, we show that there is another possible source of cross peaks. We theoretically demonstrate that for a model comprising two nondirectly interacting excitons or two-level systems (TLSs), cross peaks can arise if there is a complex-valued or quantum frequency-gap cross-correlation function between the two TLSs. Considering only real-valued or classical cross-correlation functions will result in no cross peaks. We derive and validate the mathematical expressions describing such cross peaks. We then simulate the 2D electronic spectra of an example model system comprising nondirectly interacting TLSs whose quantum cross-correlation functions arise from coupling to a common overdamped Brownian oscillator mode. We show that there are clear observational differences between such quantum correlation cross peaks with conventional direct coupling and population transfer cross peaks.

Ab initio calculations of molecular double Auger decay rates.

We report on the application of the recently developed Fano-ADC(2,2) method to compute total and partial Auger decay widths of molecular core-hole states, including explicit evaluation of double Auger decay branching ratios. The method utilizes the fast-convergent intermediate state representation to construct many-electron wave functions and is readily applicable to atoms, molecules, and clusters. The ADC(2,2) scheme describes the initial and final states of the normal Auger decay consistently up to the second order of perturbation theory. In addition, excitations with two electrons in the continuum provide access to three-electron decay modes. The method yields decay widths and the Auger electron spectra in excellent agreement with the experiment, demonstrating the high accuracy of partial widths. The average relative error of double Auger decay branching ratios compared to available experimental data is about 30%, which should be evaluated as an excellent result considering the universality of the method, the complexity of the double decay process, and the neglection of nuclear motion in this study.

The Perils of Molecular Interpretations from Vibrational Spectra of Complex Samples.

Vibrational spectroscopy is a widely used technique for chemical characterizations across various analytical sciences. Its applications are increasingly extending to the analysis of complex samples such as biofluids, providing high-throughput molecular profiling. While powerful, the technique suffers from an inherent limitation: The overlap of absorption information across different spectral domains hinders the capacity to identify individual molecular substances contributing to measured signals. Despite the awareness of this challenge, the difficulty of analyzing multi-molecular spectra is often underestimated, leading to unsubstantiated molecular interpretations. Here, we examine the prevalent overreliance on spectral band assignment and illuminate the pitfalls of correlating spectral signals to discrete molecular entities or physiological states without rigorous validation. Focusing on blood-based infrared spectroscopy, we provide examples illustrating how peak overlap among different substances, relative substance concentrations, and preprocessing steps can lead to erroneous interpretations. We advocate for a viewpoint shift towards a more careful understanding of complex spectra, which shall lead to either accepting their fingerprinting nature and leveraging machine learning analysis - or involving additional measurement modalities for robust molecular interpretations. Aiming to help translate and improve analytical practices within the field, we highlight the limitations of molecular interpretations and feature their viable applications.

Sensitive Self-Driven Single-Component Organic Photodetector Based on Vapor-Deposited Small Molecules.

Typically, organic solar cells (OSCs) and photodetectors (OPDs) comprise an electron donating and accepting material to facilitate efficient charge carrier generation. This approach has proven successful in achieving high-performance devices but has several drawbacks for upscaling and stability. This study presents a fully vacuum-deposited single-component OPD, employing the neat oligothiophene derivative DCV2-5T in the photoactive layer. Free charge carriers are generated with an internal quantum efficiency of 20 % at zero bias. By optimizing the device structure, a very low dark current of 3.4 · 10-11Acm-2 at -0.1V is achieved, comparable to the dark current of state-of-the-art bulk heterojunction OPDs. This optimization results in specific detectivities of 1· 1013Jones (based on noise measurements), accompanied by a fast photoresponse (f-3dB = 200kHz) and a broad linear dynamic range (>150dB). Ultrafast transient absorption spectroscopy unveils that charge carriers are already formed at very short time scales (<1ps). The surprisingly efficient bulk charge generation mechanism is attributed to a strong electronic coupling of the molecular exciton and charge transfer states. This work demonstrates the very high performance of single-component OPDs and proves that this novel device design is a successful strategy for highly efficient, morphological stable and easily manufacturable devices.

Molecular design of furanoacene-based singlet fission sensitizers combining diradical character and cross-substitution

Singlet fission (SF) suffers from the scarcity of available SF chromophores and sensitizers. Tetracene, pentacene and their derivatives, are two typical model systems for the study of SF photophysical phenomena. Herein, the analogues of polyacenes, furanoacenes and their derivatives, were evaluated as potential SF chromophores and sensitizers through a theoretical study. The primary photostability were evaluated by the frontier molecular orbital energy levels, diradical characters and SF relevant excited state energies according to the type-I and −II photodegradation mechanisms. The cross-substitutions by combining central substitutions with triple bonds and terminal substituents enhance photostability, tune the SF relevant excited states efficiently, and give more appropriate E(T1) to be SF sensitizers. This work helps to give better understanding of the electronic structures and SF capability of furanoacenes.

Predicting hidden bottom molecular tetraquarks with a complex scaling method

In the present work, we perform a coupled-channel analysis of B(s)(*)B¯(s)(*) systems with the one-boson-exchange potentials. We first study the I(JPC)=1(1+−) BB¯*/B*B¯* system to describe the Zb(10610) and Zb(10650) particles as molecular states and determine the reasonable range of cutoff parameter Λ. Then, other B(s)(*)B¯(s)(*) combinations with different quantum numbers are systematically investigated. Some bound states and resonances appear in the isoscalar systems, while only several shallow bound states exist for isovector systems. Far away from the excited conventional P-wave bottomium, these predicted states can be easily identified as exotic particles both theoretically and experimentally. Moreover, the ηb(nS)/ϒ(nS) plus light mesons are the excellent final states to search for the bound states, while the BB¯*+H.c. and B*B¯* channels are suitable for the resonances. We highly recommend that the LHCb and Belle II Collaborations hunt for these bottomoniumlike states in the future. Published by the American Physical Society 2024

Structural Manipulation of Spin Excitations in a Molecular Junction.

Single metallocene molecules act as sensitive spin detectors when decorating the probe of a scanning tunneling microscope (STM). However, the impact of the atomic-scale electrode details on the molecular spin state has remained elusive to date. Here, a nickelocene (Nc) STM junction is manipulated in an atomwise manner showing clearly the dependence of the spin excitation spectrum on the anchoring of Nc to Cu(111), a Cu monomer, and trimer. Moreover, while the spin state of the same Nc tip is a triplet with tunable spin excitation energies upon contacting the surface, it transitions to a Kondo-screened doublet on a Cu atom. Notably, the nontrivial magnetic exchange interaction of the molecular spin with the electron continuum of the substrate determines the spectral line shape of the spin excitations.

Artificial-intelligence-driven shot reduction in quantum measurement

Variational Quantum Eigensolver (VQE) provides a powerful solution for approximating molecular ground state energies by combining quantum circuits and classical computers. However, estimating probabilistic outcomes on quantum hardware requires repeated measurements (shots), incurring significant costs as accuracy increases. Optimizing shot allocation is thus critical for improving the efficiency of VQE. Current strategies rely heavily on hand-crafted heuristics requiring extensive expert knowledge. This paper proposes a reinforcement learning (RL)-based approach that automatically learns shot assignment policies to minimize total measurement shots while achieving convergence to the minimum of the energy expectation in VQE. The RL agent assigns measurement shots across VQE optimization iterations based on the progress of the optimization. This approach reduces VQE's dependence on static heuristics and human expertise. When the RL-enabled VQE is applied to a small molecule, a shot reduction policy is learned. The policy demonstrates transferability across systems and compatibility with other wavefunction Ansätze. In addition to these specific findings, this work highlights the potential of RL for automatically discovering efficient and scalable quantum optimization strategies.

Analysis of DD* and D¯(*)Ξcc(*) molecule by one boson exchange model based on heavy quark symmetry

Numerous exotic hadrons with heavy quarks have been reported in the experiments. In such states, symmetries of heavy quarks are considered to play a significant role. In particular, the superflavor symmetry, or also called the heavy quark antidiquark symmetry is one of the interesting symmetries, which links a quark Q to an antidiquark Q¯Q¯ having the same color representation as Q. In this paper, we study a D¯Ξcc molecular state as a superflavor partner of the doubly charmed tetraquark Tcc reported by LHCb recently. Tcc locating slightly below the DD* threshold is a candidate of the hadronic molecule. Thus by replacing the singly charmed meson D(*) with the doubly charmed baryon Ξcc(*), superflavor symmetry predicts the existence of the D¯(*)Ξcc(*) bound state. We employ the one boson exchange model respecting with the heavy quark spin symmetry where the parameter is obtained to reproduce the Tcc binding energy. We apply this model with superflavor symmetry to the D¯(*)Ξcc(*) molecule and predict a bound state with I(JP)=0(12−). If the pentaquark state corresponding to D¯(*)Ξcc(*) molecular state is observed in future experiments as predicted in this work, it is more likely that Tcc is a DD* molecular state. Published by the American Physical Society 2024

Synergistically Activated Aggregation-Induced Emission Probe for Precise In Situ Staining of Lipids in Atherosclerotic Plaques.

Visualizing the localization and distribution of lipids within arteries is crucial for studying atherosclerosis. However, existing lipid-specific probes face challenges such as strong hydrophobicity and nonspecific staining of lipophilic organelles or tissues, making them impractical for the precise identification of atherosclerotic plaques. To address this issue, we design a synergistically activated probe, Cbz-Lys-Lys-TPEB, which responds to cathepsin B (CTB) and H2O2 for the in situ generation of aggregation-induced emission luminogens (AIEgens). This enables specific staining of lipids within arteries and precise imaging of atherosclerotic plaques. The probe combines a tetraphenylethene building block with a hydrophilic peptide sequence (Cbz-Lys-Lys) and phenylboric acid module, providing excellent water solubility and fluorescence quenching in a molecular dispersion state. Upon interaction with H2O2 and CTB within plaques, the hydrophilic Cbz-Lys-Lys-TPEB probe is specifically cleaved and converted into hydrophobic AIEgens, leading to rapid aggregation and significant fluorescence enhancement. Interestingly, the in situ-liberated AIEgens display distinct lipid binding ability, effectively tracking the location and distribution of lipids in plaques. This synergistic target-activated AIEgen liberation strategy demonstrates significant feasibility for the reliable and accurate identification of atherosclerotic plaques, holding tremendous potential for clinical diagnosis and risk stratification of atherosclerosis.

Sensing the Small Change of Intermolecular Distance in Supramolecular Assembly by Using the Tunable Emission Wavelength of AIE-Active Luminogens.

The distinct molecular states - single molecule, assembly, and aggregate - of two ionic macromolecules, TPPE-APOSS and TPE-APOSS, are easily distinguishable through their tunable fluorescence emission wavelengths, which reflect variations in intermolecular distances. Both ionic macromolecules contain aggregation-induced emission (AIE) active moieties whose emission wavelengths are directly correlated to their mutual distances in solution: far away from each other as individual molecules, maintaining a tunable and relatively long distance in electrostatic interactions-controlled blackberry-type assemblies (microphase separation), or approaching van der Waals close distance in aggregates (macrophase separation). Furthermore, within the blackberry assemblies, the emission wavelength decreases monotonically with increasing assembly size, indicative of shorter intermolecular distances at nanoscale. The emission changes of TPPE-APOSS blackberry assemblies can even be visually distinguishable by eyes when their sizes and intermolecular distances are tuned. Molecular dynamics simulations further revealed that macromolecules are confined in various conformations by controllable intermolecular distances within the blackberry structure, thereby resulting in fluorescence emission with tunable wavelength.

Nonadiabatic dynamical simulations to the radiative recombination of nonfullerene acceptor molecular excited states.

Improving the radiative recombination rate of nonfullerene acceptor (NFA) molecular excited states can help to promote their photoluminescence quantum yield and thus reduce the nonradiative energy loss in NFA-based organic solar cells. In this Letter, by developing a nonadiabatic dynamical simulation method, we clarify quantitative correlations of some typical characteristics of NFA molecules with their radiative recombination rates. For a single NFA molecule, the weakening of electron-phonon coupling and the strengthening of electron-push-pull potential can each improve the radiative recombination rate. For different NFA molecular aggregates, their radiative recombination rates are all reduced compared with a single molecule, where the A-to-A and A-to-D type J-aggregates have higher rates than D-to-D type H-aggregate. To further improve the radiative recombination rate of NFA molecular J-aggregates, we should increase the intermolecular distance, such as extending the side chain length.

Phenol Photostatic Spectra and Quantum‐Classical Photodynamic Deprotonation

ABSTRACTThe spectral‐luminescence properties and photochemical conversions of phenol were analyzed for an isolated molecule as well as in water solvents in a continuum implicit model and explicit atomistic surroundings. This involved employing cut‐edge hybrid quantum‐classical methodologies to generate static optical spectra and the excited dissipative crossing potential energy curves. A combination of electronic excitations, gradient calculations, and embedding electrostatic potential fitting charges on quantum‐classical molecular dynamic propagation trajectories provided statistically averaged absorption spectra. The mixed‐reference spin‐flip multiconfigurational linear response method based on reference triplet preprocessed in the time‐dependent density‐functional theory was utilized to determine conical intersections between the lowest excited and ground states, as well as two‐stage transitions from the second excitation to the ground state. Non‐adiabatic quantum‐classical molecular dynamics defined photodissipative trajectories of excited states, their lifetimes, and crossing points through trajectory surface hopping together with the mixed‐reference spin‐flip and embedding electrostatic potential fitting approaches. Dyson orbitals of the extended Koopmans' theorem were applied to reveal the nature of molecular states at conical intersections and key points on photodynamic trajectories. Potential hydroxyl group cleavage predicted with conical intersections searching turns to “swift” OH deprotonation through |π→⟩ transition along photodynamic propagations in contrast with “long” processes leading to benzene ring deformation with stable OH bond.

Advancing Optoelectronics with Benzonitriles: Theoretical Understanding, Experimental Realization, and Single‐Molecule White Light Emission

AbstractBenzonitrile derivatives (BDs) are very promising for applications in materials science, photonics, and optoelectronics due to their intriguing electronic and optical properties. This study comprehensively investigates BDs, aiming to uncover their fundamental characteristics and potential applications. Using advanced theoretical techniques like Gaussian software, it is gained unprecedented insights into the photoinduced isomerization phenomenon for all‐optical switching in these compounds. The theoretical framework clarifies molecular transition states and explores a range of properties, providing a comprehensive understanding of BDs. Empirical data on the emission properties of BDs, from fluorescence analyses in liquid solutions to light amplification in solid‐state PMMA thin films, complement the theoretical examinations. Notably, white light emission from a single benzonitrile compound is achieved, showcasing its potential for data transmission through Li‐Fi technology. Finally, the all‐optical switching phenomenon in BDs using 3rd order nonlinear optical effects is experimentally validated. This complementary and comprehensive study advances understanding of BDs and demonstrates their potential for practical applications in emerging photonic and optoelectronic technologies.

Maximum entropy method for valence quark distributions in exotic hadrons: a study of the Zc(3900)\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$Z_c(3900)$$\\end{document} case

In this study we demonstrate the application of the maximum entropy method (MEM) to determine the valence quark distribution of exotic hadrons. Our investigation yields three key findings. Firstly, we observe a significant shift towards smaller Bjorken scale x in the peak position of the valence quark distribution for hadrons with an increasing number of valence quarks, consistent with previous results by Kawamura and Kumano. Secondly, assuming that the Zc(3900)\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$Z_c(3900)$$\\end{document} initially consists of four valence quarks, we employ MEM to determine its initial valence quark distribution, estimating a radius of rc=1.276\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$r_c=1.276$$\\end{document} fm at an extremely low resolution scale Q2\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$Q^2$$\\end{document}. Furthermore, we identify a notable discrepancy between our computed charge form factor Gc(q)\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$G_c(q)$$\\end{document} at leading order and the outcomes of hadron molecular state calculations. We propose that this form factor can be extracted from the QCD counting rule cross-section, which is grounded in generalized distribution amplitudes (GDA) linked to the multi-quark states.

Spin polarization and transport of hybrid interface states in π-stacked magnetic molecular junctions

Based on the first-principles method, the spin polarization and transport properties of the hybrid interface states (HIS) in π-stacked three-layer triphenylene molecular junction are studied theoretically. By rotating the torsion angle between two neighbor layers, the π-π interaction strength is adjusted. The results demonstrate that the HIS, generated by strong pz-d hybridization between molecules and ferromagnetic electrodes, can easily transfer to the central layer via a strong π-π interaction with similar spin polarization. A weak π-π interaction will block the transfer of the HIS, where the discrete molecular orbital states dominate the central molecule. However, an enhancement of the spin polarization for the molecular frontier orbital states is observed. The transport calculation shows that efficient transport of the HIS is realized in the π-stacked system, where the largest current appears at a strong π-π interaction but a significant current spin polarization appears at a weak π-π interaction. This work deepens our understanding of the role of the HIS in π-stacked organic systems.