-
Unsupervisedly Prompting AlphaFold2 for Few-Shot Learning of Accurate Folding Landscape and Protein Structure Prediction
Authors:
Jun Zhang,
Sirui Liu,
Mengyun Chen,
Haotian Chu,
Min Wang,
Zidong Wang,
Jialiang Yu,
Ningxi Ni,
Fan Yu,
Diqing Chen,
Yi Isaac Yang,
Boxin Xue,
Lijiang Yang,
Yuan Liu,
Yi Qin Gao
Abstract:
Data-driven predictive methods which can efficiently and accurately transform protein sequences into biologically active structures are highly valuable for scientific research and medical development. Determining accurate folding landscape using co-evolutionary information is fundamental to the success of modern protein structure prediction methods. As the state of the art, AlphaFold2 has dramatic…
▽ More
Data-driven predictive methods which can efficiently and accurately transform protein sequences into biologically active structures are highly valuable for scientific research and medical development. Determining accurate folding landscape using co-evolutionary information is fundamental to the success of modern protein structure prediction methods. As the state of the art, AlphaFold2 has dramatically raised the accuracy without performing explicit co-evolutionary analysis. Nevertheless, its performance still shows strong dependence on available sequence homologs. Based on the interrogation on the cause of such dependence, we presented EvoGen, a meta generative model, to remedy the underperformance of AlphaFold2 for poor MSA targets. By prompting the model with calibrated or virtually generated homologue sequences, EvoGen helps AlphaFold2 fold accurately in low-data regime and even achieve encouraging performance with single-sequence predictions. Being able to make accurate predictions with few-shot MSA not only generalizes AlphaFold2 better for orphan sequences, but also democratizes its use for high-throughput applications. Besides, EvoGen combined with AlphaFold2 yields a probabilistic structure generation method which could explore alternative conformations of protein sequences, and the task-aware differentiable algorithm for sequence generation will benefit other related tasks including protein design.
△ Less
Submitted 8 October, 2023; v1 submitted 20 August, 2022;
originally announced August 2022.
-
Numerical Computation of Partial Differential Equations by Hidden-Layer Concatenated Extreme Learning Machine
Authors:
Naxian Ni,
Suchuan Dong
Abstract:
The extreme learning machine (ELM) method can yield highly accurate solutions to linear/nonlinear partial differential equations (PDEs), but requires the last hidden layer of the neural network to be wide to achieve a high accuracy. If the last hidden layer is narrow, the accuracy of the existing ELM method will be poor, irrespective of the rest of the network configuration. In this paper we prese…
▽ More
The extreme learning machine (ELM) method can yield highly accurate solutions to linear/nonlinear partial differential equations (PDEs), but requires the last hidden layer of the neural network to be wide to achieve a high accuracy. If the last hidden layer is narrow, the accuracy of the existing ELM method will be poor, irrespective of the rest of the network configuration. In this paper we present a modified ELM method, termed HLConcELM (hidden-layer concatenated ELM), to overcome the above drawback of the conventional ELM method. The HLConcELM method can produce highly accurate solutions to linear/nonlinear PDEs when the last hidden layer of the network is narrow and when it is wide. The new method is based on a type of modified feedforward neural networks (FNN), termed HLConcFNN (hidden-layer concatenated FNN), which incorporates a logical concatenation of the hidden layers in the network and exposes all the hidden nodes to the output-layer nodes. HLConcFNNs have the interesting property that, given a network architecture, when additional hidden layers are appended to the network or when extra nodes are added to the existing hidden layers the representation capacity of the HLConcFNN associated with the new architecture is guaranteed to be not smaller than that of the original network architecture. Here representation capacity refers to the set of all functions that can be exactly represented by the neural network of a given architecture. We present ample benchmark tests with linear/nonlinear PDEs to demonstrate the computational accuracy and performance of the HLConcELM method and the superiority of this method to the conventional ELM from previous works.
△ Less
Submitted 15 May, 2022; v1 submitted 24 April, 2022;
originally announced April 2022.
-
Phonon modes and Raman signatures of MnBi2nTe3n+1 (n=1,2,3,4) magnetic topological heterostructures
Authors:
Yujin Cho,
Jin Ho Kang,
Liangbo Liang,
Xiangru Kong,
Subhajit Ghosh,
Fariborz Kargar,
Chaowei Hu,
Alexander A. Balandin,
Alexander A. Puretzky,
Ni Ni,
Chee Wei Wong
Abstract:
An intrinsic antiferromagnetic topological insulator $\mathrm{MnBi_2Te_4}$ can be realized by intercalating Mn-Te bilayer chain in a topological insulator, $\mathrm{Bi_2Te_3}$. $\mathrm{MnBi_2Te_4}$ provides not only a stable platform to demonstrate exotic physical phenomena, but also easy tunability of the physical properties. For example, inserting more $\mathrm{Bi_2Te_3}$ layers in between two…
▽ More
An intrinsic antiferromagnetic topological insulator $\mathrm{MnBi_2Te_4}$ can be realized by intercalating Mn-Te bilayer chain in a topological insulator, $\mathrm{Bi_2Te_3}$. $\mathrm{MnBi_2Te_4}$ provides not only a stable platform to demonstrate exotic physical phenomena, but also easy tunability of the physical properties. For example, inserting more $\mathrm{Bi_2Te_3}$ layers in between two adjacent $\mathrm{MnBi_2Te_4}$ weakens the interlayer magnetic interactions between the $\mathrm{MnBi_2Te_4}$ layers. Here we present the first observations on the inter- and intra-layer phonon modes of $\mathrm{MnBi_{2n}Te_{3n+1}}$ (n=1,2,3,4) using cryogenic low-frequency Raman spectroscopy. We experimentally and theoretically distinguish the Raman vibrational modes using various polarization configurations. The two peaks at 66 cm$^{-1}$ and 112 cm$^{-1}$ show an abnormal perturbation in the Raman linewidths below the magnetic transition temperature due to spin-phonon coupling. In $\mathrm{MnBi_4Te_7}$, the $\mathrm{Bi_2Te_3}$ layers induce Davydov splitting of the A$_{1g}$ mode around 137 cm$^{-1}$ at 5 K. Using the linear chain model, we estimate the out-of-plane interlayer force constant to be $(3.98 \pm 0.14) \times 10^{19}$ N/m$^3$ at 5 K, three times weaker than that of $\mathrm{Bi_2Te_3}$. Our work discovers the dynamics of phonon modes of the $\mathrm{MnBi_2Te_4}$ and the effect of the additional $\mathrm{Bi_2Te_3}$ layers, providing the first-principles guidance to tailor the physical properties of layered heterostructures.
△ Less
Submitted 26 July, 2021; v1 submitted 7 July, 2021;
originally announced July 2021.
-
A Method for Representing Periodic Functions and Enforcing Exactly Periodic Boundary Conditions with Deep Neural Networks
Authors:
Suchuan Dong,
Naxian Ni
Abstract:
We present a simple and effective method for representing periodic functions and enforcing exactly the periodic boundary conditions for solving differential equations with deep neural networks (DNN). The method stems from some simple properties about function compositions involving periodic functions. It essentially composes a DNN-represented arbitrary function with a set of independent periodic f…
▽ More
We present a simple and effective method for representing periodic functions and enforcing exactly the periodic boundary conditions for solving differential equations with deep neural networks (DNN). The method stems from some simple properties about function compositions involving periodic functions. It essentially composes a DNN-represented arbitrary function with a set of independent periodic functions with adjustable (training) parameters. We distinguish two types of periodic conditions: those imposing the periodicity requirement on the function and all its derivatives (to infinite order), and those imposing periodicity on the function and its derivatives up to a finite order $k$ ($k\geqslant 0$). The former will be referred to as $C^{\infty}$ periodic conditions, and the latter $C^{k}$ periodic conditions. We define operations that constitute a $C^{\infty}$ periodic layer and a $C^k$ periodic layer (for any $k\geqslant 0$). A deep neural network with a $C^{\infty}$ (or $C^k$) periodic layer incorporated as the second layer automatically and exactly satisfies the $C^{\infty}$ (or $C^k$) periodic conditions. We present extensive numerical experiments on ordinary and partial differential equations with $C^{\infty}$ and $C^k$ periodic boundary conditions to verify and demonstrate that the proposed method indeed enforces exactly, to the machine accuracy, the periodicity for the DNN solution and its derivatives.
△ Less
Submitted 14 July, 2020;
originally announced July 2020.
-
An Energy-Stable Scheme for Incompressible Navier-Stokes Equations with Periodically Updated Coefficient Matrix
Authors:
L. Lin,
N. Ni,
Z. Yang,
S. Dong
Abstract:
We present an energy-stable scheme for simulating the incompressible Navier-Stokes equations based on the generalized Positive Auxiliary Variable (gPAV) framework. In the gPAV-reformulated system the original nonlinear term is replaced by a linear term plus a correction term, where the correction term is put under control by an auxiliary variable. The proposed scheme incorporates a pressure-correc…
▽ More
We present an energy-stable scheme for simulating the incompressible Navier-Stokes equations based on the generalized Positive Auxiliary Variable (gPAV) framework. In the gPAV-reformulated system the original nonlinear term is replaced by a linear term plus a correction term, where the correction term is put under control by an auxiliary variable. The proposed scheme incorporates a pressure-correction type strategy into the gPAV procedure, and it satisfies a discrete energy stability property. The scheme entails the computation of two copies of the velocity and pressure within a time step, by solving an individual de-coupled linear equation for each of these field variables. Upon discretization the pressure linear system involves a constant coefficient matrix that can be pre-computed, while the velocity linear system involves a coefficient matrix that is updated periodically, once every $k_0$ time steps in the current work, where $k_0$ is a user-specified integer. The auxiliary variable, being a scalar-valued number, is computed by a well-defined explicit formula, which guarantees the positivity of its computed values. It is observed that the current method can produce accurate simulation results at large (or fairly large) time step sizes for the incompressible Navier-Stokes equations. The impact of the periodic coefficient-matrix update on the overall cost of the method is observed to be small in typical numerical simulations. Several flow problems have been simulated to demonstrate the accuracy and performance of the method developed herein.
△ Less
Submitted 10 September, 2019;
originally announced September 2019.
-
Optical Control of Chiral Charge Pumping in a Topological Weyl Semimetal
Authors:
M. Mehdi Jadidi,
Mehdi Kargarian,
Martin Mittendorff,
Yigit Aytac,
Bing Shen,
Jacob C. König-Otto,
Stephan Winnerl,
Ni Ni,
Alexander L. Gaeta,
Thomas E. Murphy,
H. Dennis Drew
Abstract:
Solids with topologically robust electronic states exhibit unusual electronic and optical transport properties that do not exist in other materials. A particularly interesting example is chiral charge pumping, the so-called chiral anomaly, in recently discovered topological Weyl semimetals, where simultaneous application of parallel DC electric and magnetic fields creates an imbalance in the numbe…
▽ More
Solids with topologically robust electronic states exhibit unusual electronic and optical transport properties that do not exist in other materials. A particularly interesting example is chiral charge pumping, the so-called chiral anomaly, in recently discovered topological Weyl semimetals, where simultaneous application of parallel DC electric and magnetic fields creates an imbalance in the number of carriers of opposite topological charge (chirality). Here, using time-resolved terahertz measurements on the Weyl semimetal TaAs in a magnetic field, we optically interrogate the chiral anomaly by dynamically pumping the chiral charges and monitoring their subsequent relaxation. Theory based on Boltzmann transport shows that the observed effects originate from an optical nonlinearity in the chiral charge pumping process. Our measurements reveal that the chiral population relaxation time is much greater than 1 ns. The observation of terahertz-controlled chiral carriers with long coherence times and topological protection suggests the application of Weyl semimetals for quantum optoelectronic technology.
△ Less
Submitted 6 May, 2019;
originally announced May 2019.
-
Uncovering Electron-Phonon Scattering and Phonon Dynamics in Type-I Weyl Semimetals
Authors:
Jennifer Coulter,
Gavin B. Osterhoudt,
Christina A. C. Garcia,
Yiping Wang,
Vincent Plisson,
Bing Shen,
Ni Ni,
Kenneth S. Burch,
Prineha Narang
Abstract:
Weyl semimetals are 3D phases of matter with topologically protected states that have remarkable macroscopic transport behaviors. As phonon dynamics and electron-phonon scattering play a critical role in the electrical and thermal transport, we pursue a fundamental understanding of the origin of these effects in type-I Weyl semimetals NbAs and TaAs. In the temperature-dependent Raman spectra of Nb…
▽ More
Weyl semimetals are 3D phases of matter with topologically protected states that have remarkable macroscopic transport behaviors. As phonon dynamics and electron-phonon scattering play a critical role in the electrical and thermal transport, we pursue a fundamental understanding of the origin of these effects in type-I Weyl semimetals NbAs and TaAs. In the temperature-dependent Raman spectra of NbAs, we reveal a previously unreported Fano lineshape, a signature stemming from the electron-phonon interaction. Additionally, the temperature dependence of the $A_{1}$ phonon linewidths in both NbAs and TaAs strongly deviate from the standard model of anharmonic decay. To capture the mechanisms responsible for the observed Fano asymmetry and the atypical phonon linewidth, we present first principles calculations of the phonon self-energy correction due to the electron-phonon interaction. Finally, we investigate the relationship between Fano lineshape, electron-phonon coupling, and locations of the Weyl points in these materials. Through this study of the phonon dynamics and electron-phonon interaction in these Weyl semimetals, we consider specific microscopic pathways which contribute to the nature of their macroscopic transport.
△ Less
Submitted 18 March, 2019;
originally announced March 2019.
-
Energy-Stable Boundary Conditions Based on a Quadratic Form: Applications to Outflow/Open-Boundary Problems in Incompressible Flows
Authors:
Naxian Ni,
Zhiguo Yang,
Suchuan Dong
Abstract:
We present a set of new energy-stable open boundary conditions for tackling the backflow instability in simulations of outflow/open boundary problems for incompressible flows. These boundary conditions are developed through two steps: (i) devise a general form of boundary conditions that ensure the energy stability by re-formulating the boundary contribution into a quadratic form in terms of a sym…
▽ More
We present a set of new energy-stable open boundary conditions for tackling the backflow instability in simulations of outflow/open boundary problems for incompressible flows. These boundary conditions are developed through two steps: (i) devise a general form of boundary conditions that ensure the energy stability by re-formulating the boundary contribution into a quadratic form in terms of a symmetric matrix and computing an associated eigen problem; and (ii) require that, upon imposing the boundary conditions from the previous step, the scale of boundary dissipation should match a physical scale. These open boundary conditions can be re-cast into the form of a traction-type condition, and therefore they can be implemented numerically using the splitting-type algorithm from a previous work. The current boundary conditions can effectively overcome the backflow instability typically encountered at moderate and high Reynolds numbers. These boundary conditions in general give rise to a non-zero traction on the entire open boundary, unlike previous related methods which only take effect in the backflow regions of the boundary. Extensive numerical experiments in two and three dimensions are presented to test the effectiveness and performance of the presented methods, and simulation results are compared with the available experimental data to demonstrate their accuracy.
△ Less
Submitted 18 July, 2018;
originally announced July 2018.
-
Millimeter-wave study of London penetration depth temperature dependence in Ba(Fe0.926Co0.074)2As2 single crystal
Authors:
A. A. Barannik,
N. T. Cherpak,
N. Ni,
M. A. Tanatar,
S. A. Vitusevich,
V. N. Skresanov,
P. C. Canfield,
R. Prozorov,
V. V. Glamazdin,
K. I. Torokhtii
Abstract:
In-plane surface Ka-band microwave impedance of optimally doped single crystals of the Fe-based superconductor Ba(Fe0.926Co0.074)2As2 (Tc= 22.8K) was measured. Sensitive sapphire disk quasi-optical resonator with high-Tc cuprate conducting endplates was developed specially for Fe-pnictide superconductors. It allowed finding temperature variation of London penetration depth in a form of power law,…
▽ More
In-plane surface Ka-band microwave impedance of optimally doped single crystals of the Fe-based superconductor Ba(Fe0.926Co0.074)2As2 (Tc= 22.8K) was measured. Sensitive sapphire disk quasi-optical resonator with high-Tc cuprate conducting endplates was developed specially for Fe-pnictide superconductors. It allowed finding temperature variation of London penetration depth in a form of power law, namely Δλ(T)~ Tn with n = 2.8 from low temperatures up to at least 0.6Tc consisted with radio-frequency measurements. This exponent points towards nodeless state with pairbreaking scattering, which can support one of the extended s-pairing symmetries. The dependence λ(T) at low temperatures is well described by one superconducting small-gap (Δ\cong 0.75 in kTc units, where k is Boltzman coefficient) exponential dependence.
△ Less
Submitted 30 June, 2011;
originally announced June 2011.