Simulation of inhomogeneous distributions of ultracold atoms in an optical lattice via a massively parallel implementation of nonequilibrium strong-coupling perturbation theory

We present a nonequilibrium strong-coupling approach to inhomogeneous systems of ultracold atoms in optical lattices. We demonstrate its application to the Mott-insulating phase of a two-dimensional Fermi-Hubbard model in the presence of a trap potential. Since the theory is formulated self-consistently, the numerical implementation relies on a massively parallel evaluation of the self-energy and the Green's function at each lattice site, employing thousands of CPUs. While the computation of the self-energy is straightforward to parallelize, the evaluation of the Green's function requires the inversion of a large sparse 10(d) x 10(d) matrix, with d > 6. As a crucial ingredient, our solution heavily relies on the smallness of the hopping as compared to the interaction strength and yields a widely scalable realization of a rapidly converging iterative algorithm which evaluates all elements of the Green's function. Results are validated by comparing with the homogeneous case via the local-density approximation. These calculations also show that the local-density approximation is valid in nonequilibrium setups without mass transport.

Simultaneous Perturbation Newton Algorithms for Simulation Optimization

We present a new Hessian estimator based on the simultaneous perturbation procedure, that requires three system simulations regardless of the parameter dimension. We then present two Newton-based simulation optimization algorithms that incorporate this Hessian estimator. The two algorithms differ primarily in the manner in which the Hessian estimate is used. Both our algorithms do not compute the inverse Hessian explicitly, thereby saving on computational effort. While our first algorithm directly obtains the product of the inverse Hessian with the gradient of the objective, our second algorithm makes use of the Sherman-Morrison matrix inversion lemma to recursively estimate the inverse Hessian. We provide proofs of convergence for both our algorithms. Next, we consider an interesting application of our algorithms on a problem of road traffic control. Our algorithms are seen to exhibit better performance than two Newton algorithms from a recent prior work.

Simultaneous perturbation methods for adaptive labor staffing in service systems

We consider the problem of optimizing the workforce of a service system. Adapting the staffing levels in such systems is non-trivial due to large variations in workload and the large number of system parameters do not allow for a brute force search. Further, because these parameters change on a weekly basis, the optimization should not take longer than a few hours. Our aim is to find the optimum staffing levels from a discrete high-dimensional parameter set, that minimizes the long run average of the single-stage cost function, while adhering to the constraints relating to queue stability and service-level agreement (SLA) compliance. The single-stage cost function balances the conflicting objectives of utilizing workers better and attaining the target SLAs. We formulate this problem as a constrained parameterized Markov cost process parameterized by the (discrete) staffing levels. We propose novel simultaneous perturbation stochastic approximation (SPSA)-based algorithms for solving the above problem. The algorithms include both first-order as well as second-order methods and incorporate SPSA-based gradient/Hessian estimates for primal descent, while performing dual ascent for the Lagrange multipliers. Both algorithms are online and update the staffing levels in an incremental fashion. Further, they involve a certain generalized smooth projection operator, which is essential to project the continuous-valued worker parameter tuned by our algorithms onto the discrete set. The smoothness is necessary to ensure that the underlying transition dynamics of the constrained Markov cost process is itself smooth (as a function of the continuous-valued parameter): a critical requirement to prove the convergence of both algorithms. We validate our algorithms via performance simulations based on data from five real-life service systems. For the sake of comparison, we also implement a scatter search based algorithm using state-of-the-art optimization tool-kit OptQuest. From the experiments, we observe that both our algorithms converge empirically and consistently outperform OptQuest in most of the settings considered. This finding coupled with the computational advantage of our algorithms make them amenable for adaptive labor staffing in real-life service systems.

Interaction-induced singular Fermi surface in a high-temperature oxypnictide superconductor

In the family of iron-based superconductors, LaFeAsO-type materials possess the simplest electronic structure due to their pronounced two-dimensionality. And yet they host superconductivity with the highest transition temperature T-c approximate to 55K. Early theoretical predictions of their electronic structure revealed multiple large circular portions of the Fermi surface with a very good geometrical overlap (nesting), believed to enhance the pairing interaction and thus superconductivity. The prevalence of such large circular features in the Fermi surface has since been associated with many other iron-based compounds and has grown to be generally accepted in the field. In this work we show that a prototypical compound of the 1111-type, SmFe0.92Co0.08AsO, is at odds with this description and possesses a distinctly different Fermi surface, which consists of two singular constructs formed by the edges of several bands, pulled to the Fermi level from the depths of the theoretically predicted band structure by strong electronic interactions. Such singularities dramatically affect the low-energy electronic properties of the material, including superconductivity. We further argue that occurrence of these singularities correlates with the maximum superconducting transition temperature attainable in each material class over the entire family of iron-based superconductors.

Conformal perturbation theory and higher spin entanglement entropy on the torus

We study the free fermion theory in 1+1 dimensions deformed by chemical potentials for holomorphic, conserved currents at finite temperature and on a spatial circle. For a spin-three chemical potential mu, the deformation is related at high temperatures to a higher spin black hole in hs0] theory on AdS(3) spacetime. We calculate the order mu(2) corrections to the single interval Renyi and entanglement entropies on the torus using the bosonized formulation. A consistent result, satisfying all checks, emerges upon carefully accounting for both perturbative and winding mode contributions in the bosonized language. The order mu(2) corrections involve integrals that are finite but potentially sensitive to contact term singularities. We propose and apply a prescription for defining such integrals which matches the Hamiltonian picture and passes several non-trivial checks for both thermal corrections and the Renyi entropies at this order. The thermal corrections are given by a weight six quasi-modular form, whilst the Renyi entropies are controlled by quasi-elliptic functions of the interval length with modular weight six. We also point out the well known connection between the perturbative expansion of the partition function in powers of the spin-three chemical potential and the Gross-Taylor genus expansion of large-N Yang-Mills theory on the torus. We note the absence of winding mode contributions in this connection, which suggests qualitatively different entanglement entropies for the two systems.

Impact of glycosylation on stability, structure and unfolding of soybean agglutinin (SBA): an insight from thermal perturbation molecular dynamics simulations

Glycosylation has been recognized as one of the most prevalent and complex post-translational modifications of proteins involving numerous enzymes and substrates. Its effect on the protein conformational transitions is not clearly understood yet. In this study, we have examined the effect of glycosylation on protein stability using molecular dynamics simulation of legume lectin soybean agglutinin (SBA). Its glycosylated moiety consists of high mannose type N-linked glycan (Man(9)GlcNAc(2)). To unveil the structural perturbations during thermal unfolding of these two forms, we have studied and compared them to the experimental results. From the perspective of dynamics, our simulations revealed that the nonglycosylated monomeric form is less stable than corresponding glycosylated form at normal and elevated temperatures. Moreover, at elevated temperature thermal destabilization is more prominent in solvent exposed loops, turns and ends of distinct beta sheets. SBA maintains it folded structure due to some important saltbridges, hydrogen bonds and hydrophobic interactions within the protein. The reducing terminal GlcNAc residues interact with the protein residues VAL161, PRO182 and SER225 via hydrophobic and via hydrogen bonding with ASN 9 and ASN 75. Our simulations also revealed that single glycosylation (ASN75) has no significant effect on corresponding cis peptide angle orientation. This atomistic description might have important implications for understanding the functionality and stability of Soybean agglutinin.

Ultra-sensitive pressure dependence of bandgap of rutile-GeO2 revealed by many body perturbation theory

The reported values of bandgap of rutile GeO2 calculated by the standard density functional theory within local-density approximation (LDA)/generalized gradient approximation (GGA) show a wide variation (similar to 2 eV), whose origin remains unresolved. Here, we investigate the reasons for this variation by studying the electronic structure of rutile-GeO2 using many-body perturbation theory within the GW framework. The bandgap as well as valence bandwidth at Gamma-point of rutile phase shows a strong dependence on volume change, which is independent of bandgap underestimation problem of LDA/GGA. This strong dependence originates from a change in hybridization among O-p and Ge-(s and p) orbitals. Furthermore, the parabolic nature of first conduction band along X-Gamma-M direction changes towards a linear dispersion with volume expansion. (C) 2015 AIP Publishing LLC.

Perturbation of nucleo-cytoplasmic transport affects size of nucleus and nucleolus in human cells

Size regulation of human cell nucleus and nucleolus are poorly understood subjects. 3D reconstruction of live image shows that the karyoplasmic ratio (KR) increases by 30-80% in transformed cell lines compared to their immortalized counterpart. The attenuation of nucleo-cytoplasmic transport causes the KR value to increase by 30-50% in immortalized cell lines. Nucleolus volumes are significantly increased in transformed cell lines and the attenuation of nucleo-cytoplasmic transport causes a significant increase in the nucleolus volume of immortalized cell lines. A cytosol and nuclear fraction swapping experiment emphasizes the potential role of unknown cytosolic factors in nuclear and nucleolar size regulation.

Moduli spaces of vector bundles on a singular rational ruled surface

We study moduli spaces M-X (r, c(1), c(2)) parametrizing slope semistable vector bundles of rank r and fixed Chern classes c(1), c(2) on a ruled surface whose base is a rational nodal curve. We showthat under certain conditions, these moduli spaces are irreducible, smooth and rational (when non-empty). We also prove that they are non-empty in some cases. We show that for a rational ruled surface defined over real numbers, the moduli space M-X (r, c(1), c(2)) is rational as a variety defined over R.

Complex potentials and singular integral equation for curve crack problem in antiplane elasticity

Four types of the fundamental complex potential in antiplane elasticity are introduced: (a) a point dislocation, (b) a concentrated force, (c) a dislocation doublet and (d) a concentrated force doublet. It is proven that if the axis of the concentrated force doublet is perpendicular to the direction of the dislocation doublet, the relevant complex potentials are equivalent. Using the obtained complex potentials, a singular integral equation for the curve crack problem is introduced. Some particular features of the obtained singular integral equation are discussed, and numerical solutions and examples are given.

Perturbation of the stability of amyloid fibrils through alteration of electrostatic interactions.

The self-assembly of proteins and peptides into polymeric amyloid fibrils is a process that has important implications ranging from the understanding of protein misfolding disorders to the discovery of novel nanobiomaterials. In this study, we probe the stability of fibrils prepared at pH 2.0 and composed of the protein insulin by manipulating electrostatic interactions within the fibril architecture. We demonstrate that strong electrostatic repulsion is sufficient to disrupt the hydrogen-bonded, cross-β network that links insulin molecules and ultimately results in fibril dissociation. The extent of this dissociation correlates well with predictions for colloidal models considering the net global charge of the polypeptide chain, although the kinetics of the process is regulated by the charge state of a single amino acid. We found the fibrils to be maximally stable under their formation conditions. Partial disruption of the cross-β network under conditions where the fibrils remain intact leads to a reduction in their stability. Together, these results support the contention that a major determinant of amyloid stability stems from the interactions in the structured core, and show how the control of electrostatic interactions can be used to characterize the factors that modulate fibril stability.

Del combate singular al singular combate sexual en la sátira trovadoresca medieval gallego-portuguesa

Integran este número de la revista ponencias presentadas en Studia Hispanica Medievalia VIII: Actas de las IX Jornadas Internacionales de Literatura Española Medieval, 2008, y de Homenaje al Quinto Centenario de Amadis de Gaula.

Dynamic mode I perturbation solution for a moving crack unsteadily

Rice et al. (Jounal of Mechanics and Physics of Solids 42, 813-843) analyze the propagation of a planar crack with a nominally straight front in a model elastic solid with a single displacement component. Using the form of Willis er al. (Journal of the Mechanics and Physics of Solids 43, 319-341), of dynamic mode I weight functions for a moving crack, we address that problem solved by Rice ei al. in the 3D context of elastodynamic theory. Oscillatory crack tip motion results from constructive-destructive interference of stress intensity waves. Those waves, including system of the dilatational, shear and Rayleigh waves, interact on each other and with moving edge of crack, can lead to continuing fluctuations of the crack front and propagation velocity. (C) 1997 Elsevier Science Ltd.

Using wavelet transform to study the Lipschitz local singular exponent in wall turbulence

In this paper, wavelet,transform is introduced to study the Lipschitz local singular exponent for characterising the local singularity behavior of fluctuating velocity in wall turbulence. I, is found that the local singular exponent is negative when the ejections and sweeps of coherent structures occur in a turbulent boundary layer.

Cauchy singular integral equation method for transient antiplane dynamic problems

In this paper, by use of the boundary integral equation method and the techniques of Green basic solution and singularity analysis, the dynamic problem of antiplane is investigated. The problem is reduced to solving a Cauchy singular integral equation in Laplace transform space. This equation is strictly proved to be equivalent to the dual integral equations obtained by Sih [Mechanics of Fracture, Vol. 4. Noordhoff, Leyden (1977)]. On this basis, the dynamic influence between two parallel cracks is also investigated. By use of the high precision numerical method for the singular integral equation and Laplace numerical inversion, the dynamic stress intensity factors of several typical problems are calculated in this paper. The related numerical results are compared to be consistent with those of Sih. It shows that the method of this paper is successful and can be used to solve more complicated problems. Copyright (C) 1996 Elsevier Science Ltd