Tethered Monte Carlo: computing the effective potential without critical slowing down

We present Tethered Monte Carlo, a simple, general purpose method of computing the effective potential of the order parameter (Helmholtz free energy). This formalism is based on a new statistical ensemble, closely related to the micromagnetic one, but with an extended configuration space (through Creutz-like demons). Canonical averages for arbitrary values of the external magnetic field are computed without additional simulations. The method is put to work in the two-dimensional Ising model, where the existence of exact results enables us to perform high precision checks. A rather peculiar feature of our implementation, which employs a local Metropolis algorithm, is the total absence, within errors, of critical slowing down for magnetic observables. Indeed, high accuracy results are presented for lattices as large as L = 1024.

All one-loop scalar vertices in the effective potential approach

Using the one-loop Coleman-Weinberg effective potential, we derive a general analytic expression for all the derivatives of the effective potential with respect to any number of classical scalar fields. The result is valid for a renormalisable theory in four dimensions with any number of scalars, fermions or gauge bosons. This result corresponds to the zero-external momentum contribution to a general one-loop diagram with N scalar external legs. We illustrate the use of the general result in two simple scalar singlet extensions of the Standard Model, to obtain the dominant contributions to the triple couplings of light scalar particles under the zero external momentum approximation.

Suspensions of polymer-grafted nanoparticles with added polymers-Structure and effective pair-interactions

We present the results of combined experimental and theoretical (molecular dynamics simulations and integral equation theory) studies of the structure and effective interactions of suspensions of polymer grafted nanoparticles (PGNPs) in the presence of linear polymers. Due to the absence of systematic experimental and theoretical studies of PGNPs, it is widely believed that the structure and effective interactions in such binary mixtures would be very similar to those of an analogous soft colloidal material-star polymers. In our study, polystyrene-grafted gold nanoparticles with functionality f = 70 were mixed with linear polystyrene (PS) of two different molecular weights for obtaining two PGNP: PS size ratios, xi = 0.14 and 2.76 (where, xi = M-g/M-m, M-g and M-m being the molecular weights of grafting and matrix polymers, respectively). The experimental structure factor of PGNPs could be modeled with an effective potential (Model-X), which has been found to be widely applicable for star polymers. Similarly, the structure factor of the blends with xi = 0.14 could be modeled reasonably well, while the structure of blends with xi = 2.76 could not be captured, especially for high density of added polymers. A model (Model-Y) for effective interactions between PGNPs in a melt of matrix polymers also failed to provide good agreement with the experimental data for samples with xi = 2.76 and high density of added polymers. We tentatively attribute this anomaly in modeling the structure factor of blends with xi = 2.76 to the questionable assumption of Model-X in describing the added polymers as star polymers with functionality 2, which gets manifested in both polymer-polymer and polymer-PGNP interactions especially at higher fractions of added polymers. The failure of Model-Y may be due to the neglect of possible many-body interactions among PGNPs mediated by matrix polymers when the fraction of added polymers is high. These observations point to the need for a new framework to understand not only the structural behavior of PGNPs but also possibly their dynamics and thermo-mechanical properties as well. (C) 2015 AIP Publishing LLC.

Exchange-correlation energy and potential as approximate functionals of occupied and virtual Kohn-Sham orbitals: Application to dissociating H-2

The standard local density approximation and generalized gradient approximations fail to properly describe the dissociation of an electron pair bond, yielding large errors (on the order of 50 kcal/mol) at long bond distances. To remedy this failure, a self-consistent Kohn-Sham (KS) method is proposed with the exchange-correlation (xc) energy and potential depending on both occupied and virtual KS orbitals. The xc energy functional of Buijse and Baerends [Mol. Phys. 100, 401 (2002); Phys. Rev. Lett. 87, 133004 (2001)] is employed, which, based on an ansatz for the xc-hole amplitude, is able to reproduce the important dynamical and nondynamical effects of Coulomb correlation through the efficient use of virtual orbitals. Self-consistent calculations require the corresponding xc potential to be obtained, to which end the optimized effective potential (OEP) method is used within the common energy denominator approximation for the static orbital Green's function. The problem of the asymptotic divergence of the xc potential of the OEP when a finite number of virtual orbitals is used is addressed. The self-consistent calculations reproduce very well the entire H-2 potential curve, describing correctly the gradual buildup of strong left-right correlation in stretched H-2. (C) 2003 American Institute of Physics.

Studies on gaussian effective potentials in some models

In classical field theory, the ordinary potential V is an energy density for that state in which the field assumes the value ¢. In quantum field theory, the effective potential is the expectation value of the energy density for which the expectation value of the field is ¢o. As a result, if V has several local minima, it is only the absolute minimum that corresponds to the true ground state of the theory. Perturbation theory remains to this day the main analytical tool in the study of Quantum Field Theory. However, since perturbation theory is unable to uncover the whole rich structure of Quantum Field Theory, it is desirable to have some method which, on one hand, must go beyond both perturbation theory and classical approximation in the points where these fail, and at that time, be sufficiently simple that analytical calculations could be performed in its framework During the last decade a nonperturbative variational method called Gaussian effective potential, has been discussed widely together with several applications. This concept was described as a means of formalizing our intuitive understanding of zero-point fluctuation effects in quantum mechanics in a way that carries over directly to field theory.

Bounded solutions of neutral fermions with a screened Coulomb potential

The intrinsically relativistic problem of a fermion subject to a pseudoscalar screened Coulomb plus a uniform background potential in two-dimensional space-time is mapped into a Sturm-Liouville. This mapping gives rise to an effective Morse-like potential and exact bounded solutions are found. It is shown that the uniform background potential determinates the number of bound-state solutions. The behaviour of the eigenenergies as well as of the upper and lower components of the Dirac spinor corresponding to bounded solutions is discussed in detail and some unusual results are revealed. An apparent paradox concerning the uncertainty principle is solved by recurring to the concepts of effective mass and effective Compton wavelength. (c) 2005 Elsevier B.V. All rights reserved.

Bound states of the Dirac equation for a class of effective quadratic plus inversely quadratic potentials

The Dirac equation is exactly solved for a pseudoscalar linear plus Coulomb-like potential in a two-dimensional world. This sort of potential gives rise to an effective quadratic plus inversely quadratic potential in a Sturm-Liouville problem, regardless the sign of the parameter of the linear potential, in sharp contrast with the Schrodinger case. The generalized Dirac oscillator already analyzed in a previous work is obtained as a particular case. (C) 2004 Elsevier B.V. All rights reserved.

Supersymmetry, variational method and the Lennard-Jones (12,6) potential

The formalism of supersymmetric quantum mechanics is used to determine trial functions in order to obtain eigenvalues for the Lennard-Jones (12, 6) potential from variational method. The superpotential obtained provides an effective potential which can be directly comparable to the original one.

Quadratic effective action for QED in D=2,3 dimensions

We calculate the effective action for quantum electrodynamics (QED) in D=2,3 dimensions at the quadratic approximation in the gauge fields. We analyze the analytic structure of the corresponding nonlocal boson propagators nonperturbatively in k/m. In two dimensions for any nonzero fermion mass, we end up with one massless pole for the gauge boson. We also calculate in D=2 the effective potential between two static charges separated by a distance L and find it to be a linearly increasing function of L in agreement with the bosonized theory (massive sine-Gordon model). In three dimensions we find nonperturbatively in k/m one massive pole in the effective bosonic action leading to screening. Fitting the numerical results we derive a simple expression for the functional dependence of the boson mass upon the dimensionless parameter e2/m. ©2000 The American Physical Society.

Molecular Origin of the Intrinsic Bending Force for Helical Morphology Observed in Chiral Amphiphilic Assemblies: Concentration and Size Dependence

The formation of the helical morphology in monolayers and bilayers of chiral amphiphilic assemblies is believed to be driven at least partly by the interactions at the chiral centers of the amphiphiles. However, a detailed microscopic understanding of these interactions and their relation with the helix formation is still not clear. In this article a study of the molecular origin of the chirality-driven helix formation is presented by calculating, for the first time, the effective pair potential between a pair of chiral molecules. This effective potential depends on the relative sizes of the groups attached to the two chiral centers, on the orientation of the amphiphile molecules, and also on the distance between them. We find that for the mirror-image isomers (in the racemic modification) the minimum energy conformation is a nearly parallel alignment of the molecules. On the other hand, the same for a pair of molecules of one kind of enantiomer favors a tilt angle between them, thus leading to the formation of a helical morphology of the aggregate. The tilt angle is determined by the size of the groups attached to the chiral centers of the pair of molecules considered and in many cases predicted it to be close to 45 degrees. The present study, therefore, provides a molecular origin of the intrinsic bending force, suggested by Helfrich (J. Chem. Phys. 1986, 85, 1085-1087), to be responsible for the formation of helical structure. This effective potential may explain many of the existing experimental results, such as the size and the concentration dependence of the formation of helical morphology. It is further found that the elastic forces can significantly modify the pitch predicted by the chiral interactions alone and that the modified real pitch is close to the experimentally observed value. The present study is expected to provide a starting point for future microscopic studies.

Phase transitions and gravity

Einstein's gravitational field is non-minimally coupled to a self-interacting scalar field in the presence of radiation. Such a theory can give rise to a phase transition associated with a change of sign of the gravitational “constant”. In our approach, the criterion for stability is formulated in terms of an effective potential, the phase-transition takes place due to temperature dependence of the scalar self-interaction coupling constant.

Microscopic origin of the chirality driven morphologies of the amphiphilic monolayers and bilayers

It is widely known that the compressed monolayers and bilayers of chiral lipids or fatty acids form helical morphologies, while the corresponding racemic modification gives only flat platelets without twist. No molecular explanation of this phenomenon is yet available, although subtle interactions at the chiral centers have often been proposed as the driving force behind the morphology of the aggregate to form a particular shape. In the present study, the morphologies of the chiral amphiphilic assemblies have been predicted on the basis of an effective pair potential between the molecules, which depends on the relative sizes of the groups attached to the chiral centers, the orientation of the amphiphilic molecules and also on the distance between them. It is shown that fur a pair of same kind of enantiomers, the minimum energy conformation favours a twist angle between them. This twist between the neighbouring molecules gives rise to the helicity of the aggregate. The present theory also shows from the molecular considerations that for a pair of mirror-image isomers (i.e. the racemic modification) the minimum energy conformation corresponds to the zero angle between the molecules, thus giving rise to flat platelets as observed in experiments. Another fascinating aspect of such chirality driven helical structures is that the sense (or the handedness) of the helix is highly specific about the chirality of the monomer concerned. The molecular theory shows, for the first time, that the sense of the helical structures in many cases is determined by the sizes of the groups attached to the chiral centers and the effective potential between them. The predicted senses of the helical structures are in complete agreement with the experimental results.

Globally coupled stochastic two-state oscillators: Fluctuations due to finite numbers

Infinite arrays of coupled two-state stochastic oscillators exhibit well-defined steady states. We study the fluctuations that occur when the number N of oscillators in the array is finite. We choose a particular form of global coupling that in the infinite array leads to a pitchfork bifurcation from a monostable to a bistable steady state, the latter with two equally probable stationary states. The control parameter for this bifurcation is the coupling strength. In finite arrays these states become metastable: The fluctuations lead to distributions around the most probable states, with one maximum in the monostable regime and two maxima in the bistable regime. In the latter regime, the fluctuations lead to transitions between the two peak regions of the distribution. Also, we find that the fluctuations break the symmetry in the bimodal regime, that is, one metastable state becomes more probable than the other, increasingly so with increasing array size. To arrive at these results, we start from microscopic dynamical evolution equations from which we derive a Langevin equation that exhibits an interesting multiplicative noise structure. We also present a master equation description of the dynamics. Both of these equations lead to the same Fokker-Planck equation, the master equation via a 1/N expansion and the Langevin equation via standard methods of Ito calculus for multiplicative noise. From the Fokker-Planck equation we obtain an effective potential that reflects the transition from the monomodal to the bimodal distribution as a function of a control parameter. We present a variety of numerical and analytic results that illustrate the strong effects of the fluctuations. We also show that the limits N -> infinity and t -> infinity(t is the time) do not commute. In fact, the two orders of implementation lead to drastically different results.

Theoretical Insights to Niobium-Doped Monolayer MoS2-Gold Contact

We report a first principles study of the electronic properties for a contact formed between Nb-doped monolayer MoS2 and gold for different doping concentrations. We first focus on the shift of energy levels in band structure and the density of states with respect to the Fermi level for a geometrically optimized 5 x 5 MoS2 supercell for both pristine and Nb-doped structures. The doping is achieved by substituting Mo atoms with Nb atoms at random positions. It is observed that for an experimentally reported sheet hole doping concentration of (rho(2D)) 1.8 x 10(14) cm(-2), the pristine MoS2 converts to degenerate p-type semiconductor. Next, we interface this supercell with six layers of < 111 > cleaved surface of gold to investigate the contact nature of MoS2-Au system. By careful examination of projected band structure, projected density of states, effective potential and charge density difference, we demonstrate that the Schottky barrier nature observed for pure MoS2-Au contact can be converted from n-type to p-type by efficient Nb doping.

Lead selection and characterization of antitubercular compounds using the Nested Chemical Library

Discovering new drugs to treat tuberculosis more efficiently and to overcome multidrug resistance is a world health priority. To find novel antitubercular agents several approaches have been used in various institutions worldwide, including target-based approaches against several validated mycobacterial enzymes and phenotypic screens. We screened more than 17,000 compounds from Vichem's Nested Chemical Library(TM) using an integrated strategy involving whole cell-based assays with Corynebacterium glutamicum and Mycobacterium tuberculosis, and target-based assays with protein kinases PknA, PknB and PknG as well as other targets such as PimA and bacterial topoisomerases simultaneously. With the help of the target-based approach we have found very potent hits inhibiting the selected target enzymes, but good minimal inhibitory concentrations (MIC) against M. tuberculosis were not achieved. Focussing on the whole cell-based approach several potent hits were found which displayed minimal inhibitory concentrations (MIC) against M. tuberculosis below 10 mu M and were non-mutagenic, non-cytotoxic and the targets of some of the hits were also identified. The most active hits represented various scaffolds. Medicinal chemistry-based lead optimization was performed applying various strategies and, as a consequence, a series of novel potent compounds were synthesized. These efforts resulted in some effective potential antitubercular lead compounds which were confirmed in phenotypic assays. (C) 2015 Elsevier Ltd. All rights reserved.