Calculation of molecular response properties with the second‐order coupled perturbed electron propagator

We present the theory of the electron propagator perturbed by an external electric field and show how it can be used to calculate a variety of one‐electron linear response properties that are accurate through second order in electron correlation. Some illustrative calculations are discussed. © 2000 John Wiley & Sons, Inc. Int J Quant Chem 77: 625–640, 2000     

Theoretical studies on geometry,solvent effect,and photochromic mechanism of two bis‐heterocyclic compounds containing pyrazolone ring

Two bis‐heterocyclic compounds containing pyrazolone ring, 1‐phenyl‐3‐methyl‐4‐(6‐hydro‐4‐amino‐5‐sulfo‐2,3‐pyrazine)‐pyrazole‐5‐one and 1‐phenyl‐3‐methyl‐4‐(6‐hydro‐4‐methylamino‐5‐sulfo‐2,3‐pyrazine)‐pyrazole‐5‐one, are investigated to gain a deeper insight into their geometries and photochromic mechanism by applying density functional theory. The solvent effects are simulated using the polarizable continuum model of the self‐consistent reaction field theory. Bader's atom‐in‐molecule theory is used to investigate the nature of hydrogen bonds. The data of energy, dipole moments, and the condensed Fukui functions have been calculated to assess the stability and reactivity of the title compounds. © 2009 Wiley Periodicals, Inc. Int J Quantum Chem, 2010     

Multiresponse problems: desirability and other optimization approaches

A diversity of multiresponse optimization methods has been introduced in the literature; however, their performance has not been thoroughly explored, and only a classical desirability‐based criterion has been commonly used. With the aim of contributing to help practitioners in selecting an effective criterion for solving multiresponse optimization problems developed under the response surface methodology framework, and thus to find compromise solutions that are technically and economically more favorable, the working ability of several easy‐to‐use criteria is evaluated and compared with that of a theoretically sound method. Four case studies with different numbers and types of responses are considered. Less‐sophisticated criteria were able to generate solutions similar to those generated by sophisticated methods, even when the objective is to depict the Pareto frontier in problems with conflicting responses. Two easy‐to‐use criteria that require less‐subjective information from the user yielded solutions similar to those of a classical desirability‐based criterion. Preference parameters range and increment impact on optimal solutions were also evaluated.     

A statistical approach to delta layer depth profiling

We present a simple statistical model describing the removal and relocation of material during a sputter depth profiling experiment. All input parameters are determined from low‐fluence molecular dynamics simulations, making the model de facto parameter free. The model can be used to extrapolate data from the molecular dynamics simulations to projectile fluences relevant to sputter depth profiling experiments. As a result, the erosion of the surface is calculated in terms of fluence‐dependent filling factors of different sample layers. Using input data determined for the 20‐keV C₆₀ cluster bombardment of silicon, it is found that a steady‐state erosion profile is reached after removal of approximately 20 monolayer equivalents of material. Plotting the contribution of particles from a specific layer to the instantaneous sputtered flux, one can directly determine the delta layer response function predicted from such a model. It is shown that this function can be parameterized by the semiempirical Dowsett response function, and the resulting fitting parameters are compared with published depth profile data. The model is then used to study the role of different processes influencing the observed depth resolution. We find that the statistical nature of the sputtering process suffices to explain many features of experimentally measured delta layer depth profiles. Copyright © 2012 John Wiley & Sons, Ltd.     

Optimization of the ultrasound‐assisted extraction of antioxidant phloridzin from Lithocarpus polystachyus Rehd. using response surface methodology

The purpose of this study was to optimize the extraction process of phloridzin from Lithocarpus polystachyus Rehd. leaves using response surface methodology and to determine the antioxidant capacity of the extract. A Box–Behnken design was used to analyze the effects of ethanol concentration, liquid–solid ratio, soak time and extraction time on the extraction yield of phloridzin. The content of phloridzin was determined by high‐performance liquid chromatography. To assess the antioxidant capacity of the extract, three in vitro test systems were used (1,1‐,diphenyl‐2‐picrylhydrazyl, hydroxyl radical scavenging test and reduction force). The optimal parameters obtained by response surface methodology were a volume fraction of ethanol of 64%, a liquid–solid ratio of 37:1, a soaking time of 35 h and a sonication time of 38 min. The proportion of the extraction of phloridzin from L. polystachyus under these industrial process conditions was 3.83%. According to the obtained results, response surface methodology could be suggested as an adequate model for optimizing the extraction process of phloridzin from L. polystachyus . Ultrasound extraction significantly increased the extraction rate of phloridzin, which could be used as an antioxidant in pharmaceutical and food products.     

Desirability functions as response in a d‐optimal design for evaluating the extraction and purification steps of six tranquillizers and an anti‐adrenergic by liquid chromatography‐tandem mass spectrometry

Internal standards can be added at different stages of an analytical procedure. When they are added at the beginning of a multiresidue method and their behavior is not exactly the same as that of the analytes, the intended correction for small variations within the analytical process could not be achieved. Because of this, in the present work, the use of d ‐optimal designs together with desirability functions is proposed to state the experimental response under study. The overall desirability function used relates two analytical criteria: to assess a similar chemical behavior of each analyte in relation to its internal standard and to avoid a significant reduction of the absolute peak area of the internal standards. This strategy has been applied to the analysis of the effect of four factors related to the extraction and purification steps of six tranquillizers and a β‐blocker from pig muscle analyzed by liquid chromatography–tandem mass spectrometry. The effect of those factors has been evaluated by means of an ad hoc d ‐optimal design consisting of only 11 experiments. The resulting levels of the four factors that enable to achieve the greatest overall desirability have also been compared with those obtained when either the standardized or absolute peak area has been considered as response. Differences in both the significant factors and their optimum levels have been observed. It is noticeable that the experimental effort necessary to study the effect of the factors has been reduced by more than 50% thanks to the d ‐optimal design. Copyright © 2015 John Wiley & Sons, Ltd.     

Fluorescence Modulation in Tribranched Switchable [4]Rotaxanes

Two novel tribranched [4]rotaxanes with a 1,3,5‐triphenylene core and three rotaxane arms have been designed, synthesized, and characterized by ¹H and ¹³C NMR spectroscopies and HR‐ESI mass spectrometry. [4]Rotaxanes 1 and 2 each possess the same three‐armed skeleton. Each arm incorporates two distinguishable binding sites for a dibenzo[24]crown‐8 ring, namely a dibenzylammonium site and an N‐methyltriazolium site, and is terminated by a 4‐morpholino‐naphthalimide fluorophore as a stopper. [4]Rotaxane 1 has three di‐ferrocene‐functionalized dibenzo[24]crown‐8 rings whereas 2 has three simple dibenzo[24]crown‐8 rings interlocked with the thread component. Uniform shuttling motions of the three macrocycles in both 1 and 2 can be driven by external acid–base stimuli, which were confirmed by ¹H NMR spectroscopy. However, [4]rotaxanes 1 and 2 show distinct modes of fluorescence modulation in response to external acid–base stimuli. [4]Rotaxane 1 exhibits a remarkable fluorescence decrease in response to the addition of 1,8‐diazabicyclo[5.4.0]undec‐7‐ene (DBU) as a base, which can displace the ferrocene‐functionalized macrocycle from the dibenzylammonium station to the N‐methyltriazolium station. In contrast, the fluorescence intensity of [4]rotaxane 2 showed an enhancement with the addition of DBU. Time‐resolved fluorescence measurements have been performed. The different photoinduced electron‐transfer processes responsible for the fluorescence changes in the two molecular systems are discussed. Topological structures of this kind have significant potential for the design and construction of large and complex assemblies with controllable functions.     

A Refined Intermolecular Interaction Potential for Methane: Spectral Analysis and Molecular Dynamics Simulations

The previously constructed methane interaction potential energy surface calculated at the second‐order Møller‐Plesset (MP2) perturbation theory has been significantly improved in two aspects. First, all ab initio potential energy data are calculated by the supermolecule counterpoise corrected coupled cluster with single and double and perturbative triple excitations [CCSD(T)] method with Dunning’s correlation‐consistent aug‐cc‐pVXZ, X=D, T, Q, 5, basis sets and extrapolated to the complete basis set (CBS) limits with a convergence precision of 0.01 kcal/mol. Second, instead of the simple 4‐site model proposed in the previous study, a 5‐site model has been used to represent the ab initio potential data. The simulated infrared spectrum using the potential energy surface seems to be broadly in line with the spectral features observed in experiments. Molecular dynamics simulations using the ab initio force field show quantitative agreements with experiments. The properties examined in this paper include the atom‐to‐atom radial distribution functions in liquid and supercritical phases and the self‐diffusion coefficients over a wide range of thermodynamic conditions. It is shown that the refined ab initio force field can be applied to study fluid properties in different phases.     

Large Solvation Effect in the Optical Rotatory Dispersion of Norbornenone

The anomalously large chiroptical response of (1R,4R)‐norbornenone has been probed under complementary vapor‐phase and solution‐phase conditions to assess the putative roles of environmental perturbations. Measurements of the specific rotation for isolated (gas‐phase) molecules could not be reproduced quantitatively by comprehensive quantum‐chemical calculations based on density‐functional or coupled‐cluster levels of linear‐response theory, which suggests that higher‐order treatments may be needed to accurately predict such intrinsic behavior. A substantial, yet unexpected, dependence of the dispersive optical activity on the nature (phase) of the surrounding medium has been uncovered, with the venerable Lorentz local‐field correction reproducing solvent‐mediated trends in rotatory dispersion surprisingly well, whereas more modern polarizable continuum models for implicit solvation performed less satisfactorily.     

High‐resolution capillary zone electrophoresis for transferrin glycoform analysis associated with congenital disorders of glycosylation

High‐resolution capillary zone electrophoresis is used to assess the transferrin profile in serum of patients with eight different congenital disorders of glycosylation that represent type I, type II, and mixed type I/II disorders. Capillary zone electrophoresis data are compared to patterns obtained by gel isoelectric focusing. The high‐resolution capillary zone electrophoresis method is shown to represent an effective tool to assess the diversity of transferrin patterns. Hypoglycosylated disialo‐, monosialo‐, and asialo‐transferrin in type I cases can be distinguished from the corresponding underdesialylated transferrin glycoforms present in type II disorders. The latter can be separated from and detected ahead of their corresponding hypoglycosylated forms of type I patients. Both types of glycoforms are detected in sera of mixed type I/II patients. The assay has the potential to be used as screening method for congenital disorders of glycosylation. It can be run with a few microliters of serum when microvials are used.     

Atomic forces for geometry‐dependent point multipole and Gaussian multipole models

In standard treatments of atomic multipole models, interaction energies, total molecular forces, and total molecular torques are given for multipolar interactions between rigid molecules. However, if the molecules are assumed to be flexible, two additional multipolar atomic forces arise because of (1) the transfer of torque between neighboring atoms and (2) the dependence of multipole moment on internal geometry (bond lengths, bond angles, etc.) for geometry‐dependent multipole models. In this study, atomic force expressions for geometry‐dependent multipoles are presented for use in simulations of flexible molecules. The atomic forces are derived by first proposing a new general expression for Wigner function derivatives . The force equations can be applied to electrostatic models based on atomic point multipoles or Gaussian multipole charge density. Hydrogen‐bonded dimers are used to test the intermolecular electrostatic energies and atomic forces calculated by geometry‐dependent multipoles fit to the ab initio electrostatic potential. The electrostatic energies and forces are compared with their reference ab initio values. It is shown that both static and geometry‐dependent multipole models are able to reproduce total molecular forces and torques with respect to ab initio, whereas geometry‐dependent multipoles are needed to reproduce ab initio atomic forces. The expressions for atomic force can be used in simulations of flexible molecules with atomic multipoles. In addition, the results presented in this work should lead to further development of next generation force fields composed of geometry‐dependent multipole models. © 2010 Wiley Periodicals, Inc. J Comput Chem, 2010     

Optimization and validation of a nonaqueous micellar electrokinetic chromatography method for determination of polycyclic musks in perfumes

A nonaqueous micellar electrokinetic chromatography method was developed for determination of Tonalide®, Galaxolide®, and Traseolide® polycyclic musks (PCMs). These compounds are widely used as fragrance ingredients in cosmetics. The method was optimized by using a three variable Box‐Behnken experimental design and response surface methodology. A modified chromatographic response function was defined in order to adequately weigh the terms in the response function. After optimization of experimental conditions, an electrolyte solution of 195 mM SDS and 40 mM NaH₂PO₄ in formamide was selected for the separation of the three PCMs, and the applied voltage was fixed at 30 kV. The nonaqueous MEKC method was then checked in terms of linearity, limits of detection and quantification, repeatability, intermediate precision and accuracy, providing appropriate values (i.e. RSD values for precision never exceeding 7%, and accuracy 96–107%). Nonaqueous MEKC for determination of the selected compounds was successfully applied to the analysis of commercial perfume samples.     

Moment expansion of the linear density‐density response function

We present a low rank moment expansion of the linear density‐density response function. The general interacting (fully nonlocal) density‐density response function is calculated by means of its spectral decomposition via an iterative Lanczos diagonalization technique within linear density functional perturbation theory. We derive a unitary transformation in the space of the eigenfunctions yielding subspaces with well‐defined moments. This transformation generates the irreducible representations of the density‐density response function with respect to rotations within SO(3). This allows to separate the contributions to the electronic response density from different multipole moments of the perturbation. Our representation maximally condenses the physically relevant information of the density‐density response function required for intermolecular interactions, yielding a considerable reduction in dimensionality. We illustrate the performance and accuracy of our scheme by computing the electronic response density of a water molecule to a complex interaction potential. © 2015 Wiley Periodicals, Inc.     

Using multivariate statistical methods to model the electrospray ionization response of GXG tripeptides based on multiple physicochemical parameters

Response factors were determined for twelve GXG peptides (where G stands for glycine and X is any of alanine [A], arginine [R], asparagine [N], aspartic acid [D], glycine [G], histidine [H], leucine [L], lysine [K], phenylalanine [F], serine [S], tyrosine [Y], valine [V]) by electrospray ionization mass spectrometry (ESI‐MS). The response factors were measured using a novel flow injection method. This new method is based on the Gaussian distribution of analyte concentration resulting from band‐broadening dispersion experienced by the analyte upon passage through an extended volume of PEEK tubing. This method removes the need for preparing a discrete series of standard solutions to assess concentration‐dependent response. Relative response factors were calculated for each peptide with reference to GGG. The observed trends in the relative response factors were correlated with several analyte physicochemical parameters, chosen based on current understanding of ion release from charged droplets during the ESI process. These include analyte properties: nonpolar surface area; polar surface area; gas‐phase basicity; proton affinity; and Log D. Multivariate statistical analysis using multiple linear regression, decision tree, and support vector regression models were investigated to assess their potential for predicting ESI response based on the analyte properties. The support vector regression model was more versatile and produced the least predictive error following 12‐fold cross‐validation. The effect of variation in solution pH on the relative response factors is highlighted, as evidenced by the different predictive models obtained for peptide response at two pH values (pH = 6.0 and 9.0). The relationship between physicochemical parameters and associated ionization efficiencies for GXG tripeptides is discussed based on the equilibrium partitioning model. Copyright © 2009 John Wiley & Sons, Ltd.     

Green chromatographic fingerprinting: An environmentally friendly approach for the development of separation methods for fingerprinting complex matrices

A chromatographic fingerprint is a comprehensive method that reveals the distinctive pattern of peaks across the chromatogram for a given sample. It is considered an effective strategy to assess the identity and quality of herbal materials, as well as for the control of the quality of their derived products. HPLC is the most employed technique for these purposes and it is used routinely for quality control in industry. Hence, its impact on the environment should not be neglected. This work provides a rational and generic procedure to qualitatively fingerprint complex matrices. Resource‐ and time‐saving experimental designs were selected; an alternative safer organic solvent was tested and a time‐saving and innovative response entitled the green chromatographic fingerprinting response was developed and employed. This procedure was applied in the development of chromatographic fingerprints for extracts of Bauhinia forficata and Casearia sylvestris. Moreover, the response proposed here can be combined with a complementary metric available in the literature to compare methods using different solvents. According to this, the chromatographic fingerprints developed here using ethanol as the organic solvent provided a performance better than that of reference methods in which more harmful acetonitrile or methanol were employed.     

Role of resonances in building cross sections: Comparison between the Mittag–Leffler and the T‐matrix Green function expansion approaches

Peaks in collision cross sections are often interpreted as resonances. The complex dilation method, as well as other methods relying on analytic continuation of the scattering formalism, can be used to clarify whether these structures are true resonances in the sense that they are poles of the S‐matrix and the associated Green function. The performance of the Mittag–Leffler expansion and T‐matrix Green function expansion methods are formally and computationally compared. The two methods are applied to two model potentials. Eigenenergies, s‐wave residues, and cross sections are computed with both methods. The resonance contributions to the cross sections are further analyzed by removing the residue contributions from the Mittag–Leffler and Green function expansion sums, respectively. It is suggested that the contribution of a resonance to a cross section should be defined through its S‐matrix residue. © 2006 Wiley Periodicals, Inc. Int J Quantum Chem, 2007     

Evaluation of the van Deemter equation in terms of open‐ended flow to chromatography

Plate theory and adsorption theory are the main tools available for understanding chromatographic experiments. Both theories predict a Gaussian distribution of solute molecules within the tubular system. However, Gaussian concentration distributions are observed predominantly at slow linear flow rates, while asymmetric concentration distributions are observed at the linear flow rates most used in chromatography. Allegedly, this asymmetry originates from an inhomogeneous distribution of grain sizes in the column and column overload. However, it is an experimental fact that the distribution of chemicals within an injected volume of solute changes as a function of time, while the response is measured simultaneously. Accordingly, the obtained signal cannot be compared to the theory before some type of time‐based deconvolution of the data has been performed. Adjustments to high‐performance liquid chromatography data were thus proposed through empirical equations that describe the relevant time values, peak height, peak area, and parameters of the van Deemter equation. It was proposed that the transfer of solute from the front to the rear part of the pulse during laminar open‐ended flow occurs at rate that depends on the linear flow rate, and to a lesser extent, on properties of the response function.     

Inflammation Due to Voriconazole‐induced Photosensitivity Enhanced Skin Phototumorigenesis in Xpa‐knockout Mice

Voriconazole is an antifungal agent and used as a prophylactic measure, especially in immunocompromised patients. However, there have been several reports of its adverse reactions, namely photosensitivity with intense inflammatory rashes and subsequent skin cancer development. To assess the effects of photosensitizing drugs voriconazole and hydrochlorothiazide (HCTZ ) on the enhancement of UV ‐induced inflammatory responses and UV ‐induced tumorigenesis, we utilized Xpa ‐knockout mice, which is DNA repair‐deficient and more susceptible to UV ‐induced inflammation and tumor development than wild‐type mice. Administration of voriconazole prior to broadband UVB exposure significantly upregulated multiple inflammatory cytokines compared with the vehicle‐ or HCTZ ‐administered groups. Voriconazole administration along with chronic UVB exposure produced significantly higher number of skin tumors than HCTZ or vehicle in Xpa ‐knockout mice. Furthermore, the investigation of UVB ‐induced DNA damage using embryonic fibroblasts of Xpa ‐knockout mice revealed a significantly higher 8‐oxo‐7,8‐dihydroguanine level in cells treated with voriconazole N‐oxide, a voriconazole‐metabolite during UV exposure. The data suggest that voriconazole plus UVB ‐induced inflammatory response may be related to voriconazole‐induced skin phototumorigenesis.