Hybrid quantum-classical polarizability model for single molecule biosensing

Optical whispering gallery mode biosensors are able to detect single molecules through effects of their polarizability. We address the factors that affect the polarizability of amino acids, which are the building blocks of life, via electronic structure theory. Amino acids are detected in aqueous environments, where their polarizability is different compared to the gasphase due to solvent effects. Solvent effects include structural changes, protonation and the local field enhancement through the solvent (water). We analyse the impact of these effects and find that all contribute to an increased effective polarizability in the solvent. We also address the excess polarizability relative to the displaced water cavity and develop a hybrid quantum-classical model that is in good agreement with self-consistent calculations. We apply our model to calculate the excess polarizability of 20 proteinogenic amino acids and determine the minimum resolution required to distinguish the different molecules and their ionised conformers based on their polarizability.

A Curious Case of Two-Photon Absorption in n-Helicene and n-Phenylene, n=6-10: Why n=7 is Different?

n-Helicenes and n-Phenylenes are interesting examples of twisted molecules, where although the atoms are connected through conjugated π ${\pi }$ -bonds, the π ${\pi }$ -conjugation is largely hindered by the twisted nature of the bonds. Such structures provide a unique opportunity to study the effect of twisted π ${\pi }$ -system on non-linear optical properties. In this work, we studied the two-photon absorption in donor-acceptor substituted n-helicenes and n-phenylenes employing the state-of-the-art RI-CC2 method and reported a unique feature we observed in n=7 systems. We found that both 7-helicene and 7-phenylene systems exhibit largest two-photon absorption than other members in their respective classes. Furthermore, using generalized few-state model, we provided a detailed microscopic mechanism of this unique observation involving participation of different transition dipole moment vectors and their relative orientations.

Chromophore Planarity, -BH Bridge Effect, and Two-Photon Activity: Bi- and Ter-Phenyl Derivatives as a Case Study

In this work, we have employed electronic structure theories to explore the effect of the planarity of the chromophore on the two-photon absorption properties of bi- and ter-phenyl systems. To that end, we have considered 11 bi- and 7 ter-phenyl-based chromophores presenting a donor-π-acceptor architecture. In some cases, the planarity has been enforced by bridging the rings at ortho-positions by -CH2 and/or -BH, -O, -S, and -NH moieties. The results presented herein demonstrate that in bi- and ter-phenyl systems, the planarity achieved via a -CH2 bridge increases the 2PA activity. However, the introduction of a bridge with the -BH moiety perturbs the electronic structure to a large extent, thus diminishing the two-photon transition strength to the lowest electronic excited state. As far as two-photon absorption activity is concerned, this work hints toward avoiding -BH bridge(s) to enforce planarity in bi- and ter-phenyl systems; however, one may use -CH2 bridge(s) to achieve the enhancement of the property in question. All of these conclusions have been supported by in-depth analyses based on generalized few-state models.

First Hyperpolarizabilities of Intramolecular Charge-Transfer Architectures Based on Acenaphthene Derivatives in Gas, Solution, and Solid States

Constructing charge transfer (CT) systems and packing arrangement are common and effective methods to control the efficiency of nonlinear optical (NLO) materials. Apart from the traditional through-bond CT (TBCT) systems, through-space CT (TSCT) also leads to distinctive optical and electronic properties. Meanwhile, corresponding theoretical investigations of the aggregation effect are highly desired. In this work, some TSCT and model compounds incorporating acenaphthene as a scaffold and triphenylamine (TPA) as the donor are theoretically performed to systematically reveal the effect of both solvent and solid environments on their static first hyperpolarizabilities (β_tot) by using the polarizable continuum model (PCM) and the combined quantum mechanics and molecular mechanics (QM/MM) method. Results indicate that the dichloromethane solvent effect within the PCM approach causes an almost 2 times increase of the β_tot values. Besides, the different packing modes and intermolecular interactions have remarkable influence on the second-order NLO properties. For the case of TPA-ace-CN in the crystal state, the parallel arrangement will lead to large NLO responses (4.9 × 10^-30 esu) compared to the correspondingly isolated molecule (3.4 × 10^-30 esu). However, for the TPA-ace-TRZ compound with the TSCT architecture, selection of the molecular arrangement may make the aggregate ineffective due to the offset of the through-space dipole and charge transfer between D-A groups, which lead to the β_tot values decreasing from 15.2 × 10^-30 to 7.7 × 10^-30 esu. We believe that our calculation will serve as a guide for the exploration of more efficient NLO materials wherein the molecules are oriented in their most favorable arrangements.

The orbital picture of the first dipole hyperpolarizability from many-body response theory

We present an approach for obtaining a molecular orbital picture of the first dipole hyperpolarizability (β) from correlated many-body electronic structure methods. Ab initio calculations of β rely on quadratic response theory, which recasts the sum-over-all-states expression of β into a closed-form expression by calculating a handful of first- and second-order response states; for resonantly enhanced β, damped response theory is used. These response states are then used to construct second-order response reduced one-particle density matrices (1PDMs), which, upon visualization in terms of natural orbitals (NOs), facilitate a rigorous and black-box mapping of the underlying electronic structure with β. We explain the interpretation of different components of the response 1PDMs and the corresponding NOs within both the undamped and damped response theory framework. We illustrate the utility of this new tool by deconstructing β for cis-difluoroethene, para-nitroaniline, and hemibonded OH^· + H₂O complex, computed within the framework of coupled-cluster singles and doubles response theory, in terms of the underlying response 1PDMs and NOs for a range of frequencies.

Benchmarking Correlated Methods for Static and Dynamic Polarizabilities: The T145 Data Set Evaluated with RPA, RPA(D), HRPA, HRPA(D), SOPPA, SOPPA(CC2), SOPPA(CCSD), CC2, and CCSD

Due to the importance of predicting static and dynamic polarizabilities, the performance of various correlated linear response methods including random phase approximation (RPA), RPA(D), higher-order random phase approximation (HRPA), HRPA(D), second-order polarization propagator approximation (SOPPA), SOPPA(CC2), SOPPA(CCSD), CC2, and CCSD has been evaluated against CCSD(T) (static case) and CCSD (dynamic cases) for the T145 set of 145 organic molecules. The benchmark reveals that the HRPA(D) method has the best performance for both static and dynamic polarizabilities apart from CCSD. RPA(D) ranks second for the dynamic cases and third for the static case. Using coupled-cluster amplitudes in SOPPA(CCSD) and SOPPA(CC2), the SOPPA results are significantly improved. The HRPA method has the largest deviations from the reference values for both cases. In general, according to the performance and computational cost of the methods, the HRPA(D) and RPA(D) methods are proposed for calculations of static and dynamic polarizabilities of this and similar sets of molecules.

A generalized few-state model for the first hyperpolarizability

The properties of molecules depend on their chemical structure, and thus, structure-property relations help design molecules with desired properties. Few-state models are often used to interpret experimental observations of non-linear optical properties. Not only the magnitude but also the relative orientation of the transition dipole moment vectors is needed for few-state models of the non-linear optical properties. The effect of the relative orientation of the transition dipole moment vectors is called dipole alignment, and this effect has previously been studied for multiphoton absorption properties. However, so far, no such studies are reported for the first hyperpolarizability. Here, we present a generalized few-state model for the static and dynamic first hyperpolarizability β, accounting for the effect of dipole alignment. The formulas derived in this work are general in the sense that they can be used for any few-state model, i.e., a two-state model, a three-state model, or, in general, an n-state model. Based on the formulas, we formulate minimization and maximization criteria for the alignment of transition dipole moment vectors. We demonstrate the importance of dipole alignment by applying the formulas to the static first hyperpolarizability of ortho-, meta-, and para-nitroaniline. The formulas and the analysis provide new ways to understand the structure-property relationship for β and can hence be used to fine-tune the magnitude of β in a molecule.

Tuning of hyperpolarizability, and one- and two-photon absorption of donor-acceptor and donor-acceptor-acceptor-type intramolecular charge transfer-based sensors

The present work aims to study the effect of solvent as well as arrangement of donor-acceptor groups on linear and non-linear optical (NLO) response properties of two experimentally studied intramolecular charge-transfer (ICT)-based fluorescent sensors. One of them (molecule 1) is a donor-acceptor (D-A) system with hemicyanine and dimethylanilino as electron withdrawing and donating groups, respectively, while the other one (molecule 3) is molecule 1 fused with a boron-dipyrromethene (BODIPY) moiety. BODIPY acts as the electron acceptor group of molecule 2 that as well consists of dimethylanilino as the electron donor. Density functional theory (DFT) as well as time-dependent DFT has been employed to optimize the geometry of the molecules, followed by computation of dipole moment (μ), static first hyperpolarizability (β_total), and one- and two-photon absorption (TPA) strengths. The results reveal that dipole moment as well as total static first hyperpolarizability (β_total) of the studied molecules is dominated by the respective components in the direction of charge transfer. The ratio of vector component of first hyperpolarizability (β_vec) to β_total also supports the unidirectional charge transfer in the studied systems. In molecule 3, which is a donor-acceptor-acceptor (D-A-A)-type system, the BODIPY moiety is found to play a major role in controlling the NLO response over the other acceptor group. Solvents are also found to play an important role in controlling the linear as well as NLO response of the studied systems. A significant increase in the first hyperpolarizability as well as TPA cross-section of the studied molecules is predicted due to an increase in the dielectric constant of the medium. The results presented are expected to provide a clue in tuning the NLO response of many ICT-based chromophores, especially those with D-A-A arrangements.