Extrapolation of the linear and nonlinear polarizabilities from ab initio finite oligomer calculations

Different definitions of property per unit cell and different fitting functions are employed to obtain the asymptotic limit values per unit cell of the polarizability (α), the first (β), and the second (γ) hyperpolarizabilities of an infinite oligomer. A 1/n power series function is found to be suitable for the average value and logarithmic average value per unit cell definition, and an exponentially decreasing function is found to be suitable for the difference value per unit cell definition. These conclusions are derived based on an equation expressing the total energy per unit cell of a finite linear oligomer as a power series of 1/n, presented from a perturbation treatment. Several calculations of long chain systems have been carried out to reach our conclusions. An equation of p(n)/n = a + b/n + c/n² is strongly recommended for a least‐squares fitting of the properties per unit cell to achieve a stabilization behavior when the chain length is increased. © 2006 Wiley Periodicals, Inc. Int J Quantum Chem, 2006     

Evaluation of the molecular static and dynamic first hyperpolarizabilities

By considering two prototypical π‐conjugated compounds, several technical aspects associated with the evaluation of the first hyperpolarizabilities have been addressed in this article, that is, (i) the automatization of the Romberg's scheme to improve the numerical accuracy in the finite field method, (ii) the evaluation of the frequency dispersion at correlated levels using approximate schemes, and (iii) the deviations from Kleinman's symmetry conditions. It results from this study that accurate numerical derivatives can be obtained by resorting to the Romberg's method and by analyzing the Romberg's table in terms of two quantities, the field error and the iteration error. Indeed, the resulting first hyperpolarizability values are in close agreement with those obtained using an analytical differentiation procedure. The reliability of the multiplicative and additive approximate schemes to describe the frequency dispersion at correlated levels from using HF (Hartree‐Fock) frequency dispersion has been confirmed to be limited to large wavelengths or far‐from‐resonance wavelength regions. Kleinman's symmetry conditions have been assessed, showing that for off‐diagonal components of these two π‐conjugated compounds, the deviations could be substantial. Nevertheless, good accuracy can be achieved for experimentally related quantities like β_HRS because the diagonal tensor components are dominant. © 2014 Wiley Periodicals, Inc.     

Intensity-carrying modes important for vibrational polarizabilities and hyperpolarizabilities of molecules: derivation from the algebraic properties of formulas and applications

The intensity-carrying mode (ICM) theory is developed for analyzing the vibrational motions that mainly contribute to vibrational polarizabilities and hyperpolarizabilities, which are important for describing intermolecular electrostatic interactions and nonlinear optical properties of molecules. The ICMs are derived from dipole derivatives, polarizability derivatives, and first hyperpolarizability derivatives by using algebraic properties of intensity formulas. The way to obtain explicit forms of ICMs, including the optimization method of the basis of the ICM vector space, is discussed in detail. One- and two-dimensional models are constructed on the basis of the ICMs. The theory is applied to three molecules (a push-pull type polyene, a streptocyanine dye cation, and a symmetric neutral polyene) taken as typical examples. It is shown that the ICM theory provides a reasonable picture on the vibrational polarization properties of these molecules. On the basis of this result, the validity of the valence-bond charge transfer (VB-CT) model, which is a one-dimensional model and is widely used to describe the electronic and vibrational properties of dye molecules, is also discussed.     

Implementation of the finite field perturbation method in the CRYSTAL program for calculating the dielectric constant of periodic systems

The finite field approach has been implemented in the periodic ab initio CRYSTAL program and been used for calculating the dielectric constants of crystalline LiF and MgO (FCC structure) and BeO (wurtzite structure). To maintain the periodicity along the applied field direction, a "sawtooth" potential is used in conjunction with a supercell scheme. Supercells four to five times longer than the primitive cell in the direction of the applied field provide well-converged results. The influence of the computational parameters is discussed. An alternative scheme has also been implemented, for inner check, that consists of applying a static electric field to a slab of increasing thickness in the direction orthogonal to the surface; the dielectric response at the center of the slab is shown to converge rapidly to the bulk value evaluated with the sawtooth field. The method is accurate and permits the determination of nonlinear corrections to the dielectric constant. When used in conjunction with the local density approximation (LDA) scheme, it provides for the dielectric constant of the three above-mentioned compounds values close to those recently obtained with a time-dependent density functional theory approach.     

香草醛晶体光学非线性系数的有限场方法研究

采用有限场方法计算了有机非线性材料3-甲氧基-4-羟基-苯甲醛(分子式C~8H~8O~3,简称MHBA)的二阶非线性光学系数d~2~1,d~2~2,d~2~3,d~2~5。研究结果表明在MP2近似和6-31G+(p,d)基组水平上理论值能较好地与实验值符合,同时也表明拉电子基团HC＝O和推电子基团HO对MHBA分子的光学非线性性能的影响很大,而晶体中分子间的氢键作用对d值的贡献很小。     

Finite‐field method with unbiased polarizable continuum model for evaluation of the second hyperpolarizability of an open‐shell singlet molecule in solvents

The static second hyperpolarizability γ of the complexes composed of open‐shell singlet 1,3‐dipole molecule involving a boron atom and a water molecule in aqueous phase are investigated by the finite‐field (FF) method combined with a standard polarized continuum model (PCM) and with a newly proposed unbiased PCM (UBPCM). On the basis of the comparison with the results calculated by the FF method using the full quantum and the quantum‐mechanical/molecular‐mechanical and molecular‐dynamics (QM/MM‐MD) treatments, the present FF‐UBPCM method is demonstrated to remedy the artificial overestimation of the γ caused by standard FF‐PCM calculations and to well reproduce the FF‐QM/MM‐MD and FF‐full‐QM results with much lower costs. © 2013 Wiley Periodicals, Inc.     

Numerical simulation of deformation behavior of droplet in gas under the electric field and flow field coupling

The water droplets in the process of electrostatic coalescence are important when studying electrohydrodynamics. In the present study, the electric field and flow field are coupled through the phase field method based on the Cahn–Hilliard formulation. A numerical simulation model of single droplet deformation under the coupling field was established. It simulated the deformation behavior of the movement of a droplet in the continuous phase and took the impact of droplet deformation into consideration which is affected by two-phase flow velocity, electric field strength, the droplet diameter, and the interfacial tension. The results indicated that under the single action of the flow field, when the flow velocity was lower, the droplet diameter was greater as was the droplet deformation degree. When the flow velocity was increased, the droplet deformation degree of a small-diameter droplet was at its maximum size, the large-diameter droplet had a smaller deformation degree, and the middle-diameter droplet was at a minimum deformation degree. When the flow velocity was further increased, the droplet diameter was smaller, and the droplet deformation degree was greater. Under the coupled effect of the electric field and flow field, the two-phase flow velocity and the electric field strength were greater, and the degree of droplet deformation was greater. While the droplet diameter and interfacial tension were smaller, the degree of droplet deformation was greater. Droplet deformation degree increased along with the two-phase flow velocity. The research results provided a theoretical basis for gas–liquid separation with electrostatic coalescence technology.     

Numerical simulation of droplet coalescence behavior in gas phase under the coupling of electric field and flow field

The coalescence behavior of droplets in an electric field belongs to the important research contents of electrohydrodynamics. Based on the phase field method of the Cahn–Hilliard equation, the electric field and the flow field are coupled to establish the numerical model of twin droplet coalescence in a coupled field. The effects of flow rate, electric field strength, droplet diameter, and interfacial tension on the coalescence behavior of droplets during the coalescence process were investigated. The results show that the dynamic behavior of the droplets is divided into coalescence, after coalescence rupture, and no coalescence under the coupling of electric field and flow field. The proper increase of the electric field strength will accelerate the coalescence of the droplets, and the high electric field strength causes the droplets to burst after coalescence. Excessive flow rates make droplets less prone to coalescence. Under the coupling field, the larger the droplet interface tension, the smaller the droplet diameter, the smaller the flow rate, and the shorter the droplet coalescence time. The results provide a theoretical basis for the application of electrostatic coalescence in gas–liquid separation technology.     

Global optimization of parameters in the reactive force field ReaxFF for SiOH

We have used unbiased global optimization to fit a reactive force field to a given set of reference data. Specifically, we have employed genetic algorithms (GA) to fit ReaxFF to SiOH data, using an in‐house GA code that is parallelized across reference data items via the message‐passing interface (MPI). Details of GA tuning turn‐ed out to be far less important for global optimization efficiency than using suitable ranges within which the parameters are varied. To establish these ranges, either prior knowledge can be used or successive stages of GA optimizations, each building upon the best parameter vectors and ranges found in the previous stage. We have finally arrive‐ed at optimized force fields with smaller error measures than those published previously. Hence, this optimization approach will contribute to converting force‐field fitting from a specialist task to an everyday commodity, even for the more difficult case of reactive force fields. © 2013 Wiley Periodicals, Inc.     

Tunneling through a fluctuating barrier in the presence of a periodically driving field

The zero and finite temperature tunneling dynamics of a periodically driven particle moving in a bistable potential with a fluctuating barrier is studied. We have focused on the influence of barrier fluctuation and thermal modulation on the tunneling processes in the presence of a driving field. At zero temperature, for a fixed strength of the driving field, both the tunneling probability and rate passes through a well-defined minimum when plotted as a function of fluctuation frequency while it reveals a clear maximum as a function of driving frequency. However, at T > 0 the tunneling probability and rate show two maxima as a function of both fluctuation frequency and driving frequency. In both zero and finite temperature, the tunneling rate constant decreases with increasing fluctuation strength. So, the barrier fluctuation may enhance the stability of a periodically driven system. © 2003 Wiley Periodicals, Inc. Int J Quantum Chem, 2004     

Charge distribution from a simple molecular orbital type calculation and non-bonding interaction terms in the force field MAB

A simple and fast method to calculate charge distributions in organic molecules is presented. The method is based on charge shifts within the saturated -system, driven by orbital electronegativities, coupled to a modified Hückel treatment of the unsaturated -systems. Experimental molecular dipole moments of a set of 119 molecules are reproduced with a root mean square deviation of 0.36 Debye units. Furthermore, the obtained charge distribution is used to describe hydration free energies in terms of hydrogen-bonding donor and acceptor strengths of polar groups. Least square fitting to experimental data of 281 compounds leads to values for these strengths with accuracy limits of ±4.3% and ±2.5%, respectively. Properly normalized values are taken to parametrize the hydrogen bonding terms in our MAB force field. The method is sufficiently fast to be used in the preparatory phase of interactive force-field calculations.     

Numerical comparison between Maxwell stress method and equivalent multipole approach for calculation of the dielectrophoretic force in single-cell traps

This paper presents detailed numerical calculations of the dielectrophoretic force in traps designed for single-cell trapping. A trap with eight planar electrodes is studied for spherical and ellipsoidal particles using the boundary element method (BEM). Multipolar approximations of orders one to three are compared with the full Maxwell stress tensor (MST) calculation of the electrical force on spherical particles. Ellipsoidal particles are also studied, but in their case only the dipolar approximation is available for comparison with the MST solution. The results show that a small number of multipolar terms need to be considered in order to obtain accurate results for spheres, even in the proximity of the electrodes, and that the full MST calculation is only required in the study of non-spherical particles.     

Comparison of interpolation polynomials with divided differences,interpolation polynomials with finite differences,and quadratic functions obtained by the least squares method in modeling of chromatographic responses

A novel approach to mathematical modeling of chromatographic responses based on interpolation polynomials with divided differences and with finite differences is discussed. These interpolational techniques as well as traditionally applied second‐order polynomial models obtained by least squares are compared. Interpolation techniques can be useful in situations where commonly used linear or quadratic models are not applicable: when the nature of dependence is complex or the investigated factor intervals are broad. The three analyzed modeling techniques are incorporated in a design of experiments methodology for systematic development and optimization of liquid chromatographic methods. The direct modeling of retention factors is carried out first, while the objective function for final quality measurement is calculated last. An interpolation polynomial with divided differences resulted in a high quality fit compared with the results obtained by the other two modeling approaches and succeeded in locating the desired optimum. It is shown that this modeling technique can be a useful alternative for modeling of chromatographic responses. Copyright © 2013 John Wiley & Sons, Ltd.     

Metabolic phenotyping of human plasma by 1H‐NMR at high and medium magnetic field strengths: a case study for lung cancer

Accurate identification and quantification of human plasma metabolites can be challenging in crowded regions of the NMR spectrum with severe signal overlap. Therefore, this study describes metabolite spiking experiments on the basis of which the NMR spectrum can be rationally segmented into well‐defined integration regions, and this for spectrometers having magnetic field strengths corresponding to ¹H resonance frequencies of 400 MHz and 900 MHz. Subsequently, the integration data of a case–control dataset of 69 lung cancer patients and 74 controls were used to train a multivariate statistical classification model for both field strengths. In this way, the advantages/disadvantages of high versus medium magnetic field strength were evaluated. The discriminative power obtained from the data collected at the two magnetic field strengths is rather similar, i.e. a sensitivity and specificity of respectively 90 and 97% for the 400 MHz data versus 88 and 96% for the 900 MHz data. This shows that a medium‐field NMR spectrometer (400–600 MHz) is already sufficient to perform clinical metabolomics. However, the improved spectral resolution (reduced signal overlap) and signal‐to‐noise ratio of 900 MHz spectra yield more integration regions that represent a single metabolite. This will simplify the unraveling and understanding of the related, disease disturbed, biochemical pathways. Copyright © 2017 John Wiley & Sons, Ltd.     

New basis sets for the evaluation of interaction‐induced electric properties in hydrogen‐bonded complexes

Interaction‐induced static electric properties, that is, dipole moment, polarizability, and first hyperpolarizability, of the CO? (HF)_n and N₂? (HF)_n, n = 1–9 hydrogen‐bonded complexes are evaluated within the finite field approach using the Hartree–Fock, density functional theory, Møller–Plesset second‐order perturbation theory, and coupled cluster methods, and the LPol‐n (n = ds, dl, fs, fl) basis sets. To compare the performance of the different methods with respect to the increase of the complex size, we consider as model systems linear chains of the complexes. We analyze the results in terms of the many‐body and cooperative effects. © 2012 Wiley Periodicals, Inc.     

Boundary element solution of the linear Poisson-Boltzmann equation and a multipole method for the rapid calculation of forces on macromolecules in solution

The Poisson-Boltzmann equation is widely used to describe the electrostatic potential of molecules in an ionic solution that is treated as a continuous dielectric medium. The linearized form of this equation, applicable to many biologic macromolecules, may be solved using the boundary element method. A single-layer formulation of the boundary element method, which yields simpler integral equations than the direct formulations previously discussed in the literature, is given. It is shown that the electrostatic force and torque on a molecule may be calculated using its boundary element representation and also the polarization charge for two rigid molecules may be rapidly calculated using a noniterative scheme. An algorithm based on a fast adaptive multipole method is introduced to further increase the speed of the calculation. This method is particularly suited for Brownian dynamics or molecular dynamics simulations of large molecules, in which the electrostatic forces must be calculated for many different relative positions and orientations of the molecules. It has been implemented as a set of programs in C++, which are used to study the accuracy and speed of this method for two actin monomers.