The dielectric function profile across the water interface through surface-specific vibrational spectroscopy and simulations.

The dielectric properties of interfacial water on subnanometer length scales govern chemical reactions, carrier transfer, and ion transport at interfaces. Yet, the nature of the interfacial dielectric function has remained under debate as it is challenging to access the interfacial dielectric with subnanometer resolution. Here we use the vibrational response of interfacial water molecules probed using surface-specific sum-frequency generation (SFG) spectra to obtain exquisite depth resolution. Different responses originate from water molecules at different depths and report back on the local interfacial dielectric environment via their spectral amplitudes. From experimental and simulated SFG spectra at the air/water interface, we find that the interfacial dielectric constant changes drastically across an ∼1 Šthin interfacial water region. The strong gradient of the interfacial dielectric constant leads, at charged planar interfaces, to the formation of an electric triple layer that goes beyond the standard double-layer model.

Accurate Determination of the Uniaxial Complex Refractive Index and the Optical Band Gap of Polymer Thin Films to Correlate Their Absorption Strength and Onset of Absorption.

The advanced development of optoelectronic devices requires a methodical knowledge of the fundamental material properties of the key active components. Systematic investigations and correlations of such basic optical properties can lead to new insights for the design of more potent materials. In this perspective, we provide a systematic overview of the uniaxial anisotropic complex refractive indices and the absorption coefficients obtained by ellipsometry as well as the optical band gap energies derived from Tauc plots of six selected solution-processed polymer thin films. While the optical band gap energies are intentionally distributed over the visible spectral range, we found that the absorption strength of all polymer samples are grouped in a random distribution within a rather uniform range of values.

Ultimate Limit for Optical Losses in Gold, Revealed by Quantitative Near-Field Microscopy.

We report thorough measurements of surface plasmon polaritons (SPPs) running along nearly perfect air-gold interfaces formed by atomically flat surfaces of chemically synthesized gold monocrystals. By means of amplitude- and phase-resolved near-field microscopy, we obtain their propagation length and effective mode index at visible wavelengths (532, 594, 632.8, 729, and 800 nm). The measured values are compared with the values obtained from the dielectric functions of gold that are reported in literature. Importantly, a reported dielectric function of monocrystalline gold implies ∼1.5 times shorter propagation lengths than those observed in our experiments, whereas a dielectric function reported for properly fabricated polycrystalline gold leads to SPP propagation lengths matching our results. We argue that the SPP propagation lengths measured in our experiments signify the ultimate limit of optical losses in gold, encouraging further comprehensive characterization of optical material properties of pure gold as well as other plasmonic materials.

Analysis of Dielectric Parameters of Fe2O3-Doped Polyvinylidene Fluoride/Poly(methyl methacrylate) Blend Composites.

In this paper, we report the effect of metal oxide (Fe2O3) loading in different weight ratios (0.5%, 1%, 2%, and 4%) on the structural and electrical parameters, viz., the complex dielectric constant, electric modulus spectra, and the AC conductivity, of polymeric composites of PVDF/PMMA (30/70 weight ratio) blend. The structural and geometric measurements have been analyzed with the help of peak location, peak intensity, and peak shape obtained from XRD as well as from FTIR spectra. The electrical properties have been investigated using an impedance analyzer in the frequency range 100 Hz to 1 MHz. The real parts of the complex permittivity and the dielectric loss tangent of these materials are found to be frequency independent in the range from 20 KHz to 1 MHz, but they increase with the increase in the concentration of nano-Fe2O3. The conductivity also increases with an increased loading of Fe2O3 in PVDF/PMMA polymer blends. The electric modulus spectra were used to analyze the relaxation processes associated with the Maxwell-Wagner-Sillars mechanism and chain segmental motion in the polymer mix.

Dual-Mode Infrared Absorption by Segregating Dopants within Plasmonic Semiconductor Nanocrystals.

When aliovalent dopants are sufficiently segregated to the core or near the surface of semiconductor nanocrystals, charge carriers donated by the dopants are also segregated to the core or near the surface, respectively. In Sn-doped indium oxide nanocrystals, we find that this contrast in free charge carrier concentration creates a core and shell with differing dielectric properties and results in two distinctly observable plasmonic extinction peaks. The trends in this dual-mode optical response with shell growth differ from core/shell nanoparticles composed of traditional plasmonic metals such as Au and Ag. We developed a model employing a core/shell effective medium approximation that can fit the dual-mode spectra and explain the trends in the extinction response. Lastly, we show that dopant segregation can improve sensitivity of plasmon spectra to changes in refractive index of the surrounding environment.

A new model dielectric function for loss functions and electron damping.

Trends in the zeroth frequency moment of the imaginary part of the dielectric function are studied for a wide range of metals, semiconductors and insulators. These results are combined with estimates for the inverse-first moment (related by Kramers-Kronig relations to the static dielectric function) and knowledge of the first moment from the f-sum rule. Matching all three moments allows for construction of a model dielectric function that reasonably predicts the loss function at different values of momentum and lifetime damping effects on occupied and unoccupied electron states. This is demonstrated by comparing model results and results of detailed, first-principles calculations.

Opto-Mechanics Driven Fast Martensitic Transition in Two-Dimensional Materials.

Diffusional phase-change materials, such as Ge-Sb-Te alloys, are used in rewritable nonvolatile memory devices. But the continuous pursuit of readout/write speed and reduced energy consumption in miniaturized devices calls for an optically driven, diffusionless phase change scheme in ultrathin materials. Inspired by optical tweezers, in this work, we illustrate theoretically and computationally that a linearly polarized laser pulse with selected frequency can drive an ultrafast diffusionless martensitic phase transition of two-dimensional ferroelastic materials such as SnO and SnSe monolayers, where the unit-cell strain is tweezed as a generalized coordinate that affects the anisotropic dielectric function and electromagnetic energy density. At laser power of 2.0 × 1010 and 7.7 × 109 W/cm2, the transition potential energy barrier vanishes between two 90°-orientation variants of ferroelastic SnO and SnSe monolayer, respectively, so displacive domain switching can occur within picoseconds. The estimated adiabatic thermal limit of energy input in such optomechanical martensitic transition (OMT) is at least 2 orders of magnitude lower than that in Ge-Sb-Te alloy.

Nanoscale Hydrogenography on Single Magnesium Nanoparticles.

Active plasmonics is enabling novel devices such as switchable metasurfaces for active beam steering or dynamic holography. Magnesium with its particle plasmon resonances in the visible spectral range is an ideal material for this technology. Upon hydrogenation, metallic magnesium switches reversibly into dielectric magnesium hydride (MgH2), turning the plasmonic resonances off and on. However, up until now, it has been unknown how exactly the hydrogenation process progresses in the individual plasmonic nanoparticles. Here, we introduce a new method, namely nanoscale hydrogenography, that combines near-field scattering microscopy, atomic force microscopy, and single-particle far-field spectroscopy to visualize the hydrogen absorption process in single Mg nanodisks. Using this method, we reveal that hydrogen progresses along individual single-crystalline nanocrystallites within the nanostructure. We are able to monitor the spatially resolved forward and backward switching of the phase transitions of several individual nanoparticles, demonstrating differences and similarities of that process. Our method lays the foundations for gaining a better understanding of hydrogen diffusion in metal nanoparticles and for improving future active nano-optical switching devices.

DelPhiPKa: Including salt in the calculations and enabling polar residues to titrate.

DelPhiPKa is a widely used and unique approach to compute pKa 's of ionizable groups that does not require molecular surface to be defined. Instead, it uses smooth Gaussian-based dielectric function to treat computational space via Poisson-Boltzmann equation (PBE). Here, we report an expansion of DelPhiPKa functionality to enable inclusion of salt in the modeling protocol. The method considers the salt mobile ions in solvent phase without defining solute-solvent boundary. Instead, the ions are penalized to enter solute interior via a desolvation penalty term in the Boltzmann factor in the framework of PBE. Hence, the concentration of ions near the protein is balanced by the desolvation penalty and electrostatic interactions. The study reveals that correlation between experimental and calculated pKa 's is improved significantly by taking into consideration the presence of salt. Furthermore, it is demonstrated that DelphiPKa reproduces the salt sensitivity of experimentally measured pKa 's. Another new development of DelPhiPKa allows for computing the pKa 's of polar residues such as cysteine, serine, threonine and tyrosine. With this regard, DelPhiPKa is benchmarked against experimentally measured cysteine and tyrosine pKa 's and for cysteine it is shown to outperform other existing methods (DelPhiPKa RMSD of 1.73 vs RMSD between 2.40 and 4.72 obtained by other existing pKa prediction methods).

Surface Plasmon Resonance Based Measurement of the Dielectric Function of a Thin Metal Film.

A spectral method based on surface plasmon resonance (SPR) in air is used to measure the dielectric function of a thin metal film. The method utilizes the spectral dependence of the ratio of the reflectances of p- and s-polarized waves measured in the Kretschmann configuration at different angles of incidence. By processing these dependences in the vicinity of a dip, or equivalently near the resonance wavelength, and using the dispersion characteristics of a metal film according to a proposed physical model, the real and imaginary parts of the dielectric function of the metal can be determined. The corresponding dielectric function of the metal is obtained by a least squares method for such a thickness minimizing the difference between the measured and theoretical dependence of the resonance wavelength on the the angle of incidence. The feasibility of the method is demonstrated in measuring the dielectric function of a gold film of an SPR structure comprising an SF10 glass prism and a gold coated SF10 slide with an adhesion film of chromium. The dielectric function according to the Drude⁻Lorentz model with two additional Lorentzian terms was determined in a wavelength range from 534 to 908 nm, and the results show that the gold film is composed of homogenous and rough layers with thicknesses 42.8 nm and 2.0 nm, respectively. This method is particularly useful in measuring the thickness and dielectric function of a thin metal film of SPR structures, directly in the Kretschmann configuration.

Electrostatic component of binding energy: Interpreting predictions from poisson-boltzmann equation and modeling protocols.

Macromolecular interactions are essential for understanding numerous biological processes and are typically characterized by the binding free energy. Important component of the binding free energy is the electrostatics, which is frequently modeled via the solutions of the Poisson-Boltzmann Equations (PBE). However, numerous works have shown that the electrostatic component (ΔΔGelec ) of binding free energy is very sensitive to the parameters used and modeling protocol. This prompted some researchers to question the robustness of PBE in predicting ΔΔGelec . We argue that the sensitivity of the absolute ΔΔGelec calculated with PBE using different input parameters and definitions does not indicate PBE deficiency, rather this is what should be expected. We show how the apparent sensitivity should be interpreted in terms of the underlying changes in several numerous and physical parameters. We demonstrate that PBE approach is robust within each considered force field (CHARMM-27, AMBER-94, and OPLS-AA) once the corresponding structures are energy minimized. This observation holds despite of using two different molecular surface definitions, pointing again that PBE delivers consistent results within particular force field. The fact that PBE delivered ΔΔGelec values may differ if calculated with different modeling protocols is not a deficiency of PBE, but natural results of the differences of the force field parameters and potential functions for energy minimization. In addition, while the absolute ΔΔGelec values calculated with different force field differ, their ordering remains practically the same allowing for consistent ranking despite of the force field used. © 2016 Wiley Periodicals, Inc.

Dielectric Genome of van der Waals Heterostructures.

Vertical stacking of two-dimensional (2D) crystals, such as graphene and hexagonal boron nitride, has recently lead to a new class of materials known as van der Waals heterostructures (vdWHs) with unique and highly tunable electronic properties. Ab initio calculations should in principle provide a powerful tool for modeling and guiding the design of vdWHs, but in their traditional form such calculations are only feasible for commensurable structures with a few layers. Here we show that the dielectric properties of realistic, incommensurable vdWHs comprising hundreds of layers can be efficiently calculated using a multiscale approach where the dielectric functions of the individual layers (the dielectric building blocks) are computed ab initio and coupled together via the long-range Coulomb interaction. We use the method to illustrate the 2D-3D transition of the dielectric function of multilayer MoS2 crystals, the hybridization of quantum plasmons in thick graphene/hBN heterostructures, and to demonstrate the intricate effect of substrate screening on the non-Rydberg exciton series in supported WS2. The dielectric building blocks for a variety of 2D crystals are available in an open database together with the software for solving the coupled electrodynamic equations.

pKa predictions for proteins, RNAs, and DNAs with the Gaussian dielectric function using DelPhi pKa.

We developed a Poisson-Boltzmann based approach to calculate the pKa values of protein ionizable residues (Glu, Asp, His, Lys and Arg), nucleotides of RNA and single stranded DNA. Two novel features were utilized: the dielectric properties of the macromolecules and water phase were modeled via the smooth Gaussian-based dielectric function in DelPhi and the corresponding electrostatic energies were calculated without defining the molecular surface. We tested the algorithm by calculating pKa values for more than 300 residues from 32 proteins from the PPD dataset and achieved an overall RMSD of 0.77. Particularly, the RMSD of 0.55 was achieved for surface residues, while the RMSD of 1.1 for buried residues. The approach was also found capable of capturing the large pKa shifts of various single point mutations in staphylococcal nuclease (SNase) from pKa-cooperative dataset, resulting in an overall RMSD of 1.6 for this set of pKa's. Investigations showed that predictions for most of buried mutant residues of SNase could be improved by using higher dielectric constant values. Furthermore, an option to generate different hydrogen positions also improves pKa predictions for buried carboxyl residues. Finally, the pKa calculations on two RNAs demonstrated the capability of this approach for other types of biomolecules.

Detailed study of the dielectric function of a lysozyme solution studied with molecular dynamics simulations.

The spread of microwave technology and new microwave applications in medicine have revitalized interest in the dielectric behavior of biological systems. In this work, the Fröhlich-Kirkwood approach and the linear response theory have been applied in conjunction with molecular dynamics simulations to study the dielectric response of a lysozyme solution as a model. The overall experimental dielectric behavior of a 9.88 mM lysozyme solution has been reproduced in a quantitative manner by employing a method based on the decomposition of the hydration shells close to the solute. Detailed analysis of the calculated spectra identified two δ-processes located at 200 MHz (δ1) and about 1 GHz (δ2), respectively. δ1 is associated mainly with the first hydration shell, while δ2 mainly with bulk water and the second hydration shell. Moreover, indications for the existence of an even faster relaxation in the 10(11)-Hz frequency range were found for the first time. Finally, the static dielectric constants of lysozyme and its first and second hydration shells were calculated based on the Fröhlich-Kirkwood and the linear response theory approaches.

Unveiling Influence of Dielectric Losses on the Localized Surface Plasmon Resonance in (Al,Ga)As:Sb Metamaterials.

We perform numerical modeling of the optical absorption spectra of metamaterials composed of systems of semimetal antimony nanoparticles embedded into AlxGa1-xAs semiconductor matrices. We reveal a localized surface plasmon resonance (LSPR) in these metamaterials, which results in a strong optical extinction band below, near, or above the direct band gap of the semiconductor matrices, depending on the chemical composition of the solid solutions. We elucidate the role of dielectric losses in AlxGa1-xAs, which impact the LSPR and cause non-plasmonic optical absorption. It appears that even a dilute system of plasmonic Sb nanoinclusions can substantially change the optical absorption spectra of the medium.

Dispersion Analysis of Perpendicular Modes Using a Hybrid Two-Trace Two-Dimensional (2T2D) Smart Error Sum.

Using linear dichroism theory, one would assume that a z-cut of a uniaxial crystal is equivalent to an x-cut to determine the perpendicular component of the dielectric tensor and the corresponding oscillator parameters. However, Fresnel's equations show that the effect of interfaces in the form of the continuity relations of the different components of the electric field must be considered. A consequence of the continuity relations is that perpendicular modes increase less significantly in strength with increasing angle of incidence than expected. This is a consequence of the fact that it is the inverse of the perpendicular component of the dielectric function that increasingly becomes important with a growing angle of incidence. An inverse dielectric function, however, has typically much smaller values than the dielectric function. An additional consequence is that perpendicular modes are blueshifted and coupled in such a way that oscillator strength is transferred to the higher wavenumber mode. Thus, the spectral signatures of perpendicular modes are often weak and masked by the parallel modes when two modes overlap. Accordingly, to enable dispersion analysis, it is suggested to use a hybrid of the conventional residual sum of squares and the two-trace two-dimensional (2T2D) smart error sum, which can correct systematic multiplicable errors in the experimental spectrum. As demonstrated for fresnoite (Ba2TiSi2O8), this is an important step toward determining the perpendicular component of the dielectric tensor and the corresponding oscillator parameters using dispersion analysis, since asynchronous 2T2D correlation spectra are, in particular, sensitive to perpendicular modes.

A Systematic Study of the Temperature Dependence of the Dielectric Function of GaSe Uniaxial Crystals from 27 to 300 K.

We report the temperature dependences of the dielectric function ε = ε1 + iε2 and critical point (CP) energies of the uniaxial crystal GaSe in the spectral energy region from 0.74 to 6.42 eV and at temperatures from 27 to 300 K using spectroscopic ellipsometry. The fundamental bandgap and strong exciton effect near 2.1 eV are detected only in the c-direction, which is perpendicular to the cleavage plane of the crystal. The temperature dependences of the CP energies were determined by fitting the data to the phenomenological expression that incorporates the Bose-Einstein statistical factor and the temperature coefficient to describe the electron-phonon interaction. To determine the origin of this anisotropy, we perform first-principles calculations using the mBJ method for bandgap correction. The results clearly demonstrate that the anisotropic dielectric characteristics can be directly attributed to the inherent anisotropy of p orbitals. More specifically, this prominent excitonic feature and fundamental bandgap are derived from the band-to-band transition between s and pz orbitals at the Γ-point.

Plasmon Resonance in a System of Bi Nanoparticles Embedded into (Al,Ga)As Matrix.

We reveal the feasibility of the localized surface plasmon resonance in a system of Bi nanoparticles embedded into an AlxGa1-xAs semiconductor matrix. With an ab initio determined dielectric function for bismuth and well-known dielectric properties of AlxGa1-xAs solid solution, we performed calculations of the optical extinction spectra for such metamaterial using Mie's theory. The calculations demonstrate a strong band of the optical extinction using the localized surface plasmons near a photon energy of 2.5 eV. For the semiconducting matrices with a high aluminum content x>0.7, the extinction by plasmonic nanoparticles plays the dominant role in the optical properties of the medium near the resonance photon energy.

Prediction of Frequency-Dependent Optical Spectrum for Solid Materials: A Multioutput and Multifidelity Machine Learning Approach.

The frequency-dependent optical spectrum is pivotal for a broad range of applications from material characterization to optoelectronics and energy harvesting. Data-driven surrogate models, trained on density functional theory (DFT) data, have effectively alleviated the scalability limitations of DFT while preserving its chemical accuracy, expediting material discovery. However, prevailing machine learning (ML) efforts often focus on scalar properties such as the band gap, overlooking the complexities of optical spectra. In this work, we employ deep graph neural networks (GNNs) to predict the frequency-dependent complex-valued dielectric function across the infrared, visible, and ultraviolet spectra directly from the crystal structures. We explore multiple architectures for the spectral multioutput representation of the dielectric function and utilize various multifidelity learning strategies, such as transfer learning and fidelity embedding, to address the challenges associated with the scarcity of high-fidelity DFT data. Additionally, we model key solar cell absorption efficiency metrics, demonstrating that learning these parameters is enhanced when integrated through a learning bias within the learning of the frequency-dependent absorption coefficient. This study demonstrates that leveraging multioutput and multifidelity ML techniques enables accurate predictions of optical spectra from crystal structures, providing a versatile tool for rapidly screening materials for optoelectronics, optical sensing, and solar energy applications across an extensive frequency spectrum.

Absorption in narrow and wide gap materials.

The absorption spectrum of a material reveals the absorbed light frequencies, characteristic peaks, and the line width of absorption bands. This information is critical for understanding the energy levels involved in the absorption process as well as the material's electronic structure. In this study, an equation connecting the absorption line width with the static dielectric function is derived for narrow and wide gap materials. It is then compared with the Penn model. It has been found that the constant in the Penn model has a value that is restricted to the range of 0.5 to 1. Application of this equation to various narrow and wide gap materials is then discussed.