Calculation of macroscopic first-, second-, and third-order optical susceptibilities for the urea crystal

A Multi-Scale Approach to Simulate the Nonlinear Optical Response of Molecular Nanomaterials

Nonlinear optics is essential for many recent photonic technologies. Here, a novel multi-scale approach is introduced to simulate the nonlinear optical response of molecular nanomaterials combining ab initio quantum-chemical and classical Maxwell-scattering computations. In this approach, the first hyperpolarizability tensor is computed with time-dependent density-functional theory and incorporated into a multi-scattering formalism that considers the optical interaction between neighboring molecules. Such incorporation is achieved by a novel object: the Hyper-Transition(T)-matrix. With this object at hand, the nonlinear optical response from single molecules and also from entire photonic devices can be computed, including the full tensorial and dispersive nature of the optical response of the molecules, as well as the optical interaction between different molecules as, for example, in the lattice of a molecular crystal. To demonstrate the applicability of the novel approach, the generation of a second-harmonic signal from a thin film of an Urea molecular crystal is computed and compared to more traditional simulations. Furthermore, an optical cavity is designed, which enhances the second-harmonic response of the molecular film up to more than two orders of magnitude. This approach is highly versatile and accurate and can be the working horse for the future exploration of nonlinear photonic molecular materials in structured photonic environments.

Theoretical investigation on the linear and nonlinear optical properties of DAPSH crystal

The linear polarizability, first and second hyperpolarizabilities of the asymmetric unit of DAPSH crystal are studied and compared with available experimental results. The polarization effects are included using an iterative polarization procedure, which ensures the convergence of the dipole moment of DAPSH embedded within a polarization field generated by the surrounding asymmetric units whose atomic sites are considered as point charges. We estimate macroscopic susceptibilities from the results of the polarized asymmetric units in the unit cell, considering the significant contribution of the electrostatic interactions in crystal packing. The results show that the influence of the polarization effects leads to a marked decrease of the first hyperpolarizability, compared with the respective isolated counterpart, which improves the concordance with the experiment. There is a minor influence of polarization effects on the second hyperpolarizability but our estimated result for the third-order susceptibility, related to the NLO process of the intensity dependent refractive index, is significant as compared with the results for other organic crystals, such as chalcone-derivatives. In addition, supermolecule calculations are conducted for explicit dimers in presence of the electrostatic embedding to illustrate the role played by the electrostatic interactions in the hyperpolarizabilities of the DAPSH crystal.

Molecular modeling and nonlinear optical properties of new isostructural halogenated dihydroquinolinones

Two new isostructural halogenated dihydroquinolinones were synthesized and characterized by single crystal X-ray diffraction.

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.

An Efficient Synthesis, Spectroscopic Characterization, and Optical Nonlinearity Response of Novel Salicylaldehyde Thiosemicarbazone Derivatives

In this study, seven derivatives of salicylaldehyde thiosemicarbazones (1-7) were synthesized by refluxing substituted thiosemicarbazide and salicylaldehyde in an ethanol solvent. Different spectral techniques (UV-vis, IR, and NMR) were used to analyze the prepared compounds (1-7). Accompanied by the experimental study, quantum chemical studies were also carried out at the M06/6-311G(d,p) level. A comparative analysis of the UV-visible spectra and vibrational frequencies between computational and experimental findings was also performed. These comparative data disclosed that both studies were observed to be in excellent agreement. Furthermore, natural bond orbital investigations revealed that nonbonding transitions were significant for the stability of prepared molecules. In addition, frontier molecular orbital (FMO) findings described that a promising charge transfer phenomenon was found in 1-7. The energies of FMOs were further used to determine global reactivity parameters (GRPs). These GRP factors revealed that all synthesized compounds (1-7) contain a greater hardness value (η = 2.1 eV) and a lower softness value (σ = 0.24 eV), which indicated that these compounds were less reactive and more stable. Nonlinear optical (NLO) evaluation displayed that compound 5 consisted of greater values of linear polarizability ⟨α⟩ and third-order polarizability ⟨γ⟩ of 324.93 and 1.69 × 105 a.u., respectively, while compound 3 exhibited a larger value of second-order polarizability (βtotal) of 508.41 a.u. The NLO behavior of these prepared compounds may be significant for the hi-tech NLO applications.

2D materials towards ultrafast photonic applications

Two-dimensional materials are now excelling in yet another arena of ultrafast photonics, including optical modulation through optical limiting/mode-locking, photodetectors, optical communications, integrated miniaturized all-optical devices, etc.

First principles crystal engineering of nonlinear optical materials. I. Prototypical case of urea

The crystalline materials with nonlinear optical (NLO) properties are critically important for several technological applications, including nanophotonic and second harmonic generation devices. Urea is often considered to be a standard NLO material, due to the combination of non-centrosymmetric crystal packing and capacity for intramolecular charge transfer. Various approaches to crystal engineering of non-centrosymmetric molecular materials were reported in the literature. Here we propose using global lattice energy minimization to predict the crystal packing from the first principles. We developed a methodology that includes the following: (1) parameter derivation for polarizable force field AMOEBA; (2) local minimizations of crystal structures with these parameters, combined with the evolutionary algorithm for a global minimum search, implemented in program USPEX; (3) filtering out duplicate polymorphs produced; (4) reoptimization and final ranking based on density functional theory (DFT) with many-body dispersion (MBD) correction; and (5) prediction of the second-order susceptibility tensor by finite field approach. This methodology was applied to predict virtual urea polymorphs. After filtering based on packing similarity, only two distinct packing modes were predicted: one experimental and one hypothetical. DFT + MBD ranking established non-centrosymmetric crystal packing as the global minimum, in agreement with the experiment. Finite field approach was used to predict nonlinear susceptibility, and H-bonding was found to account for a 2.5-fold increase in molecular hyperpolarizability to the bulk value.

Synthesis, spectroscopic, dielectric, molecular docking and DFT studies of (3E)-3-(4-methylbenzylidene)-3,4-dihydro-2H-chromen-2-one: an anticancer agent

BACKGROUND

Coumarin (2H-chromen-2-one) and its derivatives have a wide range of biological and pharmaceutical activities. They possess antitumor, anti-HIV, anticoagulant, antimicrobial, antioxidant, and anti-inflammatory activities. Synthesis and isolation of coumarins from different species have attracted the attention of medicinal chemists. Herein, we report the synthesis, molecular structure, dielectric, anticancer activity and docking studies with the potential target protein tankyrase.

RESULTS

Molecular structure of (3E)-3-(4-methylbenzylidene)-3,4-dihydro-2H-chromen-2-one (MBDC) is derived from quantum chemical calculations and compared with the experimental results. Intramolecular interactions, stabilization energies, and charge delocalization are calculated by NBO analysis. NLO property and dielectric quantities have also been determined. It indicates the formation of a hydrogen bonding between -OH group of alcohol and C=O of coumarin. The relaxation time increases with the increase of bond length confirming the degree of cooperation and depends upon the shape and size of the molecules. The molecule under study has shown good anticancer activity against MCF-7 and HT-29 cell lines. Molecular docking studies indicate that the MBDC binds with protein.

CONCLUSIONS

In this study, the compound (3E)-3-(4-methylbenzylidene)-3,4-dihydro-2H-chromen-2-one was synthesized and characterized by spectroscopic studies. The computed and experimental results of NMR study are tabulated. The dielectric relaxation studies show the existence of molecular interactions between MBDC and alcohol. Theoretical results of MBDC molecules provide the way to predict various binding sites through molecular modeling and these results also support that the chromen substitution is more active in the entire molecule. Molecular docking study shows that MBDC binds well in the active site of tankyrase and interact with the amino acid residues. These results are compared with the anti cancer drug molecule warfarin derivative. The results suggest that both molecules have comparable interactions and better docking scores. The results of the antiproliferative activity of MBDC and Warfarin derivative against MCF-7 breast cancer and HT-29 colon cancer cell lines at different concentrations exhibited significant cytotoxicity. The estimated half maximal inhibitory concentration (IC 50) value for MBDC and Warfarin derivative was 15.6 and 31.2 μg/ml, respectively. This enhanced cytotoxicity of MBDC in MCF-7 breast cancer and HT-29 colon cancer cell lines may be due to their efficient targeted binding and eventual uptake by the cells. Hence the compound MBDC may be considered as a drug molecule for cancer.Graphical abstractThe binding mode of the ligand MBDC at active site of protein and the graphical representation of cell inhibition for MCF-7 and HT-29 cell lines.

Polarization effects on the electric properties of urea and thiourea molecules in solid phase

We present theoretical results for the dipole moment, linear polarizability, and first hyperpolarizability of the urea and thiourea molecules in solid phase. The in-crystal electric properties were determined by applying a supermolecule approach in combination with an iterative electrostatic scheme, in which the surrounding molecules are represented by point charges. It is found for both urea and thiourea molecules that the influence of the polarization effects is mild for the linear polarizability, but it is marked for the dipole moment and first hyperpolarizability. The replacement of oxygen atoms by sulfur atoms increases, in general, the electric responses. Our second-order Møller-Plesset perturbation theory based iterative scheme predicts for the in-crystal dipole moment of urea and thiourea the values of 7.54 and 9.19 D which are, respectively, increased by 61% and 58%, in comparison with the corresponding isolated values. The result for urea is in agreement with the available experimental result of 6.56 D. In addition, we present an estimate of macroscopic quantities considering explicit unit cells of urea and thiourea crystals including environment polarization effects. These supermolecule calculations take into account partially the exchange and dispersion effects. The results illustrate the role played by the electrostatic interactions on the static second-order nonlinear susceptibility of the urea crystal.

Which charge definition for describing the crystal polarizing field and the χ(1) and χ(2) of organic crystals?

Crystal optical susceptibilities are probes to assess the performance of the charge definition employed to describe the crystal polarizing field.

A theoretical investigation of electric properties of L-arginine phosphate monohydrate including environment polarization effects

The dipole moment (μ), linear polarizability (α), and first hyperpolarizability (β(tot)) of the asymmetric unit of L-arginine phosphate (LAP) monohydrate crystal are investigated using the supermolecule approach in combination with an iterative electrostatic polarization scheme. Environment polarization effects are attained by assuring the convergence of the dipole moment of LAP embedded in the polarization field of the surrounding molecules whose atomic sites are treated as point charges. The results obtained show that in the presence of the embedding charges, the value of μ is increased by 9% but the static values of α and β(tot) are decreased, respectively, by 3% and 13%, as compared with the isolated situation. The MP2/6-311+G(d) model predicts for the in-crystal dipole moment the converged value of 33 D, in good concordance with the available experimental result of 32 D. Our estimates for the converged results of α and β(tot) are, respectively, 22.51×10(-24) and 5.01×10(-30) esu. Dispersion effects are found to have a small impact on the nonlinear optical responses of LAP in the visible region. In addition, MP2/6-311G results obtained for β(tot) by using isolated and embedded LAP dimers show that crystal packing effects have a significant contribution of the electrostatic interactions. Our results suggest that the role of the crystal environment is to minimize the effects of the intermolecular interactions in the electric properties. That is, μ and β(tot) gain a more additive character in the presence of the field of the embedding charges. This is specially marked for β(tot).

Self-assembly from the gas-phase: design and implementation of small-molecule chromophore precursors with large nonlinear optical responses

Efficiently organizing molecular nonlinear optical (NLO) chromophores having large first-order hyperpolarizabilities (beta) into acentric microstructures for electro-optic (EO) applications represents a significant materials synthesis and processing challenge, in part due to interchromophore dipolar interactions that promote centrosymmetric organization. Here we report the computational modeling, synthesis, and characterization of a series of eight heteroaromatic organic chromophores, designed to self-organize from the vapor phase via directed hydrogen-bond networks, into acentric thin films. Introduction of alpha,omega-donor-acceptor hydrogen-bonding substituents along the molecular long axes tunes properties such as hyperpolarizability, volatility, thermal stability, film-forming properties, and macroscopic NLO response (chi((2))). DFT-level molecular modeling, INDO/S optical property analysis, and sum-overstates computation indicate that molecular-core fluorination and hydrogen-bond donor incorporation can increase beta(vec) up to 40x versus that of typical fluorine-free chromophores. Furthermore, inclusion of sterically induced biphenyl conjugative decoupling between chromophore pi-donor substituents and the hydrogen-bonding donor sites increases beta by approximately 50%. Experimental thin-film second harmonic generation (SHG) spectroscopy confirms these trends in calculated responses, with chi((2)) increasing 7.5x upon chromophore core fluorination and 15x with hydrogen-bonding donor substitution, thereby achieving macroscopic responses as high as 302 pm/V at omega(o) = 1064 nm. In addition to response trends, cluster calculations also reveal linear additivity in beta(vec) with catenation for all benzoic acid-containing chromophores up to longitudinally aligned trimers. Linear scaling of SHG response with film thickness is observed for benzoic acid-containing chromophores up to 1.0 microm film thickness.

Refractive indices for molecular crystals from the response of X-ray constrained Hartree-Fock wavefunctions

Refractive indices for molecular crystals are obtained from Hartree-Fock wavefunctions constrained to reproduce a set of experimental X-ray structure factors. Coupled-perturbed Hartree-Fock theory is used to calculate the in-crystal effective polarizabilities from which the refractive indices are obtained, thus eliminating the need for the calibration procedure used in earlier work by Whitten et al. [J. Chem. Phys. 2006, 125, 174505]. The results clearly demonstrate that these X-ray constrained Hartree-Fock (XCHF) wavefunctions reflect genuine effects of intermolecular interactions in crystals. Molecular dipole moments are consistently in excellent agreement with ab initio MP2 estimates that incorporate the effects of the crystal field. Consistent agreement of the XCHF refractive indices with experimental measurements at optical frequencies confirms that this approach can provide both meaningful results and considerable insight into the relative importance of molecular properties and crystal field effects in determining the detailed nature of the refractivity tensor.

Calculation of the Microscopic and Macroscopic Linear and Nonlinear Optical Properties of Liquid Acetonitrile. II. Local Fields and Linear and Nonlinear Susceptibilities in Quadrupolar Approximation

A discrete model based on the multipolar expansion including terms up to hexadecapoles was employed to describe the electrostatic interactions in liquid acetonitrile. Liquid structures obtained form molecular dynamics simulations with different classical, nonpolarizable potentials were used to analyze the electrostatic interactions. The computed average local field was employed for the determination of the environmental effects on the linear and nonlinear electrical molecular properties. Dipole-dipole interactions yield the dominant contribution to the local field, whereas higher multipolar contributions are small but not negligible. Using the effective in-phase properties, macroscopic linear and nonlinear susceptibilities of the liquid were computed. Depending on the partial charges describing the Coulomb interactions of the force field employed, either the linear properties (refractive index and dielectric constant) were reproduced in good agreement with experiment or the nonlinear properties [third-harmonic generation (THG) and electric field induced second-harmonic (EFISH) generation] and the bulk density but never both sets of properties together. It is concluded that the partial charges of the force fields investigated are not suitable for reliable dielectric properties. New methods are probably necessary for the determination of partial charges, which should take into account the collective and long-range nature of electrostatic interactions more precisely.

Effective molecular polarizabilities and crystal refractive indices estimated from x-ray diffraction data

Although it was proposed some time ago that (hyper)polarizabilities might be estimated from the results of x-ray charge density refinements, early results were unconvincing. In this work we show that the one particle density obtained from the usual multipole refinement model does not contain sufficient information to determine these response properties and instead pursue the "constrained wave function" approach of fitting to x-ray structure factors. Simplified sum-over-states expressions are derived for determining the dipole polarizability from these wave functions, and these clearly show that the earlier work ignored important two-electron expectation values for the dipole polarizability, and two- and three-electron terms for Beta, etc. Correction factors for the simplified sum-over-states polarizability tensors from the constrained wave function are obtained by calibration against coupled Hartree-Fock ab initio results to yield in-crystal effective polarizability tensors. Results obtained for benzene, urea, and 2-methyl-4-nitroaniline demonstrate that the effective molecular polarizabilities clearly include the effects of intermolecular interactions and electron correlation, especially for urea where the effects on the polarizability are known to be quite large. We also carefully consider the way in which the linear bulk susceptibility, chi((1)), and refractive indices are determined from the x-ray fitted polarizabilities, employing three models based on a rigorous treatment of the local field. Incorrect results are obtained for the sort of molecules that are of interest in nonlinear optical applications if the molecules are approximated by single point dipoles. In contrast, the use of Lorentz-factor tensors averaged over several sites yields excellent results, with refractive indices obtained using this model in remarkably good agreement with optical measurements extrapolated to zero frequency.

Computer Simulation of the Linear and Nonlinear Optical Susceptibilities of p-Nitroaniline in Cyclohexane, 1,4-Dioxane, and Tetrahydrofuran in Quadrupolar Approximation. II. Local Field Effects and Optical Susceptibilitities

This is the second part of a study to elucidate the local field effects on the nonlinear optical properties of p-nitroaniline (pNA) in three solvents of different multipolar character, that is, cyclohexane (CH), 1,4-dioxane (DI), and tetrahydrofuran (THF), employing a discrete description of the solutions. By the use of liquid structure information from molecular dynamics simulations and molecular properties computed by high-level ab initio methods, the local field and local field gradients on p-nitroaniline and the solvent molecules are computed in quadrupolar approximation. To validate the simulations and the induction model, static and dynamic (non)linear properties of the pure solvents are also computed. With the exception of the static dielectric constant of pure THF, a good agreement between computed and experimental refractive indices, dielectric constants, and third harmonic generation signals is obtained for the solvents. For the solutions, it is found that multipole moments up to two orders higher than quadrupole have a negligible influence on the local fields on pNA, if a simple distribution model is employed for the electric properties of pNA. Quadrupole effects are found to be nonnegligible in all three solvents but are especially pronounced in the 1,4-dioxane solvent, in which the local fields are similar to those in THF, although the dielectric constant of DI is 2.2 and that of the simulated THF is 5.4. The electric-field-induced second harmonic generation (EFISH) signal and the hyper-Rayleigh scattering signal of pNA in the solutions computed with the local field are in good to fair agreement with available experimental results. This confirms the effect of the "dioxane anomaly" also on nonlinear optical properties. Predictions based on an ellipsoidal Onsager model as applied by experimentalists are in very good agreement with the discrete model predictions. This is in contrast to a recent discrete reaction field calculation of pNA in 1,4-dioxane, which found that the predicted first hyperpolarizability of pNA deviated strongly from the predictions obtained using Onsager-Lorentz local field factors.

Computer Simulation of the Linear and Nonlinear Optical Susceptibilities of p-Nitroaniline in Cyclohexane, 1,4-Dioxane, and Tetrahydrofuran in Quadrupolar Approximation. I. Molecular Polarizabilities and Hyperpolarizabilities

This is the first part of a study of the local field effects on (non)linear optical susceptibilities of solutions of para-nitroaniline (pNA) in three different solvents, cyclohexane (CH), 1,4-dioxane (DI), and tetrahydrofuran (THF), using a discrete molecular representation of the condensed phase. To account for dipolar and quadrupolar effects, the latter of which are especially important for DI solution, all the electric properties necessary to compute the local fields and local field gradients in quadrupolar approximation as well as the dipolar hyperpolarizabilities for the four molecules are computed, including frequency dispersion and vibrational contributions to the dipolar properties. The convergence of the perturbation treatment for the pure vibrational (PV) contributions is examined by comparison of the values obtained at the lowest order with those of partially computed second order in mechanical and electrical anharmonicity. For pNA, for which previous computations of the hyperpolarizabilities have generally found poor agreement with experimental results, a thorough investigation of the effects of solvent-induced geometry changes, dynamic and static correlation, frequency dispersion, and classical thermal averaging over the torsional modes of the substituent groups and the inversion mode of the amino group on the dipolar properties is carried out. Computations using self-consistent continuum reaction field models show that the amino group is substantially less pyramidalized in polar solvents than in the gas phase. With all the effects taken into account, reasonable agreement with the experimental electric-field induced second harmonic generation (EFISH) result on pNA vapor of Kaatz, Donley, and Shelton is obtained.

Theoretical investigation of the linear and second-order nonlinear susceptibilities of the 3-methyl-4-nitropyridine-1-oxyde (POM) crystal

The linear and nonlinear optical properties of the 3-methyl-4-nitropyridine-1-oxyde (POM) crystal have been evaluated using semiempirical quantum chemistry techniques. The scheme includes (i) the evaluation of the polarizability and first hyperpolarizability of increasingly large one-dimensional, two-dimensional, and three-dimensional clusters of POM, (ii) the use of the time-dependent Hartree-Fock approach to determine the static and dynamic responses in combination with semiempirical Austin model 1 Hamiltonian, (iii) the assessment, for the POM monomer and dimer, of the electron correlation effects using second-order Moller-Plesset perturbation theory with several basis sets containing polarization and diffuse functions, (iv) the assessment of the validity of the multiplicative scheme and its application to get effective polarizability and first hyperpolarizability of the POM unit cell in the crystal, (v) the use of a sum-over-states approach to attribute the first hyperpolarizability to a dominant charge-transfer excited state, and (vi) comparison with experimental data as well as with calculated values obtained using the oriented gas approximation.

The Dipole Moment, Polarizabilities, and First Hyperpolarizabilities of HArF. A Computational and Comparative Study

The electronic and vibrational contributions to the dipole moment, polarizabilities, and first hyperpolarizabilities of HArF are reported. These have been computed by using a series of systematically built basis sets and a hierarchy of computational methods. HArF has a very large first hyperpolarizability along the z axis. This has been rationalized by invoking the difference in the electronic structure between the ground and the first excited state. The argon fluorohydride has been recently derived and characterized. The present study provides complementary data for the understanding of the electronic structure of this interesting argon derivative.