Performance analysis of linear and nonlinear techniques for automatic scaling of discretized control problems

This paper introduces three (one linear and two nonlinear) automatic scaling techniques for NLPs with states and constraints spread over several orders of magnitude, without requiring complex off-the-shelf external tools. All of these methods have been compared to standard techniques and applied to three problems using SNOPT and IPOPT. The results confirm that the proposed techniques significantly improve the NLP conditioning, yielding more reliable and in some cases, faster NLP solutions.     

Toward a universal quantum QSPR operator

This article presents a discussion about the formalism, which might be associated to a general Quantum quantitative structure–properties relations operator, appearing in a Boltzmann‐like exponential form, which is based in turn on the definition of the concept of thermal voltage, applied to thermally scaled electronic density functions. Three practical numerical examples are presented, corresponding to the calculation of the polarization angle in assorted chiral molecules, the estimation of fish toxicity for perchlorobenzene within the set of chlorobenzenes and a typical quantum QSAR study on the popular Cramer steroid set.     

Nondimensional frequency scaling of aerodynamically-tensioned membranes

Membrane wings have applications that involve low Reynolds number flyers such as micro air vehicles. The time-averaged and time-dependent deformations of the membrane affect the aerodynamic characteristics of the wing, primarily in the region beyond the maximum aerodynamic efficiency of the wing. This paper investigates an appropriate nondimensional vibration frequency scaling of a spanwise tensioned membrane with free (unattached) leading and trailing edges at low Reynolds numbers relative to nondimensional aeroelastic parameters. Silicone rubber membranes with varying spanwise pre-tension, aerodynamic tension (due to wing angle-of-attack and flow dynamic pressure), modulus of elasticity, span, and thickness are studied. Experimental results are compared to a proposed scaling that simplifies the aerodynamic loading as a uniform pressure distribution acting on the membrane. Data is further compared and discussed relative to previous published results of membrane wings with finite wing spans (three-dimensional flow) and fixed (rigid) leading edges.     

Single determinant N‐representability and the kernel energy method applied to water clusters

The Kernel energy method (KEM) is a quantum chemical calculation method that has been shown to provide accurate energies for large molecules. KEM performs calculations on subsets of a molecule (called kernels) and so the computational difficulty of KEM calculations scales more softly than full molecule methods. Although KEM provides accurate energies those energies are not required to satisfy the variational theorem. In this article, KEM is extended to provide a full molecule single‐determinant N‐representable one‐body density matrix. A kernel expansion for the one‐body density matrix analogous to the kernel expansion for energy is defined. This matrix is converted to a normalized projector by an algorithm due to Clinton. The resulting single‐determinant N‐representable density matrix maps to a quantum mechanically valid wavefunction which satisfies the variational theorem. The process is demonstrated on clusters of three to twenty water molecules. The resulting energies are more accurate than the straightforward KEM energy results and all violations of the variational theorem are resolved. The N‐representability studied in this article is applicable to the study of quantum crystallography. © 2017 Wiley Periodicals, Inc.     

A class of optimal iterative methods of inverting a linear bounded operator

A family of higher-order iterative methods for inverting a linear bounded operator in a Banach space is presented. This family depending on two integer parameters can be considered as hyperpower methods of [1]. A convergence Theorem is proved and an optimal class with respect to the efficiency index is determined . Several error bounds are derived and compared with those in the literature.     

A scaling limit for the length of the longest cycle in a sparse random digraph

We discuss the length of the longest directed cycle in the sparse random digraph , constant. We show that for large there exists a function such that a.s. The function where is a polynomial in . We are only able to explicitly give the values , although we could in principle compute any .     

Synthesis and Liquid Crystalline Behavior of Random Copolymer of Poly(ethylene oxide) Macromonomer and Liquid Crystalline Monomer by the Photon Transmission Technique

Random copolymers of poly(ethylene oxide) macromonomer with p‐vinylbenzyl end‐functional group (PEOVB) and liquid crystalline monomer, namely 6‐(4‐cyanobiphenyl‐4′‐oxy)hexyl acrylate (COA), were prepared by conventional free radical polymerization. A living anionic polymerization technique was employed for the synthesis of PEO macromonomers bearing p‐vinylbenzyl moiety at one end. The photon transmission method was also applied to study the phase transitions of COA monomer and its random copolymer with PEO. It was found that, for both samples, the nematic‐smectic A transition is continuous, but the critical fluctuation regions do not allow to obtain 3D XY values. Instead, we have obtained the values close to mean field regime. Scaling of thermal hystersis for random copolymer sample near the nematic‐isotropic transition was studied as well. Thermal hysteresis loops were produced under linearly varying temperature. It was shown that the areas of the hysteresis loops scale with the temperature scanning rate with an exponent being equal to 0.614 which is in good agreement with the field‐theoretical value.     

Theory of variational calculation with a scaling correct moment functional to solve the electronic schrödinger equation directly for ground state one‐electron density and electronic energy

The reduction of the electronic Schrodinger equation or its calculating algorithm from 4N‐dimensions to a nonlinear, approximate density functional of a three spatial dimension one‐electron density for an N electron system which is tractable in practice, is a long‐desired goal in electronic structure calculation. In a seminal work, Parr et al. (Phys. Rev. A 1997, 55, 1792) suggested a well behaving density functional in power series with respect to density scaling within the orbital‐free framework for kinetic and repulsion energy of electrons. The updated literature on this subject is listed, reviewed, and summarized. Using this series with some modifications, a good density functional approximation is analyzed and solved via the Lagrange multiplier device. (We call the attention that the introduction of a Lagrangian multiplier to ensure normalization is a new element in this part of the related, general theory.) Its relation to Hartree–Fock (HF) and Kohn–Sham (KS) formalism is also analyzed for the goal to replace all the analytical Gaussian based two and four center integrals (∫g_i( r ₁)g_k( r ₂)rd r ₁d r ₂, etc.) to estimate electron‐electron interactions with cheaper numerical integration. The KS method needs the numerical integration anyway for correlation estimation. © 2012 Wiley Periodicals, Inc.     

Natural bond orbital analysis in the ONETEP code: Applications to large protein systems

First principles electronic structure calculations are typically performed in terms of molecular orbitals (or bands), providing a straightforward theoretical avenue for approximations of increasing sophistication, but do not usually provide any qualitative chemical information about the system. We can derive such information via post‐processing using natural bond orbital (NBO) analysis, which produces a chemical picture of bonding in terms of localized Lewis‐type bond and lone pair orbitals that we can use to understand molecular structure and interactions. We present NBO analysis of large‐scale calculations with the ONETEP linear‐scaling density functional theory package, which we have interfaced with the NBO 5 analysis program. In ONETEP calculations involving thousands of atoms, one is typically interested in particular regions of a nanosystem whilst accounting for long‐range electronic effects from the entire system. We show that by transforming the Non‐orthogonal Generalized Wannier Functions of ONETEP to natural atomic orbitals, NBO analysis can be performed within a localized region in such a way that ensures the results are identical to an analysis on the full system. We demonstrate the capabilities of this approach by performing illustrative studies of large proteins—namely, investigating changes in charge transfer between the heme group of myoglobin and its ligands with increasing system size and between a protein and its explicit solvent, estimating the contribution of electronic delocalization to the stabilization of hydrogen bonds in the binding pocket of a drug‐receptor complex, and observing, in situ, the n → π* hyperconjugative interactions between carbonyl groups that stabilize protein backbones. © 2012 Wiley Periodicals, Inc.