Intramolecular boron-locking strategy induced remarkable first hyperpolarizability: role of torsion angles between donor and acceptor units.

In conventional strategies to design donor-acceptor (D-A) organic molecules with a large electronic contribution to the first hyperpolarizability (ß), the effects of the torsion angles (θ1 and θ2) between donor and acceptor moieties are barely considered. To address this issue, in this work, a promising and novel intramolecular boron-locking strategy combined with the different locking groups of different acceptors to control θ1 and θ2, has been proposed to make D-A organic molecules with large ß values. Intriguingly, reducing the torsion angles will make the ß value of the pyridiny thiophene triphenylamine unit (Py-Th-TPA) dramatically increase up to 94%, which is mainly ascribed to the smaller θ1 and θ2 leading to lower excited energy of the crucial excited state, and enhanced charge transfer (CT) from TPA to Py-Th moieties, and finally greatly increase the donor and acceptor part contributions to ß. Correlation between the difference, |θ1 - θ2| and ß, provides a large coefficient of determination, R2 = 0.78, which demonstrates that |θ1 - θ2| can be regarded as a potential descriptor for designing nonlinear optics (NLO) materials with D-A architecture. Clearly, we uncovered that θ1 and θ2 play a crucial role in the performance of NLO materials with D-A fragments.

On-the-fly simulation of time-resolved fluorescence spectra and anisotropy.

We combine on-the-fly trajectory surface hopping simulations and the doorway-window representation of nonlinear optical response functions to create an efficient protocol for the evaluation of time- and frequency-resolved fluorescence (TFRF) spectra and anisotropies of the realistic polyatomic systems. This approach gives the effective description of the proper (e.g., experimental) pulse envelopes, laser field polarizations, and the proper orientational averaging of TFRF signals directly from the well-established on-the-fly nonadiabatic dynamic simulations without extra computational cost. To discuss the implementation details of the developed protocol, we chose cis-azobenzene as a prototype to simulate the time evolution of the TFRF spectra governed by its nonadiabatic dynamics. The results show that the TFRF is determined by the interplay of several key factors, i.e., decays of excited-state populations, evolution of the transition dipole moments along with the dynamic propagation, and scaling factor of the TFRF signals associated with the cube of emission frequency. This work not only provides an efficient and effective approach to simulate the TFRF and anisotropies of realistic polyatomic systems but also discusses the important relationship between the TFRF signals and the underlining nonadiabatic dynamics.

Torsion Angles between Donor and Acceptor Moieties as a Descriptor for Designing Nonlinear Optics and Thermally Activated Delayed Fluorescence Materials.

The performances of nonlinear optics (NLO) and thermally activated delayed fluorescence (TADF) materials are strongly related to the torsion angles (θ) between donor (D) and acceptor (A) moieties in D-A architecture molecules. However, the underlying relationships connecting θ to the performances of NLO/TADF materials remain unclear. Herein, we present a comprehensive theoretical study on NLO/TADF materials composed of a series of D-A backbone molecules (TPAAP/TPAAQ series and AQ-DMAC/AQ-MeFAC series) to shed light on these relationships. It is found that changing θ via the intramolecular locking strategy can greatly influence values of the first hyperpolarizability (ß) and singlet-triplet energy gap (ΔEST), further leading to better/worse performances of NLO/TADF materials, respectively. Intriguingly, a more detailed analysis indicates that the variation trends between θ and ß/ΔEST are changeable in low θ regions, exhibiting volcano-like relationships. The large coefficients of determination (R2, ranging from 0.76 to 0.93) suggest that this experimentally measurable parameter (θ) can be used as a promising descriptor to evaluate the performances of related materials. Following the revealed θ-ß/θ-ΔEST correlations, the optimal/worst torsion angles for different materials are identified. These findings highlight the importance of the intrinsic structure-performance relationships, thus providing novel design strategies for high-performance NLO/TADF materials.

Kinetics and mechanism of triclocarban degradation by the chlorination process: Theoretical calculation and experimental verification.

Triclocarban (TCC) is an antimicrobial agent commonly used in many household and personal care products, and has been found persistent in the aquatic environment. Here we elucidate the kinetics and mechanism of TCC degradation during chlorination process by density functional theory (DFT) calculation and experimental verification. Results showed that hypochlorous acid (HOCl)/hypochlorite (OCl-) reacted with TCC via Cl-substitution, OH-substitution and C-N bond cleavage pathways. The reactivity of OCl- (2.80 × 10-7 M-1 s-1) with TCC was extremely low and HOCl (1.96 M-1 s-1) played the dominant role in TCC chlorination process. The N site of TCC was the most reactive site for chlorination. The second-order rate constants, which are determined using density functional theory (DFT) (kTCC-chlorineC, 1.96 M-1 s-1), can be separated into reaction rate constants related to the reactions of HOCl and OCl- with different isomers of TCC (TCC2 and TCC6). The obtained kTCC-chlorineC was consistent with the experimental determined second-order rate constant (kTCC-chlorineE, 3.70 M-1 s-1) in chlorination process. Eight transformation products (TP348, TP382, TP127, TP161, TP195, TP330, TP204, and TP296) were experimentally detected for chlorination of TCC, which could also be predicted by DFT calculation. Explicit water molecules participated in the chlorination reaction by transmitting the proton and connecting with TCC, HOCl/OCl- and other H2O molecules, and obviously reduced the energy barrier of chlorination.

Trajectory Propagation of Symmetrical Quasi-classical Dynamics with Meyer-Miller Mapping Hamiltonian Using Machine Learning.

The long short-term memory recurrent neural network (LSTM-RNN) approach is applied to realize the trajectory-based nonadiabatic dynamics within the framework of the symmetrical quasi-classical dynamics method based on the Meyer-Miller mapping Hamiltonian (MM-SQC). After construction, the LSTM-RNN model allows us to propagate the entire trajectory evolutions of all involved degrees of freedoms (DOFs) from initial conditions. The proposed idea is proven to be reliable and accurate in the simulations of the dynamics of several site-exciton electron-phonon coupling models and three Tully's scattering models. It indicates that the LSTM-RNN model perfectly captures the dynamical information on the trajectory evolution in the MM-SQC dynamics. Our work proposes a novel machine learning approach in the simulation of trajectory-based nonadiabatic dynamic of complex systems with a large number of DOFs.

The impact of different geometrical restrictions on the nonadiabatic photoisomerization of biliverdin chromophores.

The photoisomerization mechanism of the chromophore of bacterial biliverdin (BV) phytochromes is explored via nonadiabatic dynamics simulation by using the on-the-fly trajectory surface-hopping method at the semi-empirical OM2/MRCI level. Particularly, the current study focuses on the influence of geometrical constrains on the nonadiabatic photoisomerization dynamics of the BV chromophore. Here a rather simplified approach is employed in the nonadiabatic dynamics to capture the features of geometrical constrains, which adds mechanical restrictions to the specific moieties of the BV chromophore. This simplified method provides a rather quick approach to examine the influence of geometrical restrictions on photoisomerization. As expected, different constrains bring distinctive influences on the photoisomerization mechanism of the BV chromophore, giving either strong or minor modification of both involved reaction channels and excited-state lifetimes after the constrains are added in different ring moieties. These observations not only contribute to the primary understanding of the role of the spatial restriction caused by biological environments in photoinduced dynamics of the BV chromophore, but also provide useful ideas for the artificial regulation of the photoisomerization reaction channels of phytochrome proteins.

Automatic Evolution of Machine-Learning-Based Quantum Dynamics with Uncertainty Analysis.

The machine learning approaches are applied in the dynamical simulation of open quantum systems. The long short-term memory recurrent neural network (LSTM-RNN) models are used to simulate the long-time quantum dynamics, which are built based on the key information of the short-time evolution. We employ various hyperparameter optimization methods, including simulated annealing, Bayesian optimization with tree-structured parzen estimator, and random search, to achieve the automatic construction and adjustment of the LSTM-RNN models. The implementation details of three hyperparameter optimization methods are examined, and among them, the simulated annealing approach is strongly recommended due to its excellent performance. The uncertainties of the machine learning models are comprehensively analyzed by the combination of bootstrap sampling and Monte Carlo dropout approaches, which give the prediction confidence of the LSTM-RNN models in the simulation of the open quantum dynamics. This work builds an effective machine learning approach to simulate the dynamics evolution of open quantum systems. In addition, the current study provides an efficient protocol to build optimal neural networks and estimate the trustiness of the machine learning models.

An intramolecular-locked strategy for designing nonlinear optical materials with remarkable first hyperpolarizability.

To meet the expanding demands of high performance nonlinear optical (NLO) materials, an unprecedented intramolecular-locked strategy is proposed to design NLO materials with remarkable static first hyperpolarizability (ß0). This strategy means that importing a large steric hindrance group diphenylmethane (DPM) decreases the torsion angles (θ) between the donor {triphenylamine (TPA)} and acceptor {9-H-thioxanthen-9-one-10,10-dioxide (TXO)} units, as well as between the donor (TPA) and π-bridge (benzene) fragments. The decrease of θ can accelerate the intramolecular charge transfer and enhance the contributions of the TPA, TXO and quinoxaline-6,7-dicarbo-nitrile (QCN) fragments to the axial component of the ß0 value, and then the ß0 values of TPA-TXO (ß0 = 10 762 au) and TPA-QCN (ß0 = 22 495 au) are increased by 14.9% and 34.4%, respectively. Overall, the intramolecular-locked strategy is very effective for designing high performance NLO materials.

Regularized Localized Molecular Orbitals in a Divide-and-Conquer Approach for Linear Scaling Calculations.

Non-orthogonal localized molecular orbitals (NOLMOs) have been employed as building blocks for the divide-and-conquer (DC) linear scaling method. The NOLMOs are calculated from subsystems and used for constructing the density matrix (DM) of the entire system, instead of the subsystem DM in the original DC approach. Also, unlike the original DC method, the inverse electronic temperature parameter ß is not needed anymore. Furthermore, a new regularized localization approach for NOLMOs has been developed, in which the localization cost function is a sum of the spatial spread function, as in the Boys method, and the kinetic energy, as a regularization measure to limit the oscillation of the NOLMOs. The optimal weight of the kinetic energy can be determined by optimization with analytical gradients. The resulting regularized NOLMOs have enhanced smoothness and better transferability because of reduced kinetic energies. Compared with the original DC, while NOLMO-DC has a similar computational linear scaling cost, the accuracy of NOLMO-DC is better by several orders of magnitude for large conjugated systems and by about 1 order of magnitude for other systems. The NOLMO-DC method is thus a promising development of the DC approach for linear scaling calculations.

Second-Order Nonlinear Optics Response of the Boron-Dipyrromethenes-Based Mislinked Expanded Porphyrins: Revealing the Role of the -BF2 Group.

Here, the mislinked expanded porphyrins singly (labeled A) and doubly (labeled B) neo-confused [22]smaragdyrin, the boron-dipyrromethenes-based mislinked expanded porphyrins singly (labeled C) and doubly (labeled D) neo-confused [22]smaragdyrin, where both C and D include a -BF2 group, are chosen to serve as the study objects, and theoretical calculations are carried out to study the role of the -BF2 group in the second-order nonlinear optics (NLO) behaviors. Results highlighted that the -BF2 group plays an important role for the second-order behaviors in mislinked expanded porphyrins; namely, embedding the -BF2 group well enhanced the hyper-Rayleigh scattering (HRS) value {ßHRS(0;0,0)}, C{ßHRS(0;0,0)}A{ßHRS(0;0,0)} = 2.0 and D{ßHRS(0;0,0)}B{ßHRS(0;0,0)} = 2.9, main owning to the fact that installing -BF2 increases the electron delocalization degree and decreases the excited energy of the crucial excited state.

Ultrafast Internal Conversion Dynamics through the on-the-Fly Simulation of Transient Absorption Pump-Probe Spectra with Different Electronic Structure Methods.

An on-the-fly surface-hopping simulation protocol is developed for the evaluation of transient absorption (TA) pump-probe (PP) signals of molecular systems exhibiting internal conversion to the electronic ground state. We study the nonadiabatic dynamics of azomethane and the associating TA PP spectra at three levels of the electronic-structure theory, OM2/MRCI, SA-CASSCF, and XMS-CASPT2. The impact of these methods on the population dynamics and time-resolved TA PP signals is substantially different. This difference is attributed to the strong non-Condon effects that must be taken into account for the proper understanding and interpretation of time-resolved TA PP signals of nonadiabatic polyatomic systems. This shows that the combination of the dynamical and spectral simulations definitely provides more accurate and detailed information on the microscopic mechanisms of photophysical and photochemical processes. Hence the simulation of time-resolved spectroscopic signals provides another important dimension to examine the accuracy of quantum chemistry methods.

Theoretical studies on triplet-state driven dissociation of formaldehyde by quasi-classical molecular dynamics simulation on machine-learning potential energy surface.

The H-atom dissociation of formaldehyde on the lowest triplet state (T1) is studied by quasi-classical molecular dynamic simulations on the high-dimensional machine-learning potential energy surface (PES) model. An atomic-energy based deep-learning neural network (NN) is used to represent the PES function, and the weighted atom-centered symmetry functions are employed as inputs of the NN model to satisfy the translational, rotational, and permutational symmetries, and to capture the geometry features of each atom and its individual chemical environment. Several standard technical tricks are used in the construction of NN-PES, which includes the application of clustering algorithm in the formation of the training dataset, the examination of the reliability of the NN-PES model by different fitted NN models, and the detection of the out-of-confidence region by the confidence interval of the training dataset. The accuracy of the full-dimensional NN-PES model is examined by two benchmark calculations with respect to ab initio data. Both the NN and electronic-structure calculations give a similar H-atom dissociation reaction pathway on the T1 state in the intrinsic reaction coordinate analysis. The small-scaled trial dynamics simulations based on NN-PES and ab initio PES give highly consistent results. After confirming the accuracy of the NN-PES, a large number of trajectories are calculated in the quasi-classical dynamics, which allows us to get a better understanding of the T1-driven H-atom dissociation dynamics efficiently. Particularly, the dynamics simulations from different initial conditions can be easily simulated with a rather low computational cost. The influence of the mode-specific vibrational excitations on the H-atom dissociation dynamics driven by the T1 state is explored. The results show that the vibrational excitations on symmetric C-H stretching, asymmetric C-H stretching, and C=O stretching motions always enhance the H-atom dissociation probability obviously.

Simulation of Open Quantum Dynamics with Bootstrap-Based Long Short-Term Memory Recurrent Neural Network.

The recurrent neural network with the long short-term memory cell (LSTM-NN) is employed to simulate the long-time dynamics of open quantum systems. The bootstrap method is applied in the LSTM-NN construction and prediction, which provides a Monte Carlo estimation of a forecasting confidence interval. Within this approach, a large number of LSTM-NNs are constructed by resampling time-series sequences that were obtained from the early stage quantum evolution given by numerically exact multilayer multiconfigurational time-dependent Hartree method. The built LSTM-NN ensemble is used for the reliable propagation of the long-time quantum dynamics, and the simulated result is highly consistent with the exact evolution. The forecasting uncertainty that partially reflects the reliability of the LSTM-NN prediction is also given. This demonstrates the bootstrap-based LSTM-NN approach is a practical and powerful tool to propagate the long-time quantum dynamics of open systems with high accuracy and low computational cost.

Prediction of the excited-state reaction channels in photo-induced processes of nitrofurantoin using first-principle calculations and dynamics simulations.

The understanding of the photochemistry of antibiotic compounds is important because it gives the direct information on the possible environmental pollution caused by them. Due to their large size, the theoretical studies of their excited-state reactions are rather challenging. In current work, we combined the on-the-fly trajectory surface-hopping dynamics, conical-intersection optimizations and excited-state pathway calculations to study the photochemistry of the trans-isomer of nitrofurantoin, a widely-used drug to treat the urinary tract infections. The dynamics-then-pathway approach was taken. First the trajectory surface hopping dynamics at the state-averaged complete-active-space self-consistent-field (SA-CASSCF) level with small active space and small basis sets were run. Second, the minimum-energy conical-intersection optimizations were performed. Finally the excited pathways from the Frank-Condon region to different reaction channels were built at the multi-state multi-reference second-order perturbation (MS-CASPT2) level with large active space and large basis set. Several possible channels responsible for the photo-induced reaction mechanism of the trans-nitrofurantoin were obtained, including the cleavage of the NO bond of the NO2 moiety, the photoisomerization at the central CN bond, and other internal conversion channels. Our findings give some preliminary explanations on available experimental observations. It is also demonstrates that the current theoretical approach is a powerful tool to explore the excited-state reactions in the photochemistry of media-sized or large-sized drug compounds.

Hazardous impact of diclofenac on mammalian system: Mitigation strategy through green remediation approach.

Diclofenac is an anti-inflammatory drug used as an analgesic. It is often detected in various environmental sources around the world and is considered as one of the emerging contaminants (ECs). This paper reviews the distribution of diclofenac at high concentrations in diverse environments and its adverse ecological impact. Recent studies observed strong evidence of the hazardous effect of diclofenac on mammals, including humans. Diclofenac could cause gastrointestinal complications, neurotoxicity, cardiotoxicity, hepatotoxicity, nephrotoxicity, hematotoxicity, genotoxicity, teratogenicity, bone fractures, and skin allergy in mammals even at a low concentration. Collectively, this comprehensive review relates the mode of toxicity, level of exposure, and route of administration as a unique approach for addressing the destructive consequence of diclofenac in mammalian systems. Finally, the mitigation strategy to eradicate the diclofenac toxicity through green remediation is critically discussed. This review will undoubtedly shed light on the toxic effects of pseudo-persistent diclofenac on mammals as well as frame stringent guidelines against its common usage.

Polymer-capped CdSe/ZnS quantum dots for the sensitive detection of Cu2+ and Hg2+ and the quenching mechanism.

In this work, poly(styrene-co-maleic anhydride)-capped CdSe/ZnS quantum dots (QDs) aminolyzed with ethanolamine are proposed as fluorescent probes for the detection of Cu2+ and Hg2+, and two different quenching mechanisms are discussed in detail. The coordination abilities of the surface polymer of CdSe/ZnS QDs and two metal ions are calculated by density functional theory (DFT). The photoinduced electron transfer from excited QDs to Cu2+ unoccupied orbitals is enhanced due to the coordination between Cu2+ and the surface polymer of QDs. The electron transfer consumes non-radiative energy and performs fluorescence quenching. For Hg2+, the formation of HgS and the slight aggregation of polymer-coated CdSe/ZnS QDs lead to fluorescence quenching. The probe is sensitive to both Cu2+ and Hg2+, and the response can be detected within 1 min without adjusting the pH. With the addition of a masking agent, Cu2+ and Hg2+ can be exclusively detected in coexistence with another ion. For Cu2+, a linear relation in the concentration ranging from 0.02 to 0.7 µM was found between the relative fluorescence intensity (F0/F) and the concentration of Cu2+; the limit of detection (S/N = 3) is 6.94 nM. For Hg2+, a linear relation ranging from 0.1 to 1.4 µM was found between ln(F0/F) and the concentration of Hg2+; the limit of detection is 20.58 nM.

Tuning the first hyperpolarizability of hexaphyrins with different connections of mislinked pyrrole units: a theoretical study.

In the satisfactory design and synthesis of high-performance nonlinear optical (NLO) materials, for meeting the rapidly expanding demands of optoelectronic devices, a deeper understanding of the relationship between the structures and NLO properties has become a key issue. Herein, five novel mislinked expanded hexaphyrins with different connections of pyrrole units are selected to study the relationship between the structures and NLO properties. These five mislinked expanded hexaphyrins are neo-fused, neo-confused hexaphyrins, singly, doubly, and triply N-confused hexaphyrins. From theoretical results, the order of the static first hyperpolarizability (ß0) values is found to be: neo-fused hexaphyrin (ß0 = 4163 a.u.) < neo-confused hexaphyrin (ß0 = 5494 a.u.) < singly N-confused hexaphyrin (ß0 = 6510 a.u.) < doubly N-confused hexaphyrin (ß0 = 15 130 a.u.) < triply N-confused hexaphyrin (ß0 = 26 095 a.u.). Furthermore, ß0 values of the doubly and triply N-confused hexaphyrins are improved 2.1 and 3.7 times over that of their usual parent hexaphyrin (ß0 = 7120 a.u.), respectively. It is worth noting that increasing mislinked connection numbers and changing mislinked connection ways of the pyrrole units in these mislinked expanded hexaphyrins plays a crucial role in the tune of their second-order NLO responses.

Nonadiabatic molecular dynamics simulation for the ultrafast photoisomerization of dMe-OMe-NAIP based on TDDFT on-the-fly potential energy surfaces.

Global switching on-the-fly trajectory surface hopping molecular dynamics simulation was performed on the accurate TD-B3LYP/6-31G* potential energy surfaces for E-to-Z and Z-to-E photoisomerization of dMe-OMe-NAIP up to S1(ππ*) excitation. The present TD(DFT) simulation provides accurate calculation for conical intersections between the first-excited and ground states. Thus, simulated quantum yield and lifetime of 0.23 and 620 fs (0.15 and 600 fs) for E-to-Z (Z-to-E) isomerization are in good (relatively good) agreement with experimental observation of 0.25 and 480 fs (0.24 and 430 fs), respectively. Simulated results reveal that photoisomerization pathways are initially uphill to conical intersection zones on the S1 potential energy surface and then downhill to product zones. Three types of representative conical intersections are found for determining photoisomerization mechanisms: one is the rotation type responsible for reactive isomerization and the other two are close to E and Z configurations, respectively, only for nonreactive isomerization. The present conclusions can be held in general for similar large NAIP systems of photoinduced isomerization based on E and Z configurations.

Photoisomerization-mechanism-associated excited-state hydrogen transfer in 2'-hydroxychalcone revealed by on-the-fly trajectory surface-hopping molecular dynamics simulation.

By performing global-switching on-the-fly trajectory surface-hopping molecular dynamics simulation at the OM2/MRCI (14,15) quantum level, we probed the S3(ππ*) photoisomerization mechanisms associated with excited-state intramolecular hydrogen transfer for 2'-hydroxychalcone (2HC) within the interwoven conical intersection networks from four singlet electronic states (S3, S2, S1, and S0). The simulated quantum yields of 0.03 for cis-to-trans and zero for trans-to-cis photoisomerization were due to almost all the conical intersections being localized either in the cis-2HC or in trans-2HC region, and there was little chance for sampling trajectories to reach the rotation conical intersection (S1/S0) in between cis-2HC and trans-2HC that is key for reactive isomerization. The potential energy well on the S1 state in the trans-2HC region prevents trajectories from trans-to-cis photoisomerization, while the fact there is no well on S1 state in cis-2HC region opens a few chances for trajectories to reach the rotation conical intersections. The present simulation found that excited-state intramolecular hydrogen transfers in 2HC have a negative impact for reactive isomerization, and that hydrogen transfers take place on the S1 state, while back-transfer on the S0 state prevents sampling trajectories reaching rotational conical intersections. It was realized that it could be possible to enhance the quantum yield of 2HC photoisomerization by suppressing the hydrogen transfer (such as by changing an electron-donating substitution or adjusting the substitution position to decrease the acidity of the phenol group). From a perspective view of the potential energy surfaces, the theoretical design of such 2HC derivatives could enhance (control) the quantum yield by shifting the conical intersections away from the cis- and trans-region.

Boosting the Quantum Efficiency of Ultralong Organic Phosphorescence up to 52 % via Intramolecular Halogen Bonding.

Ultralong organic phosphorescence (UOP) has attracted increasing attention due to its potential applications in optoelectronics, bioelectronics, and security protection. However, achieving UOP with high quantum efficiency (QE) over 20 % is still full of challenges due to intersystem crossing (ISC) and fast non-radiative transitions in organic molecules. Here, we present a novel strategy to enhance the QE of UOP materials by modulating intramolecular halogen bonding via structural isomerism. The QE of CzS2Br reaches up to 52.10 %, which is the highest afterglow efficiency reported so far. The crucial reason for the extraordinary QE is intramolecular halogen bonding, which can not only effectively enhance ISC by promoting spin-orbit coupling, but also greatly confine motions of excited molecules to restrict non-radiative pathways. This work provides a reasonable strategy to develop highly efficient UOP materials for practical applications.