Engineering Papers⌕ Search

Engineering topics

Alexeev, Yuri

Publications and source records attributed to Alexeev, Yuri.

53 records · Page 3

Learning to Optimize Variational Quantum Circuits to Solve Combinatorial Problems

Quantum computing is a computational paradigm with the potential to outperform classical methods for a variety of problems. Proposed recently, the Quantum Approximate Optimization Algorithm (QAOA) is considered as one of the leading candidates for demonstrating quantum advantage in the near term. QAOA is a variational hybrid quantum-classical algorithm for approximately solving combinatorial optimization problems. The quality of the solution obtained by QAOA for a given problem instance depends on the performance of the classical optimizer used to optimize the variational parameters. In this paper, we formulate the problem of finding optimal QAOA parameters as a learning task in which the knowledge gained from solving training instances can be leveraged to find high-quality solutions for unseen test instances. To this end, we develop two machine-learning-based approaches. Our first approach adopts a reinforcement learning (RL) framework to learn a policy network to optimize QAOA circuits. Our second approach adopts a kernel density estimation (KDE) technique to learn a generative model of optimal QAOA parameters. In both approaches, the training procedure is performed on small-sized problem instances that can be simulated on a classical computer; yet the learned RL policy and the generative model can be used to efficiently solve larger problems. Furthermore, extensive simulations using the IBM Qiskit Aer quantum circuit simulator demonstrate that our proposed RL- and KDE-based approaches reduce the optimality gap by factors up to 30.15 when compared with other commonly used off-the-shelf optimizers.

97 MATHEMATICS AND COMPUTING↗

Quantum Divide and Compute: Hardware Demonstrations and Noisy Simulations

Noisy, intermediate-scale quantum computers come with intrinsic limitations in terms of the number of qubits (circuit "width") and decoherence time (circuit "depth") they can have. Here, for the first time, we demonstrate a recently introduced method that breaks a circuit into smaller subcircuits or fragments, and thus makes it possible to run circuits that are either too wide or too deep for a given quantum processor. We investigate the behavior of the method on one of IBM's 20-qubit superconducting quantum processors with various numbers of qubits and fragments. We build noise models that capture decoherence, readout error, and gate imperfections for this particular processor. We then carry out noisy simulations of the method in order to account for the observed experimental results. We find an agreement within 20% between the experimental and the simulated success probabilities, and we observe that recombining noisy fragments yields overall results that can outperform the results without fragmentation.

Ayral, Thomas↗

Accelerating Quantum Technologies Through Fermilab SRF Expertise

Industry forecast of exponential growth in quantum computing power still predicts that transformative applications are decades away. System performance is limited by decoherence and requires 5-10x the number of physical qubits per effective qubit.

43 PARTICLE ACCELERATORS↗

Locating Minimum Energy Crossings of Different Spin States Using the Fragment Molecular Orbital Method

Spin-dependent processes involving nonadiabatic transitions between electronic states with different spin multiplicities play important roles in the chemistry of complex systems. The rates of these processes can be predicted based on the molecular properties at the minimum energy crossing point (MECP) between electronic states. We present the development of the MECP search technique within the fragment molecular orbital (FMO) method applicable to large complex systems. The accuracy and scalability of the new method is demonstrated on several models of the metal-sulfur protein rubredoxin. The effect of the model size on the MECP geometry and relative energy is discussed. The fragment energy decomposition and spin density delocalization analyses reveal how different protein residues and solvent molecules contribute to stabilization of the spin states. The developed FMO-MECP method can help to clarify the role of nonadiabatic spin-dependent processes in complex systems and can be used for designing mutations aimed at controlling these processes in metalloproteins, including spin-dependent catalysis and electron transfer.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Use cases of lossy compression for floating-point data in scientific data sets

Architectural and technological trends of systems used for scientific computing call for a significant reduction of scientific data sets that are composed mainly of floating-point data. Here, this article surveys and presents experimental results of currently identified use cases of generic lossy compression to address the different limitations of scientific computing systems. The article shows from a collection of experiments run on parallel systems of a leadership facility that lossy data compression not only can reduce the footprint of scientific data sets on storage but also can reduce I/O and checkpoint/restart times, accelerate computation, and even allow significantly larger problems to be run than without lossy compression. In conclusion, these results suggest that lossy compression will become an important technology in many aspects of high performance scientific computing. Because the constraints for each use case are different and often conflicting, this collection of results also indicates the need for more specialization of the compression pipelines.

97 MATHEMATICS AND COMPUTING↗

Molecular Electrostatic Potential and Electron Density of Large Systems in Solution Computed with the Fragment Molecular Orbital Method

A solvent screening model for the molecular electrostatic potential is developed for the fragment molecular orbital method combined with the polarizable continuum model at the Hartree–Fock and density functional theory levels. The accuracy of the generated potentials is established in comparison to calculations without fragmentation. Solvent effects upon the molecular electrostatic potential and electron density of solute are discussed. Finally, the method is applied to two proteins: chignolin (PDB: 1UAO) and ovine prostaglandin H(2) synthase-1 (1EQG).

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Network Community Detection on Small Quantum Computers

In recent years, a number of quantum computing devices with small numbers of qubits have become available. A hybrid quantum local search (QLS) approach that combines a classical machine and a small quantum device to solve problems of practical size is presented. Herein, the proposed approach is applied to the network community detection problem. QLS is hardware-agnostic and easily extendable to new quantum computing devices as they become available. It is demonstrated to solve the 2-community detection problem on graphs of sizes of up to 410 vertices using the 16-qubit IBM quantum computer and D-Wave 2000Q, and compare their performance with the optimal solutions. The results herein demonstrate that QLS performs similarly in terms of quality of the solution and the number of iterations to convergence on both types of quantum computers and it is capable of achieving results comparable to state-of-the-art solvers in terms of quality of the solution including reaching the optimal solutions.

97 MATHEMATICS AND COMPUTING↗

Multithreaded parallelization of the energy and analytic gradient in the fragment molecular orbital method

The fragment molecular orbital method in GAMESS is parallelized in a multithreaded OpenMP implementation combined with the MPI version of the two-level generalized distributed data interface. The energy and analytic gradient in gas phase and the polarizable continuum model of solvation are parallelized in this hybrid three-level scheme, achieving a large memory footprint reduction and a high parallel efficiency on Intel Xeon Phi processors. The parallel efficiency is demonstrated on the Stampede2 and Theta supercomputers using up to 2048 nodes (262 144 threads).

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Real-time data analysis for medical diagnosis using FPGA-accelerated neural networks

Real-time analysis of patient data during medical procedures can provide vital diagnostic feedback that significantly improves chances of success. With sensors becoming increasingly fast, frameworks such as Deep Neural Networks are required to perform calculations within the strict timing constraints for real-time operation. However, traditional computing platforms responsible for running these algorithms incur a large overhead due to communication protocols, memory accesses, and static (often generic) architectures. In this work, we implement a low-latency Multi-Layer Perceptron (MLP) processor using Field Programmable Gate Arrays (FPGAs). Unlike CPUs and Graphics Processing Units (GPUs), our FPGA-based design can directly interface sensors, storage devices, display devices and even actuators, thus reducing the delays of data movement between ports and compute pipelines. Moreover, the compute pipelines themselves are tailored specifically to the application, improving resource utilization and reducing idle cycles. We demonstrate the effectiveness of our approach using mass-spectrometry data sets for real-time cancer detection. We demonstrate that correct parameter sizing, based on the application, can reduce latency by 20% on average. Furthermore, we show that in an application with tightly coupled data-path and latency constraints, having a large amount of computing resources can actually reduce performance. Using mass-spectrometry benchmarks, we show that our proposed FPGA design outperforms both CPU and GPU implementations, with an average speedup of 144x and 21x, respectively. In our work, we demonstrate the importance of application-specific optimizations in order to minimize latency and maximize resource utilization for MLP inference. By directly interfacing and processing sensor data with ultra-low latency, FPGAs can perform real-time analysis during procedures and provide diagnostic feedback that can be critical to achieving higher percentages of successful patient outcomes.

60 APPLIED LIFE SCIENCES↗

A systematic study of minima in alanine dipeptide

The alanine dipeptide is a standard system to model dihedral angles in proteins. It is shown that obtaining the Ramachandran plot accurately is a hard problem because of many local minima; depending on the details of geometry optimizations, different Ramachandran plots can be obtained. To locate all energy minima, starting from geometries from MD simulations, 250,000 geometry optimizations were performed at the level of RHF/6-31G*, followed by re-optimizations of the located 827 minima at the level of MP2/6-311++G**, yielding 30 unique minima, most of which were not previously reported in literature. Both in vacuo and solvated structures are discussed. The minima are systematically categorized based on four backbone dihedral angles. The Gibbs energies are evaluated and the structural factors determining the relative stabilities of conformers are discussed. (c) 2018 Wiley Periodicals, Inc.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Quantum Sensing for High Energy Physics

Report of the first workshop to identify approaches and techniques in the domain of quantum sensing that can be utilized by future High Energy Physics applications to further the scientific goals of High Energy Physics.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

An efficient MPI/OpenMP parallelization of the Hartree–Fock–Roothaan method for the first generation of Intel® Xeon Phi™ processor architecture

The Hartree-Fock (HF) method in the quantum chemistry package GAMESS represents one of the most irregular algorithms in computation today. Major steps in the calculation are the irregular computation of electron repulsion integrals (ERIs) and the building of the Fock matrix. These are the central components of the main Self Consistent Field (SCF) loop, the key hotspot in Electronic Structure (ES) codes. By threading the MPI ranks in the official release of the GAMESS code, we not only speed up the main SCF loop (4x to 6x for large systems), but also achieve a significant (>2x) reduction in the overall memory footprint. These improvements are a direct consequence of memory access optimizations within the MPI ranks. We benchmark our implementation against the official release of the GAMESS code on the Intel R Xeon PhiTM supercomputer. Here, scaling numbers are reported on up to 7,680 cores on Intel Xeon Phi coprocessors.

97 MATHEMATICS AND COMPUTING↗

Advanced Electronic Structure Methods for Heterogeneous Catalysis and Separation of Heavy Metals

For code optimizations on Theta, graphene sheet systems were used for benchmarking purposes. Graphene sheets are of interest to researchers in the area of microlubricants. The proposed scientific applications for this project were not achievable with the support awarded for this project. The performance improvements on Theta reported in this work will directly support the scientific applications described in our recently submitted 2018 INCITE proposal.

42 ENGINEERING↗

Ligand binding to an allergenic lipid transfer protein enhances conformational flexibility resulting in an increase in susceptibility to gastroduodenal proteolysis

Non-specific lipid transfer proteins (LTPs) are a family of lipid-binding molecules that are widely distributed across flowering plant species, many of which have been identified as allergens. They are highly resistant to simulated gastroduodenal proteolysis, a property that may play a role in determining their allergenicity and it has been suggested that lipid binding may further increase stability to proteolysis. It is demonstrated that LTPs from wheat and peach bind a range of lipids in a variety of conditions, including those found in the gastroduodenal tract. Both LTPs are initially cleaved during gastroduodenal proteolysis at three major sites between residues 39–40, 56–57 and 79–80, with wheat LTP being more resistant to cleavage than its peach ortholog. The susceptibility of wheat LTP to proteolyic cleavage increases significantly upon lipid binding. This enhanced digestibility is likely to be due to the displacement of Tyr79 and surrounding residues from the internal hydrophobic cavity upon ligand binding to the solvent exposed exterior of the LTP, facilitating proteolysis. As a result, such knowledge contributes to our understanding as to how resistance to digestion can be used in allergenicity risk assessment of novel food proteins, including GMOs.

59 BASIC BIOLOGICAL SCIENCES↗

The Importance of Three-Body Interactions in Molecular Dynamics Simulations of Water with the Fragment Molecular Orbital Method

The analytic first derivative with respect to nuclear coordinates is formulated and implemented in the framework of the three-body fragment molecular orbital (FMO) method. The gradient has been derived and implemented for restricted Hartree-Fock, second-order Møller-Plesset perturbation, and density functional theories. The importance of the three-body fully analytic gradient is illustrated through the failure of the two-body FMO method during molecular dynamics simulations of a small water cluster. The parallel implementation of the fragment molecular orbital method, its parallel efficiency, and its scalability on the Blue Gene/Q architecture up to 262,144 CPU cores, are also discussed.

Pruitt, Spencer R.↗

Effective Ion Mobility Calculations for Macromolecules by Scattering off Electron Clouds

Broad commercialization and increasing resolving power of ion mobility spectrometry/mass spectrometry (IMS/MS) platforms has engendered an explosion of IMS applications to structural characterization of gas-phase biomolecules. That has renewed interest in more accurate and rapid ion mobility calculations needed to elicit ion geometries from the measurements. An approach based on scattering off electron density isosurfaces (SEDI) that mirrors the physics of molecular collisions was proven superior to the common methods involving atomic coordinates a decade ago, but has remained impractical for large ions because of extreme computational demands. Here, we accelerate SEDI by up to ~500 times using the fragment molecular orbital (FMO) approach for surface generation and the multiplexed scattering algorithm in conjunction with the new grid extrapolation procedure for cross section evaluations. Parallelization of the code on a supercomputer has produced major further speed gains, allowing SEDI calculations for proteins (defined by over a million surface points) with the precision of <0.1% in one minute. Initial tests reveal the anticipated dependences of mobility on the ion charge state and lower cross sections in view of reduced surface roughness. Present developments are expected to lead to broad application of SEDI in IMS studies of macromolecules, enabling more accurate and reliable structural assignments.

Alexeev, Yuri↗

High Accuracy Predictions of the Bulk Properties of Water (ESP Technical Report: February, 2013)

Water is the most important liquid to humanity and a deep understanding of its behavior is of critical importance to national priority scientific issues such as global warming. Despite this, accurate theory and modeling of water is lacking, and the question of the detailed mechanism behind its bulk properties remains one of the outstanding unsolved problems in science. While many simulation studies of liquid water have been undertaken, resource and capability limitations have forced the vast majority to be performed using empirical model potentials with questionable reliability. This project is the first to simulate water with high accuracy and high precision, employing recently developed advanced techniques that allow first principles ab initio quantum mechanical methods to harness hundreds of thousands of processors and scale up to thousands of water molecules.

Mayes, Maricris L.↗