Elucidating the Photophysics and Nonlinear Optical Properties of a Novel Azo Prototype for Possible Photonic Applications: A Quantum Chemical Analysis.

The photophysics and nonlinear optical responses of a novel nitrothiazol-methoxyphenol molecule were investigated using density functional theory (DFT) and time-dependent density functional theory (TD-DFT) methods with the polarizable continuum model to take the solvent effect into account. Special attention is paid to the description of the lowest absorption band, characterized as a strong π → π* state in the visible region of the spectrum. The TD-DFT emission spectrum analysis reveals a significant Stokes shift of more than 120 nm for the π → π* state in gas phase condition.

Influence of Solvents and Halogenation on ESIPT of Benzimidazole Derivatives for Designing Turn-on Fluorescence Probes.

This work reports a theoretical investigation of the solvent polarity as well as the halogenation of benzimidazole derivatives during excited state intramolecular proton transfer (ESIPT). It details how the environment and halogen substitution may contribute to the efficiency of ESIPT upon keto-enol tautomerism and exploits this effect to design fluorescence sensing. For this purpose, we first examine the conformational equilibrium of benzimidazole derivatives containing different halogen atoms, which results in intramolecular proton transfer, using density-functional theory (DFT) combined with the polarizable continuum model (PCM).

Calculation of the one- and two-photon absorption spectra of water-soluble stilbene derivatives using a multiscale QM/MM approach.

We calculated the one- (OPA) and two-photon absorption (TPA) spectra of two large water-soluble stilbene derivatives presenting TPA cross sections of about 400 GM. However, the reported experimental TPA spectra present a spectral gap region, and a theoretical study of these promising molecules seems now timely and relevant. These molecules are composed of 200 or more atoms, becoming a challenge to obtain the TPA spectra even using density functional theory at the time-dependent quadratic response formalism.

The Solute Polarization and Structural Effects on the Nonlinear Optical Response of Based Chromone Molecules.

The solute polarization due to solvent is a an electrostatic quantum effect that impacts diverse molecular properties, including the nonlinear optical response of a material. An iterative procedure that allows updating the solute charge distribution in the presence of the solvent is combined with a sequential Monte Carlo/Quantum Mechanics methodology and Density Functional Theory methods to evaluate the nonlinear optical (NLO) response using the hyper Rayleigh scattering (HRS) of a series of chromones recently identified in Chamaecrista diphylla, an herbaceous plant abundant throughout the Americas and used in folk medicine. From this study, it is determined that from gas to solvent environment, the systems acquire low refractive index (n) and an improvement of the first hyperpolarizability (β ), signaling potential NLO uses.

Free-Energy Landscape of the S2 Reaction CHBr + Cl → CHCl + Br in Different Liquid Environments.

This work describes in detail the reaction path of the well-known S2 reaction CHBr + Cl → CHCl + Br, whose reaction rate has a huge variation with the solvent in the gas phase and in protic and aprotic liquid environments. We employed the ASEC-FEG method to optimize for minima (reactants and products) and saddle points (transition states) in the in-solution free-energy hypersurface. The method takes atomistic details of the solvent into account.

Free Energy Computation for an Isomerizing Chromophore in a Molecular Cavity via the Average Solvent Electrostatic Configuration Model: Applications in Rhodopsin and Rhodopsin-Mimicking Systems.

We present a novel technique for computing the free energy differences between two chromophore "isomers" hosted in a molecular environment (a generalized solvent). Such an environment may range from a relatively rigid protein cavity to a flexible solvent environment. The technique is characterized by the application of the previously reported "average electrostatic solvent configuration" method, and it is based on the idea of using the free energy perturbation theory along with a chromophore annihilation procedure in thermodynamic cycle calculations.

On the population of triplet states of 2-seleno-thymine.

The population and depopulation mechanisms leading to the lowest-lying triplet states of 2-Se-Thymine were studied at the MS-CASPT2/cc-pVDZ level of theory. Several critical points on different potential energy hypersurfaces were optimized, including minima, conical intersections, and singlet-triplet crossings. The accessibility of all relevant regions on the potential energy hypersurfaces was investigated by means of minimum energy paths and linear interpolation in internal coordinates techniques.

A new interpretation of the absorption and the dual fluorescence of Prodan in solution.

Remarkable interest is associated with the interpretation of the Prodan fluorescent spectrum. A sequential hybrid Quantum Mechanics/Molecular Mechanics method was used to establish that the fluorescent emission occurs from two different excited states, resulting in a broad asymmetric emission spectrum. The absorption spectra in several solvents were measured and calculated using different theoretical models presenting excellent agreement.

Solvation Structures and Deactivation Pathways of Luminescent Isothiazole-Derived Nucleobases: , , and .

The photophysical relaxation pathways of , , and luminescent nucleobases were investigated with the MS-CASPT2 quantum-chemical method and double-ζ basis sets (cc-pVDZ) in gas and condensed phases (1,4-dioxane and water) with the sequential Monte Carlo/CASPT2 and free energy gradient (FEG) methods. Solvation shell structures, in the ground and excited states, were examined with the pairwise radial distribution function (()) and solute-solvent hydrogen-bond networks. Site-specific hydrogen bonding analysis evidenced relevant changes between both electronic states.

DICE: A Monte Carlo Code for Molecular Simulation Including the Configurational Bias Monte Carlo Method.

Solute-solvent systems are an important topic of study, as the effects of the solvent on the solute can drastically change its properties. Theoretical studies of these systems are done with methods, molecular simulations, or a combination of both. The simulations of molecular systems are usually performed with either molecular dynamics (MD) or Monte Carlo (MC) methods.

Unraveling the Electric Field-Induced Second Harmonic Generation Responses of Stilbazolium Ion Pairs Complexes in Solution Using a Multiscale Simulation Method.

The electric field-induced second harmonic generation (EFISHG) response has been largely used to describe the first β and the second γ hyperpolarizabilities in solution. Although the EFISHG technique cannot be applied to charged compounds (due to the external static electric field), it can be used to describe ion pairs as neutral complexes. A multiscale computational approach is required to generate representative geometrical configurations of such kinds of complexes (using classical force fields), to compute the electronic structure of each configuration (using quantum mechanics methods), and to perform statistical analyses describing the behavior of the nonlinear optical properties.

Microscopic Origin of Different Hydration Patterns of para-Nitrophenol and Its Anion: A Study Combining Multiconfigurational Calculations and the Free-Energy Gradient Method.

A theoretical study of the solvatochromic shifts of para-nitrophenol ( pNP) and para-nitrophenolate anion ( pNP) in aqueous solution is presented using a QM/MM methodology with molecular dynamics simulation. The optimized structures in aqueous solution are obtained using both the polarizable continuum and the free-energy gradient methods. For pNP, the calculated redshifts at the CASPT2 (12,10) level are, respectively, 0.

Theoretical study of the NMR chemical shift of Xe in supercritical condition.

In this work we investigate the level of theory necessary for reproducing the non-linear variation of the Xe nuclear magnetic resonance (NMR) chemical shift with the density of Xe in supercritical conditions. In detail we study how the Xe chemical shift depends under supercritical conditions on electron correlation, relativistic and many-body effects. The latter are included using a sequential-QM/MM methodology, in which a classical MD simulation is performed first and the chemical shift is then obtained as an average of quantum calculations of 250 MD snapshots conformations carried out for Xe clusters (n = 2 - 8 depending on the density).

An Average Solvent Electrostatic Configuration Protocol for QM/MM Free Energy Optimization: Implementation and Application to Rhodopsin Systems.

A novel atomistic methodology to perform free energy geometry optimization of a retinal chromophore covalently bound to any rhodopsin-like protein cavity is presented and benchmarked by computing the absorption maxima wavelengths (λ) of distant rhodopsin systems. The optimization is achieved by computing the Nagaoka's Free Energy Gradient (FEG) within an Average Solvent Electrostatic Configuration (ASEC) atomistic representation of the thermodynamic equilibrium and minimizing such quantity via an iterative procedure based on sequential classical MD and constrained QM/MM geometry optimization steps. The performance of such an ASEC-FEG protocol is assessed at the CASPT2//CASSCF/Amber level by reproducing the λ values observed for 12 mutants of redesigned human cellular retinol binding protein II (hCRBPII) systems; a set of 10 distant wild-type rhodopsins from vertebrates, invertebrates, eubacteria, and archaea organisms; and finally a set of 10 rhodopsin mutants from an eubacterial rhodopsin.

Electronic structure and absorption spectra of fluorescent nucleoside analogues.

This work presents a systematic investigation of the electronic and conformational properties of five new fluorescent nucleobases belonging to the alphabet based on the isothiazole[4,3-d]pyrimidine molecule, very recently synthesized. This is of particular importance in the characterization of the main electronic aspects of these fluorescent nucleosides. The solvent effects of 1,4-dioxane and water were included combining the Sequential Monte Carlo/CASPT2 and the Free Energy Gradient (FEG) methods.

A First-Principles Approach to the Dynamics and Electronic Properties of p-Nitroaniline in Water.

Born-Oppenheimer molecular dynamics of p-nitroaniline (PNA) in water was carried out and the electronic structure was investigated by time-dependent density functional theory. Hydrogen bonding involving the PNA nitro and amine groups and the water molecules leads to an ∼160 cm(-1) red shift of the ν(N-O) and ν(N-H) stretching frequencies relative to the gas phase species. Our estimate for the peak position of the charge transfer (CT) band in the absorption spectrum of PNA in water (∼3.

Solvent Effects on the Dynamic Polarizability and Raman Response of Molecule-Metal Oxide Hybrid Clusters.

Currently, there is considerable interest in the properties of semiconducting metal oxide nanoparticle substrates because of their utility in surface-enhanced Raman scattering, dye-sensitized solar cells, and photocatalysis. While the enhancement of Raman activities of molecules adsorbed on these nanoparticles is due to a large increase in the polarizability, because of charge transfer from the molecule to the semiconducting nanoparticle, little is known about the factors responsible for modulating the polarizability, particularly the influence of the solvent. Consequently, we have carried out Monte Carlo simulations of several hybrids to study the solvent effect on the dynamic polarizabilities and electronic spectra.

Probing Lewis Acid-Base Interactions with Born-Oppenheimer Molecular Dynamics: The Electronic Absorption Spectrum of p-Nitroaniline in Supercritical CO2.

The structure and dynamics of p-nitroaniline (PNA) in supercritical CO2 (scCO2) at T = 315 K and ρ = 0.81 g cm(-3) are investigated by carrying out Born-Oppenheimer molecular dynamics, and the electronic absorption spectrum in scCO2 is determined by time dependent density functional theory. The structure of the PNA-scCO2 solution illustrates the role played by Lewis acid-base (LA-LB) interactions.

Behavior of the dielectric constant of Ar near the critical point.

The fundamental question of the behavior of the dielectric constant near the critical point is addressed using Ar as the probe system. The neighborhood of the liquid-vapor critical point of Ar is accessed by classical Monte Carlo simulation and then explicit quantum mechanics calculations are performed to study the behavior of the dielectric constant. The theoretical critical temperature is determined by calculating the position of the discontinuity of the specific heat and is found to be at T(c)Theor=148.