Comparative molecular dynamics simulations provided insights into the mechanisms of cold-adaption of alginate lyases from the PL7 family.

Alginate is an important polysaccharide that is abundant in the marine environments, including the Polar Regions, and bacterial alginate lyases play key roles in its degradation. Many reported alginate lyases show characteristics of cold-adapted enzymes, including relatively low temperature optimum of activities (Topt) and low thermal stabilities. However, the cold-adaption mechanisms of alginate lyases remain unclear. Here, we studied the cold-adaptation mechanisms of alginate lyases by comparing four members of the PL7 family from different environments: AlyC3 from the Arctic ocean (Psychromonas sp. C-3), AlyA1 from the temperate ocean (Zobellia galactanivorans), PA1167 from the human pathogen (Pseudomonas aeruginosa PAO1), and AlyQ from the tropic ocean (Persicobacter sp. CCB-QB2). Sequence comparison and comparative molecular dynamics (MD) simulations revealed two main strategies of cold adaptation. First, the Arctic AlyC3 and temperate AlyA1 increased the flexibility of the loops close to the catalytic center by introducing insertions at these loops. Second, the Arctic AlyC3 increased the electrostatic attractions with the negatively charged substrate by introducing a high portion of positively charged lysine at three of the insertions mentioned above. Furthermore, our study also revealed that the root mean square fluctuation (RMSF) increased greatly when the temperature was increased to Topt or higher, suggesting the RMSF increase temperature as a potential indicator of the cold adaptation level of the PL7 family. This study provided new insights into the cold-adaptation mechanisms of bacterial alginate lyases and the marine carbon cycling at low temperatures.

Annotation of 2,507 Saccharomyces cerevisiae genomes.

Saccharomyces cerevisiae (baker's yeast, budding yeast) is one of the most important model organisms for biological research and is a crucial microorganism in industry. Currently, a huge number of Saccharomyces cerevisiae genome sequences are available at the public domain. However, these genomes are distributed at different websites and a large number of them are released without annotation information. To provide one complete annotated genome data resource, we collected 2,507 Saccharomyces cerevisiae genome assemblies and re-annotated 2,506 assemblies using a custom annotation pipeline, producing a total of 15,407,164 protein-coding gene models. With a custom pipeline, all these gene sequences were clustered into families. A total of 1,506 single-copy genes were selected as marker genes, which were then used to evaluate the genome completeness and base qualities of all assemblies. Pangenomic analyses were performed based on a selected subset of 847 medium-high-quality genomes. Statistical comparisons revealed a number of gene families showing copy number variations among different organism sources. To the authors' knowledge, this study represents the largest genome annotation project of S. cerevisiae so far, providing rich genomic resources for the future studies of the model organism S. cerevisiae and its relatives.IMPORTANCESaccharomyces cerevisiae (baker's yeast, budding yeast) is one of the most important model organisms for biological research and is a crucial microorganism in industry. Though a huge number of Saccharomyces cerevisiae genome sequences are available at the public domain, these genomes are distributed at different websites and most are released without annotation, hindering the efficient reuse of these genome resources. Here, we collected 2,507 genomes for Saccharomyces cerevisiae, performed genome annotation, and evaluated the genome qualities. All the obtained data have been deposited at public repositories and are freely accessible to the community. This study represents the largest genome annotation project of S. cerevisiae so far, providing one complete annotated genome data set for S. cerevisiae, an important workhorse for fundamental biology, biotechnology, and industry.

Mechanistic photophysics of tellurium-substituted cytosine: Electronic structure calculations and nonadiabatic dynamics simulations.

Previously, the MS-CASPT2 method was performed to study the static and qualitative photophysics of tellurium-substituted cytosine (TeC). To get quantitative information, we used our recently developed QTMF-FSSH dynamics method to simulate the excited-state decay of TeC. The CASSCF method was adopted to reduce the calculation costs, which was confirmed to provide reliable structures and energies as those of MS-CASPT2. A detailed structural analysis showed that only 5% trajectories will hop to the lower triplet or singlet state via the twisted (S2 /S1 /T2 )T intersection, while 67% trajectories will choose the planar intersections of (S2 /S1 /T3 /T2 /T1 )P and (S2 /S1 /T2 /T1 )P but subsequently become twisted in other electronic states. By contrast, ~28% trajectories will maintain in a plane throughout dynamics. Electronic population revealed that the S2 population will ultrafast transfer to the lower triplet or singlet state. Later, the TeC system will populate in the spin-mixed electronic states composed of S1 , T1 and T2 . At the end of 300 fs, most trajectories (~74%) will decay to the ground state and only 17.4% will survive in the triplet states. Our dynamics simulation verified that tellurium substitution will enhance the intersystem crossings, but the very short triplet lifetime (ca. 125 fs) will make TeC a less effective photosensitizer.

Insights into the composition and assembly mechanism of microbial communities on intertidal microsand grains.

Introduction: Marine microorganisms are essential in marine ecosystems and have always been of interest. Currently, most marine microbial communities are studied at the bulk scale (millimeters to centimeters), and the composition, function and underlying assembly mechanism of microbial communities at the microscale (sub-100 micrometers) are unclear. Methods: The microbial communities on microsand grains (40-100 µm, n = 150) from marine sediment were investigated and compared with those on macrosand grains (400-1000 µm, n = 60) and bulk sediments (n = 5) using amplicon sequencing technology. Results: The results revealed a significant difference between microsand grains and macrosand grains. Microsand grains had lower numbers of operational taxonomic units (OTUs(97%)) and predicted functional genes than macrosand grains and bulk-scale samples. Microsand grains also showed greater intersample differences in the community composition and predicted functional genes than macrosand grains, suggesting a high level of heterogeneity of microbial communities at the microscale. Analyses based on ecological models indicated that stochastic processes dominated the assembly of microbial communities on sand grains. Consistently, cooccurrence network analyses showed that most microbial cooccurrence associations on sand grains were highly unstable. Metagenomic sequencing and further genome-scale metabolic modeling revealed that only a small number (1.3%) of microbe pairs showed high cooperative potential. Discussion: This study explored the microbial community of marine sediments at the sub-100 µm scale, broadening the knowledge of the structure and assembly mechanism of marine microbial communities.

Enrichment Pretreatment Expands the Microbial Diversity Cultivated from Marine Sediments.

The majority of the microbial diversity in nature has not been recovered through cultivation. Enrichment is a classical technique widely used in the selective cultivation of specific taxa. Whether enrichment is suitable for cultivation studies that aim to recover large numbers of species remains little explored. To address this issue, we evaluated the potential of enrichment pretreatment in the cultivation of bacteria from marine sediments. Upon obtaining and classifying a total of 943 pure cultures from chitin and cellulose enrichment pretreatment systems and a control system, our results showed that species obtained using enrichment pretreatment differed greatly from those without enrichment. Multiple enrichment media and different enrichment times increased the number of cultivated species in a sample. Amplicon sequencing showed that the increased relative abundance during pretreatment contributed greatly to bacterial cultivation. The testing of degradation abilities against chitin and cellulose and the whole-genome sequencing of representative strains suggested that microorganism-microorganism interactions play roles in the expanded diversity of cultivated bacteria. This study provides new insights into the abilities of enrichment in exploring cultivable diversity and mining microbial resources.

Quantum mechanics/molecular mechanics studies on the excited-state decay mechanisms of cytidine aza-analogues: 5-azacytidine and 2'-deoxy-5-azacytidine in aqueous solution.

The excited state properties and deactivation pathways of two DNA methylation inhibitors, i.e., 5-azacytidine (5ACyd) and 2'-deoxy-5-azacytidine (5AdCyd) in aqueous solution, are comprehensively explored with the QM(CASPT2//CASSCF)/MM protocol. We systematically map the feasible decay mechanisms based on the obtained excited-state decay paths involving all the identified minimum-energy structures, conical intersections, and crossing points driving the different internal conversion (IC) and intersystem crossing (ISC) routes in and between the 1ππ*, 1nπ*, 3ππ*, 3nπ*, and S0 states. Unlike the 1nπ* state below the 1ππ* state in 5ACyd, deoxyribose group substitution at the N1 position leads to the 1ππ* state becoming the S1 state in 5AdCyd. In 5ACyd and 5AdCyd, the initially populated 1ππ* state mainly deactivates to the S0 state through the direct 1ππ* → S0 IC or mediated by the 1nπ* state. The former nearly barrierless IC channel of 1ππ* → S0 occurs ultrafast via the nearby low-lying 1ππ*/S0 conical intersection. In the latter IC channel of 1ππ* → 1nπ* → S0, the initially photoexcited 1ππ* state first approaches the nearby S2/S1 conical section 1ππ*/1nπ* and then undergoes efficient IC to the 1nπ* state, followed by the further IC to the initial S0 state via the S1/S0 conical intersection 1nπ*/S0. The 1nπ*/S0 conical intersection is estimated to be located 6.0 and 4.9 kcal mol-1 above the 1nπ* state minimum in 5ACyd and 5AdCyd, respectively, at the QM(CASPT2)/MM level. In addition to the efficient singlet-mediated IC channels, the minor ISC routes would populate 1ππ* to T1(ππ*) through 1ππ* → T1 or 1ππ* → 1nπ* → T1. Relatively, the 1ππ* → 1nπ* → T1 route benefits from the spin-orbit coupling (SOC) of 1nπ*/3ππ* of 8.7 cm-1 in 5ACyd and 10.2 cm-1 in 5AdCyd, respectively. Subsequently, the T1 system will approach the nearby T1/S0 crossing point 3ππ*/S0 driving it back to the S0 state. Given the 3ππ*/S0 crossing point located above the T1 minimum and the small T1/S0 SOC, i.e., 8.4 kcal mol-1 and 2.1 cm-1 in 5ACyd and 6.8 kcal mol-1 and 1.9 cm-1 in 5AdCyd, respectively, the slow T1 → S0 would trap the system in the T1 state for a while. The present work could contribute to understanding the mechanistic photophysics and photochemistry of similar aza-nucleosides and their derivatives.

Improved Interfacial Ion Migration and Deposition through the Chain-Liquid Synergistic Effect by a Carboxylated Hydrogel Electrolyte for Stable Zinc Metal Anodes.

The large-scale applicability of Zn-metal anodes is severely impeded by the issues such as the dendrite growth, complicated hydrogen evolution, and uncontrollable passivation reaction. Herein, a negatively charged carboxylated double-network hydrogel electrolyte (Gelatin/Sodium alginate-acetate, denoted as Gel/SA-acetate) has been developed to stabilize the interfacial electrochemistry, which restructures a type of Zn2+ ion solvent sheath optimized via a chain-liquid synergistic effect. New hydrogen bonds are reconstructed with water molecules by the zincophilic functional groups, and directional migration of hydrated Zn2+ ions is therefore induced. Concomitantly, the robust chemical bonding of such hydrogel layers to the Zn slab exhibits a desirable anti-catalytic effect, thereby greatly diminishing the water activity and eliminating side reactions. Subsequently, a symmetric cell using the Gel/SA-acetate electrolyte demonstrates a reversible plating/stripping performance for 1580 h, and an asymmetric cell reaches a state-of-the-art runtime of 5600 h with a high average Coulombic efficiency of 99.9 %. The resultant zinc ion hybrid capacitors deliver exceptional properties including the capacity retention of 98.5 % over 15000 cycles, energy density of 236.8 Wh kg-1 , and high mechanical adaptability. This work is expected to pave a new avenue for the development of novel hydrogel electrolytes towards safe and stable Zn anodes.

Understanding the Excited-State Relaxation Mechanisms of Xanthophyll Lutein by Multi-configurational Electronic Structure Calculations.

The contradictory behaviors in light harvesting and non-photochemical quenching make xanthophyll lutein the most attractive functional molecule in photosynthesis. Despite several theoretical simulations on the spectral properties and excited-state dynamics, the atomic-level photophysical mechanisms need to be further studied and established, especially for an accurate description of geometric and electronic structures of conical intersections for the lowest several electronic states of lutein. In the present work, semiempirical OM2/MRCI and multi-configurational restricted active space self-consistent field methods were performed to optimize the minima and conical intersections in and between the 1Ag-, 2Ag-, 1Bu+, and 1Bu- states. Meanwhile, the relative energies were refined by MS-CASPT2(10,8)/6-31G*, which can reproduce correct electronic state properties as those in the spectroscopic experiments. Based on the above calculation results, we proposed a possible excited-state relaxation mechanism for lutein from its initially populated 1Bu+ state. Once excited to the optically bright 1Bu+ state, the system will propagate along the key reaction coordinate, i.e., the stretching vibration of the conjugated carbon chain. During this period of time, the 1Bu- state will participate in and forms a resonance state between the 1Bu- and 1Bu+ states. Later, the system will rapidly hop to the 2Ag- state via the 1Bu+/2Ag- conical intersection. Finally, the lutein molecule will survive in the 2Ag- state for a relatively long time before it internally converts to the ground state directly or via a twisted S1/S0 conical intersection. Notably, though the photophysical picture may be very different in solvents and proteins, the current theoretical study proposed a promising calculation protocol and also provided many valuable mechanistic insights for lutein and similar carotenoids.

Phosphorus-Mediated Local Charge Distribution of N-Configuration Adsorption Sites with Enhanced Zincophilicity and Hydrophilicity for High-Energy-Density Zn-Ion Hybrid Supercapacitors.

Tailor-made carbonaceous-based cathodes with zincophilicity and hydrophilicity are highly desirable for Zn-ion storage applications, but it remains a great challenge to achieve both advantages in the synthesis. In this work, a template electrospinning strategy is developed to synthesize nitrogen and phosphorous co-doped hollow porous carbon nanofibers (N, P-HPCNFs), which deliver a high capacity of 230.7 mAh g-1 at 0.2 A g-1 , superior rate capability of 131.0 mAh g-1 at 20 A g-1 , and a maximum energy density of 196.10 Wh kg-1 at the power density of 155.53 W kg-1 . Density functional theory calculations (DFT) reveal that the introduced P dopants regulate the distribution of local charge density of carbon materials and therefore facilitate the adsorption of Zn ions due to the increased electronegativity of pyridinic-N. Ab initio molecular dynamics (AIMD) simulations indicate that the doped P species induce a series of polar sites and create a hydrophilic microenvironment, which decreases the impedance between the electrode and the electrolyte and therefore accelerates the reaction kinetics. The marriage of ex situ/in situ experimental analyses and theoretical simulations uncovers the origin of the enhanced zincophilicity and hydrophilicity of N, P-HPCNFs for energy storage, which accounts for the faster ion migration and electrochemical processes.

Quantum mechanics/molecular mechanics studies on mechanistic photophysics of cytosine aza-analogues: 2,4-diamino-1,3,5-triazine and 2-amino-1,3,5-triazine in aqueous solution.

The excited-state properties and photophysics of cytosine aza-analogues, i.e., 2,4-diamino-1,3,5-triazine (2,4-DT) and 2-amino-1,3,5-triazine (2-AT) in solution have been systematically explored using the QM(MS-CASPT2//CASSCF)/MM approach. The excited-state nonradiative relaxation mechanisms for the initially photoexcited S1(ππ*) state decay back to the S0 state are proposed in terms of the present computed minima, surface crossings (conical intersections and singlet-triplet crossings), and excited-state decay paths in the S1, S2, T1, T2, and S0 states. Upon photoexcitation to the bright S1(ππ*) state, 2,4-DT quickly relaxes to its S1 minimum and then overcomes a small energy barrier of 5.1 kcal mol-1 to approach a S1/S0 conical intersection, where the S1 system hops to the S0 state through S1 → S0 internal conversion (IC). In addition, at the S1 minimum, the system could partially undergo intersystem crossing (ISC) to the T1 state, followed by further ISC to the S0 state via the T1/S0 crossing point. In the T1 state, an energy barrier of 7.9 kcal mol-1 will trap 2,4-DT for a while. In parallel, for 2-AT, the system first relaxes to the S1 minimum and then S1 → S0 IC or S1 → T1 → S0 ISCs take place to the S0 state by surmounting a large barrier of 15.3 kcal mol-1 or 11.9 kcal mol-1, respectively, which heavily suppress electronic transition to the S0 state. Different from 2,4-DT, upon photoexcitation in the Franck-Condon region, 2-AT can quickly evolve in an essentially barrierless manner to nearby S2/S1 conical intersection, where the S2 and T1 states can be populated. Once it hops to the S2 state, the system will overcome a relatively small barrier (6.6 kcal mol-1vs. 15.3 kcal mol-1) through IC to the S0 state. Similarly, an energy barrier of 11.9 kcal mol-1 heavily suppresses the T1 state transformation to the S0 state. The present work manifests that the amination/deamination of the triazine rings can affect some degree of different vertical and adiabatic excitation energies and nonradiative decay pathways in solution. It not only rationalizes excited-state decay dynamics of 2,4-DT and 2-AT in aqueous solution but could also provide insights into the understanding of the photophysics of aza-nucleobases.

Modulating Cation Migration and Deposition with Xylitol Additive and Oriented Reconstruction of Hydrogen Bonds for Stable Zinc Anodes.

Highly reversible plating/stripping in aqueous electrolytes is one of the critical processes determining the performance of Zn-ion batteries, but it is severely impeded by the parasitic side reaction and dendrite growth. Herein, a novel electrolyte engineering strategy is first proposed based on the usage of 100 mM xylitol additive, which inhibits hydrogen evolution reaction and accelerates cations migration by expelling active H2 O molecules and weakening electrostatic interaction through oriented reconstruction of hydrogen bonds. Concomitantly, xylitol molecules are preferentially adsorbed by Zn surface, which provides a shielding buffer layer to retard the sedimentation and suppress the planar diffusion of Zn2+ ions. Zn2+ transference number and cycling lifespan of Zn∥Zn cells have been significantly elevated, overwhelmingly larger than bare ZnSO4 . The cell coupled with a NaV3 O8 cathode still behaves much better than the additive-free device in terms of capacity retention.

Photoinduced boron atom insertion of benzocyclobutene forming an unprecedented fused boron heterocyclic radical.

Two novel boron heterocyclic radicals, an addition bicyclo[4.2.1]octa-1,3,5-trien-1-yl-borane radical (A) and an insertion 7-1H-borolo[1,2-a]borinine radical (B), were synthesized, and characterized in the reaction of atomic boron with benzocyclobutene. Species B involving a fused boron heterocyclic was spectroscopically characterized for the first time. This work is a new approach for boron-mediated molecular editing and the synthesis of fused boron heterocyclic compounds.

Combined QM (MS-CASPT2)/MM studies on photocyclization and photoisomerization of a fulgide derivative in toluene solution.

Photocyclization and photoisomerization of fulgides have been extensively studied experimentally and computationally due to their significant potential applications for example as photoswitches in memory devices. However, the reported excited-state decay mechanisms of fulgides do not include the effects of solvation explicitly to date. Herein, calculations using the high-level MS-CASPT2//CASSCF method were conducted to explore the photoinduced excited-state decay processes of the Eα conformer of a fulgide derivative in toluene with solvent effects treated by implicit PCM and explicit QM/MM models, respectively. Several minima and conical intersections were optimized successfully in and between the S0 and S1 states; then, two nonadiabatic excited-state decay channels that could efficiently drive the system to the ground state were proposed based on the excited-state ring-closure and isomerization paths. In addition, we also found that in the ring-closure path, the potential energy surface is essentially barrierless before approaching the conical intersection, while it needs to overcome a small energy barrier along the E → Z photoisomerization path for the nonadiabatic S1 → S0 internal conversion process. The present computational results could provide useful mechanistic insights into the photoinduced cyclization and isomerization reactions of fulgide and its derivatives.

[Spatialization and Spatio-temporal Dynamics of Energy Consumption Carbon Emissions in China].

The adverse effects of global climate change on human production and life are becoming increasingly prominent. Responding to climate change has become a severe challenge faced by human society, and the reduction in greenhouse gas emissions has gradually become a common action by all countries. Therefore, analyzing carbon emissions through scientific methods has become an important foundation for responding to the national "dual carbon" strategy. This study used provincial-level carbon emission statistics, combined with nighttime light data and population data, and assigned carbon emissions to the grid scale. It also analyzed the temporal and spatial characteristics and evolution characteristics of carbon emissions in China in 2000, 2005, 2010, 2015, and 2018, as well as the correlation between carbon emissions and the economy. The results showed that:① from 2000 to 2018, the total CO2 emissions in China continued to grow, but the growth rate slowed over time. The average annual growth rate of carbon emissions dropped from 9.9% in 2000-2010 to 7.4% in 2010-2018. From the perspective of spatial distribution, carbon-free areas were mainly distributed in the northwest uninhabited area and northeast forest and mountainous areas, low-carbon emissions were mainly distributed in the vast small and medium-sized cities and towns, and high-carbon emissions were concentrated in northern, central, eastern coastal, and western provincial capitals and urban agglomerations. ② Carbon emissions had high-value or low-value agglomerations at prefecture-level cities; this agglomeration tended to stabilize as a whole and had strengthened after 2005. Low-low agglomeration areas were mainly distributed in the western contiguous areas and Hainan Island. With economic and social development, low-low agglomeration areas began to fragment and reduce in size; high-high agglomeration areas were mainly distributed in the Beijing-Tianjin-Hebei urban agglomeration, Taiyuan urban agglomeration, Yangtze River Delta urban agglomerations, and Pearl River Delta urban agglomerations, and the scale was gradually strengthened and consolidated; high-low and low-high agglomeration areas mainly appeared in neighboring cities with large differences in economic development levels. ③ Carbon emissions in most parts of China were relatively stable. The areas where carbon emissions had changed were mainly distributed in the peripheral areas of provincial capitals and key cities, and there was a circle structure with no changes in the central urban area and changes in carbon emissions in the peripheral areas. ④ The overall process of urban development in China from 2000 to 2018 followed a shift from "low emission-low income" to "high emission-low income" to "high emission-high income" and finally to "low emission-high income." The growth rate of carbon emissions in China is slowing down. Under the background of the "dual carbon" strategy, different regions face different carbon emission reduction tasks and pressures due to different carbon emission situations. Therefore, the differentiated carbon emissions policy should be implemented by regions and industries.

Quantum mechanics/molecular mechanics studies on excited state decay pathways of 5-azacytosine in aqueous solution.

In this work, we have used the QM(CASPT2//CASSCF)/MM approach to study the photophysical properties and relaxation mechanism of 5-azacytosine (5-AC) in aqueous solution. Based on the relevant minimum-energy structures and intersection structures, and excited-state decay paths in the S1, S2, T1, T2, and S0 states, several feasible excited-state nonradiative decay channels from the initially populated S2(ππ*) state are proposed. Two major channels are singlet-mediated nonradiative pathways, in which the S2 system will internally convert (IC) to the S0 state directly or mediated by the 1nπ* state via a 1ππ*/1nπ* conical intersection. The minor ones are related to intersystem crossing (ISC) processes. The system would populate to the T1 state via the S2 → S1 → T1 or S2 → T2 → T1 ISC process, followed by further decay to the S0 state via the transition from T1 to S0. However, due to small spin-orbit couplings (SOCs) at the singlet-triplet crossing points, the related ISC would be less efficient and probably take longer. The present work rationalizes the ultrafast excited-state decay dynamics of 5-AC in aqueous solution and its low quantum yields of triplets and fluorescence. It provides important mechanistic insights into understanding 5-AC's derivatives and analogues.

Selectively Identifying Exposed-over-Unexposed C-C+ Pairs in Human Telomeric i-Motif Structures with Length-Dependent Polymorphism.

The i-motif structure (iM) has attracted much attention, because of its in vivo bioactivity and wide in vitro applications such as DNA-based switches. Herein, the length-dependent folding of cytosine-rich repeats of the human telomeric 5'-(CCCTAA)n-1CCC-3' (iM-n, where n = 2-8) was fully explored. We found that iM-4, iM-5, and iM-8 mainly form the intramolecular monomer iM structures, while a tetramolecular structure populates only for iM-3. However, iM-6 and iM-7 have the potential to fold as well into the dimeric iM structures besides the monomer ones. The natural hypericin (Hyp) was used as the polymorphism-selective probe to recognize the iM structures. Interestingly, only iM-3, iM-6, and iM-7 can efficiently switch on the Hyp fluorescence by specifically binding with the outmost C-C+ base pairs that are exposed directly to solution. However, other iM structures that fold in a way with a coverage of the outmost C-C+ pairs by loop sequences are totally unavailable for the Hyp binding. Theoretical modeling indicates that adaptive π-π and cation-π interactions contribute to the Hyp recognition toward the exposed C-C+ pairs. This specific iM recognition can be boosted by a photocatalytic DNAzyme construct. Our work provides a reliable fluorescence method to selectively explore the polymorphism of iM structures.

A theoretical insight into excited-state decay and proton transfer of p-nitrophenylphenol in the gas phase and methanol solution.

The excited-state decay (ESD) and proton transfer (EPT) of p-nitrophenylphenol (NO2-Bp-OH), especially in the triplet states, were not characterized with high-level theoretical methods to date. Herein, the MS-CASPT2//CASSCF and QM(MS-CASPT2//CASSCF)/MM methods were employed to gain an atomic-level understanding of the ESD and EPT of NO2-Bp-OH in the gas phase and its hydrogen-bonded complex in methanol. Our calculation results revealed that the S1 and S2 states of NO2-Bp-OH are of 1ππ* and 1nπ* characters at the Franck-Condon (FC) point, which correspond to the ICT-EPT and intramolecular charge-transfer (ICT) states in spectroscopic experiments. The former state has a charge-transfer property that could facilitate the EPT reaction, while the latter one might be unfavorable for EPT. The vertical excitation energies of these states are almost degenerate at the FC region and the electronic configurations of 1ππ* and 1nπ* will exchange from the S1 FC region to the S1 minimum, which means that the 1nπ* state will participate in ESD once NO2-Bp-OH departs from the S1 FC region. Besides, we found that three triplets lie below the first bright state and will play very important roles in intersystem crossing processes. In terms of several pivotal surface crossings and relevant linearly interpolated internal coordinate (LIIC) paths, three feasible but competing ESD channels that could effectively lead the system to the ground state or the lowest triplet state were put forward. Once arrived at the T1 state, the system has enough time and internal energy to undergo the EPT reaction. The methanol solvent has a certain effect on the relative energies and spin-orbit couplings, but does not qualitatively change the ESD processes of NO2-Bp-OH. By contrast, the solvent effects will remarkably stabilize the proton-transferred product by the hydrogen bond networks and assist to form the triplet anion. Our present work would pave the road to properly understand the mechanistic photochemistry of similar hydroxyaromatic compounds.

Generalized Ab Initio Nonadiabatic Dynamics Simulation Methods from Molecular to Extended Systems.

Nonadiabatic dynamics simulation has become a powerful tool to describe nonadiabatic effects involved in photophysical processes and photochemical reactions. In the past decade, our group has developed generalized trajectory-based ab initio surface-hopping (GTSH) dynamics simulation methods, which can be used to describe a series of nonadiabatic processes, such as internal conversion, intersystem crossing, excitation energy transfer and charge transfer of molecular systems, and photoinduced nonadiabatic carrier dynamics of extended systems with and without spin-orbit couplings. In this contribution, we will first give a brief introduction to our recently developed methods and related numerical implementations at different computational levels. Later, we will present some of our latest applications in realistic systems, which cover organic molecules, biological proteins, organometallic compounds, periodic organic and inorganic materials, etc. Final discussion is given to challenges and outlooks of ab initio nonadiabatic dynamics simulations.

Mechanistic photophysics and photochemistry of unnatural bases and sunscreen molecules: insights from electronic structure calculations.

Photophysics and photochemistry are basic subjects in the study of light-matter interactions and are ubiquitous in diverse fields such as biology, energy, materials, and environment. A full understanding of mechanistic photophysics and photochemistry underpins many recent advances and applications. This contribution first provides a short discussion on the theoretical calculation methods we have used in relevant studies, then we introduce our latest progress on the mechanistic photophysics and photochemistry of two classes of molecular systems, namely unnatural bases and sunscreens. For unnatural bases, we disclose the intrinsic driving forces for the ultrafast population to reactive triplet states, impacts of the position and degree of chalcogen substitutions, and the effects of complex environments. For sunscreen molecules, we reveal the photoprotection mechanisms that dissipate excess photon energy to the surroundings by ultrafast internal conversion to the ground state. Finally, relevant theoretical challenges and outlooks are discussed.