Effect of acceptor moieties on static and dynamic first hyperpolarizability of azobenzene chromophores

Multimillion Atom Simulations of Di-8-ANEPPS Chromophores Embedded in a Model Plasma Membrane: Toward the Investigation of Realistic Dyed Cell Membranes

A multistep computational approach has been employed to study a multimillion all-atom dyed plasma membrane, with no less than 42 different lipid species spanning the major head groups and a variety of fatty acids, as well as cholesterol, with the objective of investigating its structure and dynamics, as well as its impact on the embedded di-8-ANEPPS dyes. The latter are commonly used as bioimaging probes and serve as local microscopes. So, they provide information on membrane morphology via their second harmonic nonlinear optical (NLO) responses, which have the advantage of being specific to interface regions and sensitive to the chromophore environment. In previous studies, this chromophore has only been studied in simpler membrane models, far from the complexity of real lipid bilayers, while, owing to the ever-increasing computational resources, multimillion lipid bilayers have been studied, giving access to the effects of its heterogeneity. First, using molecular dynamics (MD) simulations, it is found that the combination of lipids produces a more ordered and denser membrane compared to its homogeneous model counterparts, while the local environment of the embedded dyes becomes enriched in phosphatidylcholine. Subsequently, the second harmonic first hyperpolarizability of the probes was calculated at the TDDFT level on selected frames of MD, highlighting the influence of the lipid environment. Due to the complexity of the system, machine learning (ML) tools have been employed to establish relationships between the membrane structural parameters, the orientation of the probes, and their NLO responses. These ML approaches have revealed influential features, including the presence of diacylglycerol lipids close to the dye. On the whole, this work provides a first step toward understanding the cooperation, synergy, and interactions that occur in such complex guest-host environments, which have emerged as new targets for drug design and membrane lipid therapy.

First hyperpolarizability of the di-8-ANEPPS and DR1 nonlinear optical chromophores in solution. An experimental and multi-scale theoretical chemistry study

The solvent effects on the linear and second-order nonlinear optical properties of an aminonaphtylethenylpyridinium (ANEP) dye are investigated by combining experimental and theoretical chemistry methods. On the one hand, deep near infrared (NIR) hyper-Rayleigh scattering (HRS) measurements (1840-1950 nm) are performed on solutions of di-8-ANEPPS in deuterated chloroform, dimethylformamide, and dimethylsulfoxide to determine their first hyperpolarizablity (βHRS). For the first time, these HRS experiments are carried out in the picosecond regime in the deep NIR with very moderate (≤3 mW) average input power, providing a good signal-to-noise ratio and avoiding solvent thermal effects. Moreover, the frequency dispersion of βHRS is investigated for Disperse Red 1 (DR1), a dye commonly used as HRS external reference. On the other hand, these are compared with computational chemistry results obtained by using a sequential molecular dynamics (MD) then quantum mechanics (QM) approach. The MD method allows accounting for the dynamical nature of the molecular structures. Then, the QM part is based on TDDFT/M06-2X/6-311+G* calculations using solvation models ranging from continuum to discrete ones. Measurements report a decrease of the βHRS of di-8-ANEPPS in more polar solvents and these effects are reproduced by the different solvation models. For di-8-ANEPPS and DR1, comparisons show that the use of a hybrid solvation model, combining the description of the solvent molecules around the probe by point charges with a continuum model, already achieves quasi quantitative agreement with experiment. These results are further improved by using a polarizable embedding that includes the atomic polarizabilities in the solvent description.

Unravelling the Effects of Cholesterol on the Second-Order Nonlinear Optical Responses of Di-8-ANEPPS Dye Embedded in Phosphatidylcholine Lipid Bilayers

Cholesterol is known for its role in maintaining the correct fluidity and rigidity of the animals cell membranes and thus their functions. Assessing the content and the role of cholesterol in lipid bilayers is therefore of crucial importance for a deeper understanding and control of membrane functioning. In this computational work, we investigate bilayers built from three types of glycerophospholipid phosphatidylcholine (PC) lipids, namely dipalmitoylphosphatidylcholine (DPPC), 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC), and dioleoylphosphatidylcholine (DOPC), and containing different amounts of cholesterol by analyzing the second-harmonic generation (SHG) nonlinear optical (NLO) response of a probe molecule, di-8-ANEPPS, inserted into the membranes. This molecular property presents the advantage to be specific to interfacial regions such as lipid bilayers. To unravel these effects, Molecular Dynamics (MD) simulations have been performed on both DPPC and DOPC lipids by varying the cholesterol mole fraction (from 0 to 0.66), while POPC was only considered as a pure bilayer. In the case of the structural properties of the bilayers, all the analyses converge toward the same conclusion: as the mole fraction of cholesterol increases, the systems become more rigid, confirming the condensing effect of cholesterol. In addition, the chromophore is progressively more aligned with respect to the normal to the bilayer. On the contrary, addition of unsaturation disorders the lipid bilayers, with barely no impact on the alignment of the chromophore. Then, using the frames obtained from the MD simulations, the first hyperpolarizability β of the dye in its environment has been computed at the TDDFT level. On the one hand, the addition of cholesterol induces a progressive increase of the diagonal component the β tensor parallel to the bilayer normal. On the other hand, larger β values have been calculated for the unsaturated than for the saturated lipid systems. In summary, this study illustrates the relationship between the composition and structure of the bilayers and the NLO responses of the embedded dye.

Theoretical study by DFT of organometallic complexes based on metallocenes active in NLO

This study is based on the valuation of a few model molecules. The objective of this research is focused on linear optical (LO) and nonlinear optical (NLO) enhancement of five organometallic molecules based on different metallocenes. These molecules were subjected to several calculations by different long-range functionals CAM-B3LYP, LC-BLYP, LC-wPBE, wB97X, M11, and the following three Minnesota functionals: M06-2X and M08-HX in comparison with the MP2 approach. Hence, the CAM-B3LYP functional gave the closest NLO values to the MP2 method. Second, molecule 3A based on nickelocene recorded the highest static (β_tot) value which is 76.46 Χ 10^-30 esu and 4803.4 Χ 10^-30 esu under the laser wavelength λ = 532 nm. Third, intramolecular charge transfers (ICTs) of the molecules studied are all directed in both directions (donor to acceptor and vice versa). Finally, the specific solvent for molecules 2A and 3A is acetonitrile, and the maximum wavelengths obtained for the isolated or solvated molecule are all located in the near UV; the corresponding interval is between 250 and 395 nm. Graphical abstract.

Benchmark study of the linear and nonlinear optical polarizabilities in proto-type NLO molecule of para-nitroaniline

In the present investigation, for the first time, we have performed a thorough study about different functionals and basis sets for linear and nonlinear optical (NLO) properties of para-nitroaniline ([Formula: see text]-NA), which is considered as proto-type NLO molecule, among organic NLO materials. There is a dire need of such data base for [Formula: see text]-NA because many investigators are using such values of [Formula: see text]-NA for comparative analysis. A range of different functionals including HF, BLYP, PW91, PBE, B3LYP, M06, M06-2X, PBE0, BHandHLYP, CAM-B3LYP, LC-BLYP, and B3LYP-D3 are applied in conjugation with several commonly basis sets such as 6-31G*, 6-311G*, 6-311G**, 6-311+G**, cc-pVDZ, and cc-pVTZ. A variety of functional and basis sets combinations are calculated and graphically compared with each other. The calculated total dipole moment for the [Formula: see text]-NA is found to be 6.79[Formula: see text]D which is quite closer to experimentally determined value. The lowest calculated value for linear isotropic polarizability at HF/6-31G* level of theory is [Formula: see text] esu while higher values observed with remaining all methods especially 14% polarizability increases in presence of basis set with diffuse functions and similar trend of variation is also observed in linear anisotropic polarizability. Similarly, the calculated value of frequency dependent second-order polarizability is found to be [Formula: see text] esu at PBE0/6-311+G** level of theory which is quite closer to experimental value of [Formula: see text] esu. A comparison between the calculated and experimental results shows good agreement among geometries, dipole moments and NLO polarizabilities for [Formula: see text]-NA. Moreover, the frontier molecular orbital (FMO) and electron density difference map (EDDM) analysis along with density of state (DOS) plots are also presented to get more physical intuitions into the structure–property relationship and electronic communications between terminal accepter and donor groups through [Formula: see text]-conjugation. The present investigation provides benchmark data including various commonly used functionals and basis sets for the calculation of NLO properties of [Formula: see text]-NA. Thus, the present investigation will put straight several future studies when it comes to comparative NLO study of organic materials.