1
|
Omodemi O, Kaledin M, Kaledin AL. Permutationally invariant polynomial representation of polarizability tensor surfaces for linear regression analysis. J Comput Chem 2022; 43:1495-1503. [PMID: 35737590 DOI: 10.1002/jcc.26952] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/20/2022] [Accepted: 06/09/2022] [Indexed: 11/07/2022]
Abstract
A linearly parameterized functional form for a Cartesian representation of molecular dipole polarizability tensor surfaces (PTS) is described. The proposed expression for the PTS is a linearization of the recently reported power series ansatz of the original Applequist model, which by construction is non-linear in parameter space. This new approach possesses (i) a unique solution to the least-squares fitting problem; (ii) a low level of the computational complexity of the resulting linear regression procedure, comparable to those of the potential energy and dipole moment surfaces; and (iii) a competitive level of accuracy compared to the non-linear PTS model. Calculations of CH4 PTS, with polarizabilities fitted to 9000 training set points with the energies up to 14,000 cm-1 show an impressive level of accuracy of the linear PTS model obtained with ~1600 parameters: ~1% versus 0.3% RMSE for the non-linear vs. linear model on a test set of 1000 configurations.
Collapse
Affiliation(s)
- Oluwaseun Omodemi
- Department of Chemistry & Biochemistry, Kennesaw State University, Kennesaw, Georgia, USA
| | - Martina Kaledin
- Department of Chemistry & Biochemistry, Kennesaw State University, Kennesaw, Georgia, USA
| | - Alexey L Kaledin
- Cherry L. Emerson Center for Scientific Computation and Department of Chemistry, Emory University, Atlanta, Georgia, USA
| |
Collapse
|
2
|
Omodemi O, Sprouse S, Herbert D, Kaledin M, Kaledin AL. On the Cartesian Representation of the Molecular Polarizability Tensor Surface by Polynomial Fitting to Ab Initio Data. J Chem Theory Comput 2021; 18:37-45. [PMID: 34958587 DOI: 10.1021/acs.jctc.1c01015] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
We describe an approach to constructing an analytic Cartesian representation of the molecular dipole polarizability tensor surface in terms of polynomials in interatomic distances with a training set of ab initio data points obtained from a molecular dynamics (MD) simulation or by any other available means. The proposed formulation is based on a perturbation treatment of the unmodified point dipole polarizability model of Applequist [ J. Am. Chem. Soc. 1972, 94, 2952] and is shown here to be, by construction (i) free of short-range or other singularities or discontinuities, (ii) symmetric and translationally invariant, and (iii) nonreliant on a body-fixed coordinate system. Permutational invariance of like nuclei is demonstrated to be readily applicable, making this approach useful for highly fluxional and reactive systems. Derivation of the method is described in detail, adding brief didactic numerical examples of H2 and H2O and concluding with an MD simulation of the Raman spectrum of H5O2+ at 300 K with the polarizability tensor fitted to CCSD(T)/aug-cc-pVTZ data obtained using the HBB-4B potential [ J. Chem. Phys. 2005, 122, 044308].
Collapse
Affiliation(s)
- Oluwaseun Omodemi
- Department of Chemistry & Biochemistry, Kennesaw State University, 370 Paulding Avenue NW, Box # 1203, Kennesaw, Georgia 30144, United States
| | - Sarah Sprouse
- Department of Chemistry & Biochemistry, Kennesaw State University, 370 Paulding Avenue NW, Box # 1203, Kennesaw, Georgia 30144, United States
| | - Destyni Herbert
- Department of Chemistry & Biochemistry, Kennesaw State University, 370 Paulding Avenue NW, Box # 1203, Kennesaw, Georgia 30144, United States
| | - Martina Kaledin
- Department of Chemistry & Biochemistry, Kennesaw State University, 370 Paulding Avenue NW, Box # 1203, Kennesaw, Georgia 30144, United States
| | - Alexey L Kaledin
- Cherry L. Emerson Center for Scientific Computation, Emory University, 1515 Dickey Drive, Atlanta, Georgia 30322, United States
| |
Collapse
|
3
|
Le Breton G, Bonhomme O, Brevet PF, Benichou E, Loison C. First hyperpolarizability of water at the air-vapor interface: a QM/MM study questions standard experimental approximations. Phys Chem Chem Phys 2021; 23:24932-24941. [PMID: 34726679 DOI: 10.1039/d1cp02258j] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Surface Second-Harmonic Generation (S-SHG) experiments provide a unique approach to probe interfaces. One important issue for S-SHG is how to interpret the S-SHG intensities at the molecular level. Established frameworks commonly assume that each molecule emits light according to an average molecular hyperpolarizability tensor β(-2ω,ω,ω). However, for water molecules, this first hyperpolarizability is known to be extremely sensitive to their environment. We have investigated the molecular first hyperpolarizability of water molecules within the liquid-vapor interface, using a quantum description with explicit, inhomogeneous electrostatic embedding. The resulting average molecular first hyperpolarizability tensor depends on the distance relative to the interface, and it practically respects the Kleinman symmetry everywhere in the liquid. Within this numerical approach, based on the dipolar approximation, the water layer contributing to the Surface Second Harmonic Generation (S-SHG) intensity is less than a nanometer. The results reported here question standard interpretations based on a single, averaged hyperpolarizability for all molecules at the interface. Not only the molecular first hyperpolarizability tensor significantly depends on the distance relative to the interface, but it is also correlated to the molecular orientation. Such hyperpolarizability fluctuations may impact the S-SHG intensity emitted by an aqueous interface.
Collapse
Affiliation(s)
- Guillaume Le Breton
- University of Lyon, Université Claude Bernard Lyon 1, CNRS, Institut Lumière Matière, F-69622, Villeurbanne, France.
| | - Oriane Bonhomme
- University of Lyon, Université Claude Bernard Lyon 1, CNRS, Institut Lumière Matière, F-69622, Villeurbanne, France.
| | - Pierre-François Brevet
- University of Lyon, Université Claude Bernard Lyon 1, CNRS, Institut Lumière Matière, F-69622, Villeurbanne, France.
| | - Emmanuel Benichou
- University of Lyon, Université Claude Bernard Lyon 1, CNRS, Institut Lumière Matière, F-69622, Villeurbanne, France.
| | - Claire Loison
- University of Lyon, Université Claude Bernard Lyon 1, CNRS, Institut Lumière Matière, F-69622, Villeurbanne, France.
| |
Collapse
|
4
|
Kosar N, Ayub K, Mahmood T. Surface functionalization of twisted graphene C 32H 15 and C 104H 52 derivatives with alkalis and superalkalis for NLO response; a DFT study. J Mol Graph Model 2020; 102:107794. [PMID: 33212383 DOI: 10.1016/j.jmgm.2020.107794] [Citation(s) in RCA: 15] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/04/2020] [Revised: 10/18/2020] [Accepted: 10/30/2020] [Indexed: 12/14/2022]
Abstract
Herein, we present the detailed comparative study on geometric, electronic, optical and non-linear optical response of alkalis and superalkalis doped twisted graphene. The results illustrate that alkali metals and superalkalis interact with the central ring of the twisted graphene through non-covalent interactions which demonstrate the stability of the resultant complexes. NBO charges indicate the transfer of electrons from dopant (alkali metal atoms and superalkalis) towards twisted graphene sheet. Superalkalis doped twisted graphene complexes exhibit higher first hyperpolarizability values compared to alkali metals analogues. Among superalkalis doped complexes, K3O@C104H52 shows the highest βo value of 1.68 × 105 au. In frequency dependent first hyperpolarizability analysis, strong second harmonic generation (SHG) response of K3O@C32H15 complex is observed at both selected resonance frequency values (532 nm and1064 nm) whereas EOPE value of K3O@C32H15 complex shows higher induced response at 1064 nm wavelength. The static hyperpolarizability (βo) further increases under the influence of applied electric field. Among all complexes, Li3O@C32H15 graphene complex has the highest βo value (1.40 × 105 au) under applied electric field along x axis when sheet is in y-z plane. This analysis will be an important guideline for future studies on twisted graphene based NLO materials.
Collapse
Affiliation(s)
- Naveen Kosar
- Department of Chemistry, University of Management and Technology (UMT), C11, Johar Town Lahore, Pakistan
| | - Khurshid Ayub
- Department of Chemistry, COMSATS University Islamabad, Abbottabad Campus, Abbottabad, 22060, Pakistan
| | - Tariq Mahmood
- Department of Chemistry, COMSATS University Islamabad, Abbottabad Campus, Abbottabad, 22060, Pakistan.
| |
Collapse
|
5
|
Abstract
We review recent work on property decomposition techniques using quantum chemical methods and discuss some topical applications in terms of quantum mechanics-molecular mechanics calculations and the constructing of properties of large molecules and clusters. Starting out from the so-called LoProp decomposition scheme [Gagliardi et al., J. Chem. Phys., 2004, 121, 4994] for extracting atomic and inter-atomic contributions to molecular properties we show how this method can be generalized to localized frequency-dependent polarizabilities, to localized hyperpolarizabilities and to localized dispersion coefficients. Some applications of the generalized decomposition technique are reviewed - calculations of frequency-dependent polarizabilities, Rayleigh scattering of large clusters, and calculations of hyperpolarizabilities of proteins.
Collapse
Affiliation(s)
- Hans Ågren
- KTH Royal Institute of Technology, School of Engineering Sciences in Chemistry, Biotechnology and Health, Department of Theoretical Chemistry and Biology, SE-106 91 Stockholm, Sweden.
| | | | | |
Collapse
|
6
|
Harczuk I, Nagy B, Jensen F, Vahtras O, Ågren H. Local decomposition of imaginary polarizabilities and dispersion coefficients. Phys Chem Chem Phys 2017; 19:20241-20250. [DOI: 10.1039/c7cp02399e] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/19/2022]
Abstract
We present a new way to compute the two-body contribution to the dispersion energy using ab initio theory.
Collapse
Affiliation(s)
- Ignat Harczuk
- KTH Royal Institute of Technology
- School of Biotechnology
- Division of Theoretical Chemistry and Biology
- SE-106 91 Stockholm
- Sweden
| | - Balazs Nagy
- Department of Chemistry
- Aarhus University
- Langelandsgade 140
- 8000 Aarhus C
- Denmark
| | - Frank Jensen
- Department of Chemistry
- Aarhus University
- Langelandsgade 140
- 8000 Aarhus C
- Denmark
| | - Olav Vahtras
- KTH Royal Institute of Technology
- School of Biotechnology
- Division of Theoretical Chemistry and Biology
- SE-106 91 Stockholm
- Sweden
| | - Hans Ågren
- KTH Royal Institute of Technology
- School of Biotechnology
- Division of Theoretical Chemistry and Biology
- SE-106 91 Stockholm
- Sweden
| |
Collapse
|
7
|
Boleininger M, Guilbert AAY, Horsfield AP. Gaussian polarizable-ion tight binding. J Chem Phys 2016; 145:144103. [DOI: 10.1063/1.4964391] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Affiliation(s)
- Max Boleininger
- Department of Physics and Thomas Young Centre, Imperial College London, South Kensington Campus, London SW7 2AZ, United Kingdom
| | - Anne AY Guilbert
- Department of Physics and Thomas Young Centre, Imperial College London, South Kensington Campus, London SW7 2AZ, United Kingdom
| | - Andrew P. Horsfield
- Department of Materials and Thomas Young Centre, Imperial College London, South Kensington Campus, London SW7 2AZ, United Kingdom
| |
Collapse
|
8
|
Harczuk I, Vahtras O, Ågren H. First Hyperpolarizability of Collagen Using the Point Dipole Approximation. J Phys Chem Lett 2016; 7:2132-2138. [PMID: 27203480 DOI: 10.1021/acs.jpclett.6b00721] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/05/2023]
Abstract
The application of localized hyperpolarizabilities to predict a total protein hyperpolarizability is presented for the first time, using rat-tail collagen as a demonstration example. We employ a model comprising the quadratic Applequist point-dipole approach, the so-called LoProp transformation, and a procedure with molecular fractionation using conjugate caps to determine the atomic and bond contributions to the net β tensor of the collagen [(PPG)10]3 triple-helix. By using Tholes exponential damping modification to the dyadic tensor in the Applequist equations, a correct qualitative agreement with experiment is found. The intensity of the βHRS signal and the depolarization ratios are best reproduced by decomposing the LoProp properties into the atomic positions and using Tholes exponential damping with the original damping parameter. Some ramifications of the model for general protein property optimization are briefly discussed.
Collapse
Affiliation(s)
- Ignat Harczuk
- School of Biotechnology, Division of Theoretical Chemistry and Biology, KTH Royal Institute of Technology , SE-106 91 Stockholm, Sweden
| | - Olav Vahtras
- School of Biotechnology, Division of Theoretical Chemistry and Biology, KTH Royal Institute of Technology , SE-106 91 Stockholm, Sweden
| | - Hans Ågren
- School of Biotechnology, Division of Theoretical Chemistry and Biology, KTH Royal Institute of Technology , SE-106 91 Stockholm, Sweden
| |
Collapse
|
9
|
Harczuk I, Vahtras O, Ågren H. Modeling Rayleigh Scattering of Aerosol Particles. J Phys Chem B 2016; 120:4296-301. [DOI: 10.1021/acs.jpcb.6b02278] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Ignat Harczuk
- School of Biotechnology,
Division of Theoretical Chemistry and Biology, KTH Royal Institute of Technology, SE-106 91 Stockholm, Sweden
| | - Olav Vahtras
- School of Biotechnology,
Division of Theoretical Chemistry and Biology, KTH Royal Institute of Technology, SE-106 91 Stockholm, Sweden
| | - Hans Ågren
- School of Biotechnology,
Division of Theoretical Chemistry and Biology, KTH Royal Institute of Technology, SE-106 91 Stockholm, Sweden
| |
Collapse
|
10
|
Cardenuto MH, Castet F, Champagne B. Investigating the first hyperpolarizability of liquid carbon tetrachloride. RSC Adv 2016. [DOI: 10.1039/c6ra23187j] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022] Open
Abstract
Sequential QMMM calculations have been carried out to investigate the first hyperpolarizability of liquid CCl4.
Collapse
Affiliation(s)
- Marcelo Hidalgo Cardenuto
- Laboratoire de Chimie Théorique
- Unité de Chimie Physique Théorique et Structurale
- University of Namur
- B-5000 Namur
- Belgium
| | - Frédéric Castet
- Université de Bordeaux
- Institut des Sciences Moléculaires (ISM)
- UMR 5255 CNRS
- F-33405 Talence Cedex
- France
| | - Benoît Champagne
- Laboratoire de Chimie Théorique
- Unité de Chimie Physique Théorique et Structurale
- University of Namur
- B-5000 Namur
- Belgium
| |
Collapse
|
11
|
List NH, Olsen JMH, Kongsted J. Excited states in large molecular systems through polarizable embedding. Phys Chem Chem Phys 2016; 18:20234-50. [DOI: 10.1039/c6cp03834d] [Citation(s) in RCA: 60] [Impact Index Per Article: 7.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Using the polarizable embedding model enables rational design of light-sensitive functional biological materials.
Collapse
Affiliation(s)
- Nanna Holmgaard List
- Department of Physics, Chemistry and Pharmacy
- University of Southern Denmark
- 5230 Odense M
- Denmark
| | | | - Jacob Kongsted
- Department of Physics, Chemistry and Pharmacy
- University of Southern Denmark
- 5230 Odense M
- Denmark
| |
Collapse
|