Toward Accurate Two-Photon Absorption Spectrum Simulations: Exploring the Landscape beyond the Generalized Gradient Approximation.

In this Letter, we present a pioneering analysis of the density functional approximations (DFAs) beyond the generalized gradient approximation (GGA) for predicting two-photon absorption (2PA) strengths of a set of push-pull π-conjugated molecules. In more detail, we have employed a variety of meta-generalized gradient approximation (meta-GGA) functionals, including SCAN, MN15, and M06-2X, to assess their accuracy in describing the 2PA properties of a chosen set of 48 organic molecules. Analytic quadratic response theory is employed for these functionals, and their performance is compared against the previously studied DFAs and reference data obtained at the coupled-cluster CC2 level combined with the resolution-of-identity approximation (RI-CC2). A detailed analysis of the meta-GGA functional performance is provided, demonstrating that they improve upon their predecessors in capturing the key electronic features of the π-conjugated two-photon absorbers. In particular, the Minnesota functional MN15 shows very promising results as it delivers pleasingly accurate chemical rankings for two-photon transition strengths and excited-state dipole moments.

Radiation-Driven Polymerisation of Methacrylic Acid in Aqueous Solution: A Chemical Events Monte Carlo Study.

This study employed a coarse-grained Monte Carlo (MC) simulation to investigate the radiation-induced polymerisation of methacrylic acid (MAA) in an aqueous solution. This method provides an alternative to traditional kinetic models, enabling a detailed examination of the micro-structure and growth patterns of MAA polymers, which are often not captured in other approaches. In this work, we generated multiple clones of a simulation box, each containing a specific chemical composition. In these simulations, every coarse-grained (CG) bead represents an entire monomer. The growth function, defined by the chemical behaviour of interacting substances, was determined through repeated random sampling. This approach allowed us to simulate the complex process of radiation-induced polymerisation, enhancing our understanding of the formation of poly(methacrylic acid) hydrogels at a microscopic level; while Monte Carlo simulations have been applied in various contexts of polymerisation, this study's specific approach to modelling the radiation-induced polymerisation of MAA in an aqueous environment, utilising the data obtained by quantum chemistry modelling, with an emphasis on micro-structural growth, has not been extensively explored in existing studies. This understanding is important for advancing the synthesis of these hydrogels, which have potential applications in diverse fields such as materials science and medicine.

A polarizable coarse-grained model for metal, metal oxide and composite metal/metal oxide nanoparticles: development and implementation.

We present a polarizable coarse-grained model for metal, metal oxide, and composite metal/metal oxide nanoparticles with well-defined crystalline surfaces. The developed model uses a low-resolution polarizable "surface beads" representation of the nanoparticle's geometry and pairwise cross nanoparticle potential consisting of van der Waals and electrostatic interaction terms. The electrostatic interaction term of the cross nanoparticle potential incorporates a crucial physical aspect of electrostatic interaction into the metal and metal oxide systems, such as induced surface charges, making it possible to explore the nanoparticles' behavior in complex environments as well as investigate the interplay between electrostatic and van der Waals interactions in nanoparticle systems. The iterative stability, computational scaling, and performance of the presented model was tested on selected systems of gold, titanium dioxide, and composite gold/titanium dioxide nanoparticle systems. The model exhibits robust iterative stability and is able to converge the charge equilibration equation for fluctuating induced charges and dipoles within 10-60 "tug-tow" iterations in challenging situations, like crowded nanoparticle systems or nanoparticle systems in extreme external electric fields. The computation scaling of the presented model is semi-linear with respect to the number of nanoparticles in the system. It slightly varies depending on the size distribution of nanoparticles in a specific nanoparticle system. The computation cost of the model is significantly lower than that of conventional atomistic polarizable force field models and enables the treatment of large nanoparticle systems that are beyond the reach of currently existing atomistic force field models.

A polarizable coarse-grained model for metal, metal oxide and composite metal/metal oxide nanoparticles and its applications.

We present a selected set of exemplifying applications of the novel polarizable coarse-grained model [see the first part] to various outstanding problems in the physics and chemistry of nanoparticles: electrostatic potential around silver and gold nanoparticles; spontaneous and external electric field-driven self-organization of gold and silver nanoparticle systems; and physisorption of carbon dioxide on titanium dioxide nanoparticles decorated with a gold catalyst. In the first application, the developed model has shown capabilities of predicting long-range potential with accuracy comparable to the tight-binding density functional theory methods. Furthermore, the electrostatic potential analysis in hot spot regions allowed us to identify twin defect lines in a silver nanostar as a promising candidate for an enhancer in surface-enhanced Raman spectroscopy. In the second application, the developed model has facilitated the elucidation of the microscopic mechanisms responsible for the self-organization of gold and silver nanoparticles. Analysis of Monte-Carlo simulations established that the self-organization process is driven by van der Waals interactions in the absence of an external electric field, and that it becomes gradually driven by electrostatic interactions in the presence of an external electric field with increasing strength of the external electric field. In the third application, the developed model combined with Monte-Carlo simulations has identified the dominant mechanism responsible for carbon dioxide transfer to the catalytic sites. Analysis of the obtained results indicates that surface diffusion is the dominant mechanism for the transport of carbon dioxide to the catalytic sites, and only in exceptional situations, direct physisorption becomes a competitive mechanism with the surface diffusion mechanism. These successful applications of the developed model indicate its wide range of applicability to various problems in the chemistry and physics of nanoparticles.

Antibacterial Activity of Silver and Gold Particles Formed on Titania Thin Films.

Metal-based nanoparticles with antimicrobial activity are gaining a lot of attention in recent years due to the increased antibiotics resistance. The development and the pathogenesis of oral diseases are usually associated with the formation of bacteria biofilms on the surfaces; therefore, it is crucial to investigate the materials and their properties that would reduce bacterial attachment and biofilm formation. This work provides a systematic investigation of the physical-chemical properties and the antibacterial activity of TiO2 thin films decorated by Ag and Au nanoparticles (NP) against Veillonella parvula and Neisseria sicca species associated with oral diseases. TiO2 thin films were formed using reactive magnetron sputtering by obtaining as-deposited amorphous and crystalline TiO2 thin films after annealing. Au and Ag NP were formed using a two-step process: magnetron sputtering of thin metal films and solid-state dewetting. The surface properties and crystallographic nature of TiO2/NP structures were investigated by SEM, XPS, XRD, and optical microscopy. It was found that the higher thickness of Au and Ag thin films results in the formation of the enlarged NPs and increased distance between them, influencing the antibacterial activity of the formed structures. TiO2 surface with AgNP exhibited higher antibacterial efficiency than Au nanostructured titania surfaces and effectively reduced the concentration of the bacteria. The process of the observation and identification of the presence of bacteria using the deep learning technique was realized.

Electronic g-Tensor Calculations for Dangling Bonds in Nanodiamonds.

The electronic g-tensor calculations are performed for dangling bonds (DBs) introduced into nanodiamonds (NDs) with four different functional groups on their surfaces. For hydrogenated and fluorinated NDs, it is found that g-shifts of the latter vary in a much wider range, and the same is also true for the total energy differences between the highest and the lowest energy DBs. In addition, it is shown that the shape of NDs significantly impacts the energetics and g-shifts of DBs, whereas the influence of the size is much less pronounced, as is the influence of the presence of one DB in the vicinity of the other, resulting in no substantial change on their magnetic behavior. For hydroxylated and aminated NDs, it is demonstrated that the variation range of g-shifts is larger for the former, whereas the opposite is seen regarding the total energy differences. On the whole, some of the positions of DBs can be energetically very costly in these NDs; besides, the lowest energy DBs are irregular, that is, formed by OH- and NH2-bonded C atoms, contrasting with hydrogenated and fluorinated NDs, for which irregular DBs are the most energetically unfavorable.

Size-dependent polarizabilities and van der Waals dispersion coefficients of fullerenes from large-scale complex polarization propagator calculations.

While the anomalous non-additive size-dependencies of static dipole polarizabilities and van der Waals C6 dispersion coefficients of carbon fullerenes are well established, the widespread reported scalings for the latter (ranging from N2.2 to N2.8) call for a comprehensive first-principles investigation. With a highly efficient implementation of the linear complex polarization propagator, we have performed Hartree-Fock and Kohn-Sham density functional theory calculations of the frequency-dependent polarizabilities for fullerenes consisting of up to 540 carbon atoms. Our results for the static polarizabilities and C6 coefficients show scalings of N1.2 and N2.2, respectively, thereby deviating significantly from the previously reported values obtained with the use of semi-classical/empirical methods. Arguably, our reported values are the most accurate to date as they represent the first ab initio or first-principles treatment of fullerenes up to a convincing system size.

Efficient implementation of isotropic cubic response functions for two-photon absorption cross sections within the self-consistent field approximation.

Within the self-consistent field approximation, computationally tractable expressions for the isotropic second-order hyperpolarizability have been derived and implemented for the calculation of two-photon absorption cross sections. The novel tensor average formulation presented in this work allows for the evaluation of isotropic damped cubic response functions using only ∼3.3% (one-photon off-resonance regions) and ∼10% (one-photon resonance regions) of the number of auxiliary Fock matrices required when explicitly calculating all the needed individual tensor components. Numerical examples of the two-photon absorption cross section in the one-photon off-resonance and resonance regions are provided for alanine-tryptophan and 2,5-dibromo-1,4-bis(2-(4-diphenylaminophenyl)vinyl)-benzene. Furthermore, a benchmark set of 22 additional small- and medium-sized organic molecules is considered. In all these calculations, a quantitative assessment is made of the reduced and approximate forms of the cubic response function in the one-photon off-resonance regions and results demonstrate a relative error of less than ∼5% when using the reduced expression as compared to the full form of the isotropic cubic response function.

Dalton Project: A Python platform for molecular- and electronic-structure simulations of complex systems.

The Dalton Project provides a uniform platform access to the underlying full-fledged quantum chemistry codes Dalton and LSDalton as well as the PyFraME package for automatized fragmentation and parameterization of complex molecular environments. The platform is written in Python and defines a means for library communication and interaction. Intermediate data such as integrals are exposed to the platform and made accessible to the user in the form of NumPy arrays, and the resulting data are extracted, analyzed, and visualized. Complex computational protocols that may, for instance, arise due to a need for environment fragmentation and configuration-space sampling of biochemical systems are readily assisted by the platform. The platform is designed to host additional software libraries and will serve as a hub for future modular software development efforts in the distributed Dalton community.

Plasmonic nano-shells: atomistic discrete interaction versus classic electrodynamics models.

Using the extended discrete interaction model and Mie theory, we investigate the tunability of the optical polarizability of small metallic nano-shells. We show that the spectral positions of symmetric and antisymmetric dipolar plasmon resonances vary with the ratio of particle radius to hole radius in a manner similar to one predicted for uniform metallic nano-shells using a semiclassical approach of two coupled harmonic oscillators. We show that, according to the extended discrete interaction model, the dipolar plasmon resonances are also present for nano-shells in the 2-13 nm size region and show the same functional dependence seen for larger nano-shells. Using previously fitted data from experiment, we can predict the size-dependence of the plasma frequency for nano-shells in the 1-15 nm size region. We find that Mie theory, which utilizes the electron mean free path correction for the permittivity, is not able to reproduce the same functional form of the dipolar modes for the nano-shells of the same sizes.

Optomagnetic Effect Induced by Magnetized Nanocavity Plasmon.

We propose a new type of optomagnetic effect induced by a highly confined plasmonic field in a nanocavity. It is shown that a very large dynamic magnetic field can be generated as the result of the inhomogeneity of nanocavity plasmons, which can directly activate spin-forbidden transitions in molecules. The dynamic optomagnetic effects on optical transitions between states of different spin multiplicities are illustrated by first-principles calculations for C60. Remarkably, the intensity of spin forbidden singlet-to-triplet transitions can even be stronger than that of singlet-to-singlet transitions when the spatial distribution of plasmon is comparable with the molecular size. This approach not only offers a powerful optomagnetic means to rationally fabricate molecular excited states with different multiplicities but also provides a groundbreaking concept of the light-matter interaction that could lead to the observation of new physical phenomena and the development of new techniques.

Mechanistic Studies on NaHCO3 Hydrogenation and HCOOH Dehydrogenation Reactions Catalysed by a FeII Linear Tetraphosphine Complex.

We present a theoretical extension of the previously published bicarbonate hydrogenation to formate and formic acid dehydrogenation catalysed by FeII complexes bearing the linear tetraphosphine ligand tetraphos-1. The hydrogenation reaction was found to proceed at the singlet surface with two competing pathways: A) H2 association to the Fe-H species followed by deprotonation to give a Fe(H)2 intermediate, which then reacts with CO2 to give formate. B) CO2 insertion into the Fe-H bond, followed by H2 association and subsequent deprotonation. B was found to be slightly preferred with an activation energy of 22.8 kcal mol-1 , compared to 25.3 for A. Further we have reassigned the Fe-H complex, as a Fe(H)(H2 ), which undergoes extremely rapid hydrogen exchange.

Formation of N-oxide in the third oxidation of [RuII(tpy)(L)(OH2)]2.

The widely studied water oxidation catalyst [RuII(tpy)(L)(OH2)]2+ (tpy = 2,2':6',2''-terpyridine, L = bpy = 2,2'-bipyridine or L = bpm = 2,2'-bipyrimidine) is still under scrutiny. Here we present a new suggestion for one of the key steps, results that could resolve conflicting interpretations of experimental results. Instead of the previously proposed RuV[double bond, length as m-dash]O species, a proposed key species in water oxidation catalysis, we propose that the third oxidation of RuII could lead to a RuIIIN-oxide, via rapid water coordination followed by oxidation and concerted N-O bond formation.

Origin of the Absorption Band of Bromophenol Blue in Acidic and Basic pH: Insight from a Combined Molecular Dynamics and TD-DFT/MM Study.

We study the linear and nonlinear optical properties of a well-known acid-base indicator, bromophenol blue (BPB), in aqueous solution by employing static and integrated approaches. In the static approach, optical properties have been calculated using time-dependent density functional theory (TD-DFT) on the fully relaxed geometries of the neutral and different unprotonated forms of BPB. Moreover, both closed and open forms of BPB were considered. In the integrated approach, the optical properties have been computed over many snapshots extracted from molecular dynamics simulation using a hybrid time-dependent density functional theory/molecular mechanics approach. The static approach suggests closed neutral ⇒ anionic interconversion as the dominant mechanism for the red shift in the absorption spectra of BPB due to a change from acidic to basic pH. It is found by employing an integrated approach that the two interconversions, namely open neutral ⇒ anionic and open neutral ⇒ dianionic, can contribute to the pH-dependent shift in the absorption spectra of BPB. Even though both static and integrated approaches reproduce the pH-dependent red shift in the absorption spectra of BPB, the latter one is suitable to determine both the spectra and spectral broadening. Finally, the computed static first hyperpolarizability for various protonated and deprotonated forms of BPB reveals that this molecule can be used as a nonlinear optical probe for pH sensing in addition to its highly exploited use as an optical probe.

Optical Properties of Gold Nanoclusters Functionalized with a Small Organic Compound: Modeling by an Integrated Quantum-Classical Approach.

Motivated by the growing importance of organometallic nanostructured materials and nanoparticles as microscopic devices for diagnostic and sensing applications, and by the recent considerable development in the simulation of such materials, we here choose a prototype system - para-nitroaniline (pNA) on gold nanoparticles - to demonstrate effective strategies for designing metal nanoparticles with organic conjugates from fundamental principles. We investigated the motion, adsorption mode, and physical chemistry properties of gold-pNA particles, increasing in size, through classical molecular dynamics (MD) simulations in connection with quantum chemistry (QC) calculations. We apply the quantum mechanics-capacitance molecular mechanics method [Z. Rinkevicius et al. J. Chem. Theory Comput. 2014, 10, 989] for calculations of the properties of the conjugate nanoparticles, where time dependent density functional theory is used for the QM part and a capacitance-polarizability parametrization of the MM part, where induced dipoles and charges by metallic charge transfer are considered. Dispersion and short-range repulsion forces are included as well. The scheme is applied to one- and two-photon absorption of gold-pNA clusters increasing in size toward the nanometer scale. Charge imaging of the surface introduces red-shifts both because of altered excitation energy dependence and variation of the relative intensity of the inherent states making up for the total band profile. For the smaller nanoparticles the difference in the crystal facets are important for the spectral outcome which is also influenced by the surrounding MM environment.

Hybrid Complex Polarization Propagator/Molecular Mechanics Method for Heterogeneous Environments.

We introduce a hybrid complex polarization propagator/molecular mechanics method for the calculation of near-resonant and resonant response properties of molecules in heterogeneous environments, which consist of a metallic surface, or nanoparticle, and a solvent. The applicability and performance of the method is demonstrated by computations of linear absorption spectra of p-nitroaniline physisorbed at a gold/dimethyl sulfoxide interface in the UV/vis and near-carbon-K-edge regions of the spectrum. It is shown that the shift of absorption cross-section induced by the heterogeneous environment varies significantly depending on the nature of the excited states encountered in the targeted frequency region as well as on the actual size of the resonant frequencies, and that the solvent component of the heterogeneous environment is responsible for the major part of the environmental shift, especially in the higher frequency range of the carbon K-edge region.

Toward fully nonempirical simulations of optical band shapes of molecules in solution: a case study of heterocyclic ketoimine difluoroborates.

This study demonstrates that a hybrid density functional theory/molecular mechanics approach can be successfully combined with time-dependent wavepacket approach to predict the shape of optical bands for molecules in solutions, including vibrational fine structure. A key step in this treatment is the estimation of the inhomogeneous broadening based on the hybrid approach, where the polarization between solute and atomically decomposed solvent is taken into account in a self-consistent manner. The potential of this approach is shown by predicting optical absorption bands for three heterocyclic ketoimine difluoroborates in solution.

A time-dependent density-functional theory and complete active space self-consistent field method study of vibronic absorption and emission spectra of coumarin.

Time-dependent density-functional theory (TD-DFT) and complete active space multiconfiguration self-consistent field (CASSCF) calculations have been used to determine equilibrium structures and vibrational frequencies of the ground state and several singlet low-lying excited states of coumarin. Vertical and adiabatic transition energies of S1, S2, and S3 have been estimated by TD-B3LYP and CASSCF/PT2. Calculations predict that the dipole-allowed S1 and S3 states have a character of (1)(ππ*), while the dipole-forbidden (1)(nπ*) state is responsible for S2. The vibronic absorption and emission spectra of coumarin have been simulated by TD-B3LYP and CASSCF calculations within the Franck-Condon approximation, respectively. The simulated vibronic spectra show good agreement with the experimental observations available, which allow us to reasonably interpret vibronic features in the S0→S1 and S0→S3 absorption and the S0←S1 emission spectra. Based on the calculated results, activity, intensity, and density of the vibronic transitions and their contribution to the experimental spectrum profile have been discussed.

Non-linear optical properties of molecules in heterogeneous environments: a quadratic density functional/molecular mechanics response theory.

We generalize a density functional theory/molecular mechanics approach for heterogeneous environments with an implementation of quadratic response theory. The updated methodology allows us to address a variety of non-linear optical, magnetic and mixed properties of molecular species in complex environments, such as combined metallic, solvent and confined organic environments. Illustrating calculations of para-nitroaniline on gold surfaces and in solution reveals a number of aspects that come into play when analyzing second harmonic generation of such systems--such as surface charge flow, coupled surface-solvent dynamics and induced geometric and electronic structure effects of the adsorbate. Some ramifications of the methodology for applied studies are discussed.

A Hybrid Density Functional Theory/Molecular Mechanics Approach for Linear Response Properties in Heterogeneous Environments.

We introduce a density functional theory/molecular mechanical approach for computation of linear response properties of molecules in heterogeneous environments, such as metal surfaces or nanoparticles embedded in solvents. The heterogeneous embedding environment, consisting from metallic and nonmetallic parts, is described by combined force fields, where conventional force fields are used for the nonmetallic part and capacitance-polarization-based force fields are used for the metallic part. The presented approach enables studies of properties and spectra of systems embedded in or placed at arbitrary shaped metallic surfaces, clusters, or nanoparticles. The capability and performance of the proposed approach is illustrated by sample calculations of optical absorption spectra of thymidine absorbed on gold surfaces in an aqueous environment, where we study how different organizations of the gold surface and how the combined, nonadditive effect of the two environments is reflected in the optical absorption spectrum.