Recent Advances in Wave Function-Based Methods of Molecular-Property Calculations

SCC-DFTB calculation of the static first hyperpolarizability: from gas phase molecules to functionalized surfaces

The performance of the self-consistent charge density functional tight binding (SCC-DFTB) method for calculating the first hyperpolarizability of π-conjugated compounds has been assessed with respect to results obtained with high-level ab initio methods and density functional theory (DFT). The SCC-DFTB method performs similarly or better than DFT with the PBE XC functional. Thus, if for small π-conjugated linkers SCC-DFTB can reproduce trends, for longer chains the first hyperpolarizabilities are overestimated. In the case of push-pull thiophenes, the β values are strongly overestimated, as it is also the case with the B3LYP and PBE XC functionals. On the other hand, the SCC-DFTB method closely reproduces the evolution of β in p-disubstituted benzenes as a function of the donor and acceptor groups, as estimated at the MP2 level. The reliability of SCC-DFTB to determine the bond length alternation and the dihedral angles between the aromatic rings has also been tackled, demonstrating that both are underestimated. Overall, the SCC-DFTB calculations are of the same quality as those performed with the conventional PBE XC functional on which the method was parameterized but the SCC-DFTB calculations are computationally very little demanding, and it can therefore be adopted for very large systems for screening nonlinear optical materials as well as for assessing structure-property relationships. This is illustrated with an application on the first hyperpolarizability of an indolino-oxazolidine molecular switch grafted on a SiO2 surface. This has enabled to pinpoint (i) the effect of the surface on the donor/acceptor character of the linking substituent, (ii) the impact of molecular orientation, (iii) the role of a spacer between the π-conjugated switch and the surface, (iv) the global effect of the surface on the β contrast, and also (v) the fact that the molecular switches can maintain this contrast when adsorbed.

The calculation of 29Si NMR chemical shifts of tetracoordinated silicon compounds in the gas phase and in solution

²⁹Si NMR chemical shifts have been calculated for a variety of silicon-based species relevant in organic synthesis.

On the origin of the very strong two-photon activity of squaraine dyes – a standard/damped response theory study

In the present work, we report the mechanism of a very large increase in the two-photon (TP) activity of squaraine based molecules upon changing the substituents.

Design and characterization of molecular nonlinear optical switches

Nanoscale structures, including molecules, supramolecules, polymers, functionalized surfaces, and crystalline/amorphous solids, can commute between two or more forms, displaying contrasts in their nonlinear optical (NLO) properties. Because of this property, they have high potential for applications in data storage, signal processing, and sensing. As potential candidates for integration into responsive materials, scientists have been intensely studying organic and organometallic molecules with switchable first hyperpolarizability over the past two decades. As a result of this, researchers have been able to synthesize and characterize several families of molecular NLO switches that differ by the stimulus used to trigger the commutation. These stimuli can include light irradiation, pH variation, redox reaction, and ion recognition, among others. The design of multistate (including several switchable units) and multifunctional (triggered with different stimuli) systems has also motivated a large amount of work, aiming at the improvement of the storage capacity of optical memories or the diversification of the addressability of the devices. In complement to the synthesis of the compounds and the characterization of their NLO responses by means of hyper-Rayleigh scattering, quantum chemical calculations play a key role in the design of molecular switches with high first hyperpolarizability contrasts. Through the latter, we can gain a fundamental understanding of the various factors governing the efficiency of the switches. These are not easily accessible experimentally, and include donor/acceptor contributions, frequency dispersion, and solvent effects. In this Account, we illustrate the similarities of the experimental and theoretical tools to design and characterize highly efficient NLO switches but also the difficulties in comparing them. After providing a critical overview of the different theoretical approaches used for evaluating the first hyperpolarizabilities, we report two case studies in which theoretical simulations have provided guidelines to design NLO switches with improved efficiencies. The first example presents the joint theoretical/experimental characterization of a new family of multi-addressable NLO switches based on benzazolo-oxazolidine derivatives. The second focuses on the photoinduced commutation in merocyanine-spiropyran systems, where the significant NLO contrast could be exploited for metal cation identification in a new generation of multiusage sensing devices. Finally, we illustrate the impact of environment on the NLO switching properties, with examples based on the keto-enol equilibrium in anil derivatives. Through these representative examples, we demonstrate that the rational design of molecular NLO switches, which combines experimental and theoretical approaches, has reached maturity. Future challenges consist in extending the investigated objects to supramolecular architectures involving several NLO-responsive units, in order to exploit their cooperative effects for enhancing the NLO responses and contrasts.

Halogen effect on structure and 13C NMR chemical shift of 3,6-disubstituted-N-alkyl carbazoles

Structures of selected 3,6-dihalogeno-N-alkyl carbazole derivatives were calculated at the B3LYP/6-311++G(3df,2pd) level of theory, and their (13) C nuclear magnetic resonance (NMR) isotropic shieldings were predicted using density functional theory (DFT). The model compounds contained 9H, N-methyl and N-ethyl derivatives. The relativistic effect of Br and I atoms on nuclear shieldings was modeled using the spin-orbit zeroth-order regular approximation (ZORA) method. Significant heavy atom shielding effects for the carbon atom directly bonded with Br and I were observed (~-10 and ~-30 ppm while the other carbon shifts were practically unaffected). The decreasing electronegativity of the halogen substituent (F, Cl, Br, and I) was reflected in both nonrelativistic and relativistic NMR results as decreased values of chemical shifts of carbon atoms attached to halogen (C3 and C6) leading to a strong sensitivity to halogen atom type at 3 and 6 positions of the carbazole ring. The predicted NMR data correctly reproduce the available experimental data for unsubstituted N-alkylcarbazoles.

The Dalton quantum chemistry program system

Dalton is a powerful general-purpose program system for the study of molecular electronic structure at the Hartree-Fock, Kohn-Sham, multiconfigurational self-consistent-field, Møller-Plesset, configuration-interaction, and coupled-cluster levels of theory. Apart from the total energy, a wide variety of molecular properties may be calculated using these electronic-structure models. Molecular gradients and Hessians are available for geometry optimizations, molecular dynamics, and vibrational studies, whereas magnetic resonance and optical activity can be studied in a gauge-origin-invariant manner. Frequency-dependent molecular properties can be calculated using linear, quadratic, and cubic response theory. A large number of singlet and triplet perturbation operators are available for the study of one-, two-, and three-photon processes. Environmental effects may be included using various dielectric-medium and quantum-mechanics/molecular-mechanics models. Large molecules may be studied using linear-scaling and massively parallel algorithms. Dalton is distributed at no cost from http://www.daltonprogram.org for a number of UNIX platforms.

One-bond 29Si-1H spin-spin coupling constants in the series of halosilanes: benchmark SOPPA and DFT calculations, relativistic effects, and vibrational corrections

A number of most representative second order polarization propagator approach (SOPPA) based wavefunction methods, SOPPA, SOPPA(CC2) and SOPPA(CCSD), and density functional theory (DFT) based methods, B3LYP, PBE0, KT2, and KT3, have been benchmarked in the calculation of the one-bond (29)Si-(1)H spin-spin coupling constants in the series of halosilanes SiH(n)X(4-n) (X = F, Cl, Br, I), both at the non-relativistic and full four-parameter Dirac's relativistic levels taking into account vibrational corrections. At the non-relativistic level, the wavefunction methods showed much better results as compared with those of DFT. At the DFT level, out of four tested functionals, the Perdew, Burke, and Ernzerhof's PBE0 showed best performance. Taking into account, relativistic effects and vibrational corrections noticeably improves wavefunction methods results, but generally worsens DFT results.

Estimation of isotropic nuclear magnetic shieldings in the CCSD(T) and MP2 complete basis set limit using affordable correlation calculations

A linear correlation between isotropic nuclear magnetic shielding constants for seven model molecules (CH2 O, H2 O, HF, F2 , HCN, SiH4 and H2 S) calculated with 37 methods (34 density functionals, RHF, MP2 and CCSD(T)), with affordable pcS-2 basis set and corresponding complete basis set results, estimated from calculations with the family of polarization-consistent pcS-n basis sets is reported. This dependence was also supported by inspection of profiles of deviation between CBS estimated nuclear shieldings and shieldings obtained with the significantly smaller basis sets pcS-2 and aug-cc-pVTZ-J for the selected set of 37 calculation methods. It was possible to formulate a practical approach of estimating the values of isotropic nuclear magnetic shielding constants at the CCSD(T)/CBS and MP2/CBS levels from affordable CCSD(T)/pcS-2, MP2/pcS-2 and DFT/CBS calculations with pcS-n basis sets. The proposed method leads to a fairly accurate estimation of nuclear magnetic shieldings and considerable saving of computational efforts.

Multiconfigurational Self-Consistent Field Calculations of the Magnetically Induced Current Density Using Gauge-Including Atomic Orbitals

Nondynamical electron correlation based on a genuine multiconfigurational theory is of considerable importance for a balanced ab initio calculation of aromatic and antiaromatic molecules either with open-shell character or quasi-degeneracy in the electronic states. Among the various aromaticity indices, the calculation of magnetically induced ring current densities (MICD) has emerged as a strong contender, providing both a qualitative and a quantitative description of the effect. We report here the first implementation of MICD at the multiconfigurational self-consistent field (MCSCF) level of theory together with example calculations. This extension makes the method applicable to systems that cannot be appropriately handled with earlier implementations based on a single-reference starting function. We present the formulation of the MCSCF MICD theory along with applications to a prototypical antiaromatic (cyclobutadiene) and an aromatic (benzyne) system, both systems that require a multiconfigurational description. We compare the MCSCF results to those obtained using Hartree-Fock and Kohn-Sham density functional theory and discuss the effects of static correlation on the aromaticity.

The Absolute Shielding Constants of Heavy Nuclei: Resolving the Enigma of the (119)Sn Absolute Shielding

We demonstrate that the apparent disagreement between experimental determinations and four-component relativistic calculations of the absolute shielding constants of heavy nuclei is due to the breakdown of the commonly assumed relation between the electronic contribution to the nuclear spin-rotation constants and the paramagnetic contribution to the NMR shielding constants. We demonstrate that this breakdown has significant consequences for the absolute shielding constant of (119)Sn, leading to errors of about 1000 ppm. As a consequence, we expect that many absolute shielding constants of heavy nuclei will be in need of revision.