Catalog footprint

What is connected

95works
26topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

95 published item(s)

preprint2026arXiv

Chiral Discrimination on Gate-Based Quantum Computers

We present a novel approach to chiral discrimination using gate-based quantum processors, addressing a key challenge in adapting conventional control techniques using modern quantum computing. Schemes such as stimulated rapid adiabatic passage (STIRAP) and shortcuts to adiabaticity (STAP) have shown strong potential for enantiomer discrimination; their reliance on analog and continuous-time control makes them incompatible with digital gate-based quantum computing architectures. Here, we adapt these protocols for quantum computers by discretizing their Gaussian-shaped pulses through Trotterization. We simulate the chiral molecule 1,2-propanediol and experimentally validate this gate-based implementation on IBM quantum hardware. Our results demonstrate that this approach is a viable foundation for advancing chiral discrimination protocols, preparing the way for quantum-level manipulation of molecular chirality on accessible quantum architectures.

preprint2022arXiv

A general quantum algorithm for open quantum dynamics demonstrated with the Fenna-Matthews-Olson complex

Using quantum algorithms to simulate complex physical processes and correlations in quantum matter has been a major direction of quantum computing research, towards the promise of a quantum advantage over classical approaches. In this work we develop a generalized quantum algorithm to simulate any dynamical process represented by either the operator sum representation or the Lindblad master equation. We then demonstrate the quantum algorithm by simulating the dynamics of the Fenna-Matthews-Olson (FMO) complex on the IBM QASM quantum simulator. This work represents a first demonstration of a quantum algorithm for open quantum dynamics with a moderately sophisticated dynamical process involving a realistic biological structure. We discuss the complexity of the quantum algorithm relative to the classical method for the same purpose, presenting a decisive query complexity advantage of the quantum approach based on the unique property of quantum measurement.

preprint2022arXiv

Characterizing quantum circuits with qubit functional configurations

We propose a theory of characterizing quantum circuits with qubit functional configurations. Any quantum circuit can be decomposed into alternating sequences of 1-qubit unitary gates and CNOT gates. Each CNOT sequence prepares the current quantum state into a layer of qubit functional configuration to specify the rule for the next 1-qubit unitary sequence on how to collectively modify the state vector entries. All the functional configuration layers on a quantum circuit define its type which can include many other circuits sharing the same configuration layers. Studying the functional configuration types allows us to collectively characterize the properties and behaviors of many quantum circuits. We demonstrate the theory's application to the hardware-efficient ansatzes of variational quantum algorithms. For potential applications, the functional configuration theory may allow systematic understanding and development of quantum algorithms based on their functional configuration types.

preprint2022arXiv

Finite-Size Scaling on a Digital Quantum Simulator using Quantum Restricted Boltzmann Machine

The critical point and the critical exponents for a phase transition can be determined using the Finite-Size Scaling (FSS) analysis. This method assumes that the phase transition occurs only in the infinite size limit. However, there has been a lot of interest recently in quantum phase transitions occuring in finite size systems such as a single two-level system interacting with a single bosonic mode e.g. in the Quantum Rabi Model (QRM). Since, these phase transitions occur at a finite system size, the traditional FSS method is rendered inapplicable for these cases. For cases like this, we propose an alternative FSS method in which the truncation of the system is done in the Hilbert space instead of the physical space. This approach has previously been used to calculate the critical parameters for stability and symmetry breaking of electronic structure configurations of atomic and molecular systems. We calculate the critical point for the quantum phase transition of the QRM using this approach. We also provide a protocol to implement this method on a digital quantum simulator using the Quantum Restricted Boltzmann Machine algorithm. Our work opens up a new direction in the study of quantum phase transitions on quantum devices.

preprint2022arXiv

Geometrical picture of the electron-electron correlation at the large-D limit

In electronic structure calculations, the correlation energy is defined as the difference between the mean field and the exact solution of the non relativistic Schrödinger equation. Such an error in the different calculations is not directly observable as there is no simple quantum mechanical operator, apart from correlation functions, that correspond to such quantity. Here, we use the dimensional scaling approach, in which the electrons are localized at the large-dimensional scaled space, to describe a geometric picture of the electronic correlation. Both, the mean field, and the exact solutions at the large-D limit have distinct geometries. Thus, the difference might be used to describe the correlation effect. Moreover, correlations can be also described and quantified by the entanglement between the electrons, which is a strong correlation without a classical analog. Entanglement is directly observable and it is one of the most striking properties of quantum mechanics and bounded by the area law for local gapped Hamiltonians of interacting many-body systems. This study opens up the possibility of presenting a geometrical picture of the electron-electron correlations and might give a bound on the correlation energy. The results at the large-D limit and at D=3 indicate the feasibility of using the geometrical picture to get a bound on the electron-electron correlations.

preprint2022arXiv

Quantum Machine Learning for Chemistry and Physics

Machine learning (ML) has emerged into formidable force for identifying hidden but pertinent patterns within a given data set with the objective of subsequent generation of automated predictive behavior. In the recent years, it is safe to conclude that ML and its close cousin deep learning (DL) have ushered unprecedented developments in all areas of physical sciences especially chemistry. Not only the classical variants of ML , even those trainable on near-term quantum hardwares have been developed with promising outcomes. Such algorithms have revolutionzed material design and performance of photo-voltaics, electronic structure calculations of ground and excited states of correlated matter, computation of force-fields and potential energy surfaces informing chemical reaction dynamics, reactivity inspired rational strategies of drug designing and even classification of phases of matter with accurate identification of emergent criticality. In this review we shall explicate a subset of such topics and delineate the contributions made by both classical and quantum computing enhanced machine learning algorithms over the past few years. We shall not only present a brief overview of the well-known techniques but also highlight their learning strategies using statistical physical insight. The objective of the review is to not only to foster exposition to the aforesaid techniques but also to empower and promote cross-pollination among future-research in all areas of chemistry which can benefit from ML and in turn can potentially accelerate the growth of such algorithms.

preprint2022arXiv

The unitary dependence theory for understanding quantum circuits and states

We develop a unitary dependence theory to characterize the behaviors of quantum circuits and states in terms of how quantum gates manipulate qubits and determine their measurement probabilities. A qubit has dependence on a 1-qubit unitary gate if its measurement probabilities depend on the parameters of the gate. A 1-qubit unitary creates such a dependence onto the target qubit, and a CNOT gate copies all the dependences from the control qubit to the target qubit. The complete dependence picture of the output state details the connections qubits may have when being manipulated or measured. Compared to the conventional entanglement description of quantum circuits and states, the dependence picture offers more practical information, easier generalization to many-qubit systems, and better robustness upon partitioning of the system. The unitary dependence theory is a useful tool for understanding quantum circuits and states that is based on the practical ideas of how manipulations and measurements are affected by elementary gates.

preprint2021arXiv

A quantum encryption design featuring confusion, diffusion, and mode of operation

Quantum cryptography -- the application of quantum computing techniques to cryptography has been extensively investigated. Two major directions of quantum cryptography are quantum key distribution (QKD) and quantum encryption, with the former focusing on secure key distribution and the latter focusing on encryption using quantum algorithms. In contrast to the success of the QKD, the development of quantum encryption algorithms is limited to designs of mostly one-time pads (OTP) that are unsuitable for most communication needs. In this work we propose a non-OTP quantum encryption scheme utilizing a quantum state creation process to encrypt messages. As essentially a non-OTP quantum block cipher the method stands out against existing methods with the following features: 1. complex key-ciphertext relation (i.e. confusion) and complex plaintext-ciphertext relation (i.e. diffusion); 2. mode of operation design for practical encryption on multiple blocks. These features provide key reusability and protection against eavesdropping and standard cryptanalytic attacks.

preprint2021arXiv

Non-Randomness of Google's Quantum Supremacy Benchmark

The first achievement of quantum supremacy has been claimed recently by Google for the random quantum circuit benchmark with 53 superconducting qubits. Here, we analyze the randomness of Google's quantum random-bit sampling. The heat maps of Google's random bit-strings show stripe patterns at specific qubits in contrast to the Haar-measure or classical random-bit strings. Google's data contains more bit 0 than bit 1, i.e., about 2.8\% difference, and fail to pass the NIST random number tests, while the Haar-measure or classical random-bit samples pass. Their difference is also illustrated by the Marchenko-Pastur distribution and the Girko circular law of random matrices of random bit-strings. The calculation of the Wasserstein distances shows that Google's random bit-strings are farther away from the Haar-measure random bit-strings than the classical random bit-strings. Our results imply that random matrices and the Wasserstein distance could be new tools for analyzing the performance of quantum computers.

preprint2021arXiv

Relaxation of a Stationary State on a Quantum Computer Yields Unique Spectroscopic Fingerprint of the Computer's Noise

Quantum computing has the potential to revolutionize computing for certain classes of problems with exponential scaling, and yet this potential is accompanied by significant sensitivity to noise, requiring sophisticated error correction and mitigation strategies. Here we simulate the relaxations of stationary states at different frequencies on several quantum computers to obtain unique spectroscopic fingerprints of their noise. Response functions generated from the data reveal a clear signature of non-Markovian dynamics, demonstrating that each of the quantum computers acts as a non-Markovian bath with a unique colored noise profile. The study suggest that noisy intermediate-scale quantum computers (NISQ) provide a built-in noisy bath that can be analyzed from their simulation of closed quantum systems with the results potentially being harnessed for error mitigation or open-system simulation.

preprint2021arXiv

Statistical Approach to Quantum Phase Estimation

We introduce a new statistical and variational approach to the phase estimation algorithm (PEA). Unlike the traditional and iterative PEAs which return only an eigenphase estimate, the proposed method can determine any unknown eigenstate-eigenphase pair from a given unitary matrix utilizing a simplified version of the hardware intended for the Iterative PEA (IPEA). This is achieved by treating the probabilistic output of an IPEA-like circuit as an eigenstate-eigenphase proximity metric, using this metric to estimate the proximity of the input state and input phase to the nearest eigenstate-eigenphase pair and approaching this pair via a variational process on the input state and phase. This method may search over the entire computational space, or can efficiently search for eigenphases (eigenstates) within some specified range (directions), allowing those with some prior knowledge of their system to search for particular solutions. We show the simulation results of the method with the Qiskit package on the IBM Q platform and on a local computer.

preprint2020arXiv

Characterization of quantum states based on creation complexity

The creation complexity of a quantum state is the minimum number of elementary gates required to create it from a basic initial state. The creation complexity of quantum states is closely related to the complexity of quantum circuits, which is crucial in developing efficient quantum algorithms that can outperform classical algorithms. A major question unanswered so far is what quantum states can be created with a number of elementary gates that scales polynomially with the number of qubits. In this work we first show for an entirely general quantum state it is exponentially hard (requires a number of steps that scales exponentially with the number of qubits) to determine if the creation complexity is polynomial. We then show it is possible for a large class of quantum states with polynomial creation complexity to have common coefficient features such that given any candidate quantum state we can design an efficient coefficient sampling procedure to determine if it belongs to the class or not with arbitrarily high success probability. Consequently partial knowledge of a quantum state's creation complexity is obtained, which can be useful for designing quantum circuits and algorithms involving such a state.

preprint2020arXiv

Gap states and valley-spin filtering in transition metal dichalcogenide monolayers

The magnetically-induced valley-spin filtering in transition metal dichalcogenide monolayers ($MX_{2}$, where $M$=Mo, W and $X$=S, Se, Te) promises new paradigm in information processing. However, the detailed understanding of this effect is still limited, regarding its underlying transport processes. Herein, it is suggested that the filtering mechanism can be greately elucidated by the concept of metal-induced gap states (MIGS), appearing in the electrode-terminated $MX_{2}$ materials {\it i.e.} the referential filter setup. In particular, the gap states are predicted here to mediate valley- and spin-resolved charge transport near the ideal electrode/$MX_{2}$ interface, and therefore to initiate filtering. It is also argued that the role of MIGS increases when the channel length is diminished, as they begin to govern the overall valley-spin transport in the tunneling regime. In what follows, the presented study yields fundamental scaling trends for the valley-spin selectivity with respect to the intrinsic physics of the filter materials. As a result, it facilitates insight into the analyzed effects and provide design guidelines toward efficient valley-spin filter devices, that base on the discussed materials or other hexagonal monolayers with a broken inversion symmetry.

preprint2020arXiv

Hybrid quantum-classical neural network for calculating ground state energies of molecules

We present a hybrid quantum classical neural network that can be trained to perform electronic structure calculation and generate potential energy curves of simple molecules. The method is based on the combination of parameterized quantum circuits and measurements. With unsupervised training, the neural network can generate electronic potential energy curves based on training at certain bond lengths. To demonstrate the power of the proposed new method, we present results of using the quantum-classical hybrid neural network to calculate ground state potential energy curves of simple molecules such as H$_2$, LiH and BeH$_2$. The results are very accurate and the approach could potentially be used to generate complex molecular potential energy surfaces.

preprint2020arXiv

Introducing Open boundary conditions in modeling nonperiodic materials and interfaces: the impact of the periodic assumption

Simulations are essential to accelerate the discovery of new materials and to gain full understanding of known ones. Although hard to realize experimentally, periodic boundary conditions are omnipresent in material simulations. In this work, we intro-duce ROBIN (recursive open boundary and interfaces), the first method allowing open boundary conditions in material and interface modeling. The computational costs are limited to solving quantum properties in a focus area which allows explicitly discretizing millions of atoms in real space and to consider virtually any type of environment (be it periodic, regular, or ran-dom). The impact of the periodicity assumption is assessed in detail with silicon dopants in graphene. Graphene was con-firmed to produce a band gap with periodic substitution of 3% carbon with silicon in agreement with published periodic boundary condition calculations. Instead, 3% randomly distributed silicon in graphene only shifts the energy spectrum. The predicted shift agrees quantitatively with published experimental data. Key insight of this assessment is, assuming periodici-ty elevates a small perturbation of a periodic cell into a strong impact on the material property prediction. Periodic boundary conditions can be applied on truly periodic systems only. More general systems should apply an open boundary method for reliable predictions.

preprint2020arXiv

Magnetic flux noise in superconducting qubits and the gap states continuum

In the present study we investigate the selected local aspects of the metal-induced gap states (MIGSs) at the disordered metal-insulator interface, that were previously proposed to produce magnetic moments responsible for the magnetic flux noise in some of the superconducting qubit modalities. Our analysis attempts to supplement the available studies and provide new theoretical contribution toward their validation. In particular, we explicitly discuss the behavior of the MIGSs in the momentum space as a function of the local onsite energy deviation, that mimics random potential disorder at the interface. It is found, that when the difference between the characteristic electronic potentials in the insulator increases, the corresponding MIGSs become more localized. This effect is associated with the increasing degree of the potential disorder that was earlier observed to produce highly localized MIGSs in the superconducting qubits. At the same time, the presented findings show that the disorder-induced localization of the MIGSs can be related directly to the decay characteristics of these states as well as to the bulk electronic properties of the insulator. As a result, our study reinforces plausibility of the previous corresponding investigations on the origin of the flux noise, but also allows to draw future directions toward their better verification.

preprint2020arXiv

Maximal entropy approach for quantum state tomography

Quantum computation has been growing rapidly in both theory and experiments. In particular, quantum computing devices with a large number of qubits have been developed by IBM, Google, IonQ, and others. The current quantum computing devices are noisy intermediate-scale quantum $($NISQ$)$ devices, and so approaches to validate quantum processing on these quantum devices are needed. One of the most common ways of validation for an n-qubit quantum system is quantum tomography, which tries to reconstruct a quantum system's density matrix by a complete set of observables. However, the inherent noise in the quantum systems and the intrinsic limitations poses a critical challenge to precisely know the actual measurement operators which make quantum tomography impractical in experiments. Here, we propose an alternative approach to quantum tomography, based on the maximal information entropy, that can predict the values of unknown observables based on the available mean measurement data. This can then be used to reconstruct the density matrix with high fidelity even though the results for some observables are missing. Of additional contexts, a practical approach to the inference of the quantum mechanical state using only partial information is also needed.

preprint2020arXiv

Unorthodox dimensional interpolations for He, Li, Be atoms and hydrogen molecule

We present a simple interpolation formula using dimensional limits $D=1$ and $D=\infty$ to obtain the $D=3$ ground-state energies of atoms and molecules. For atoms, these limits are linked by first-order perturbation terms of electron-electron interactions. This unorthodox approach is illustrated by ground-states for two, three, and four electron atoms, with modest effort to obtain fairly accurate results. Also, we treat the ground-state of H$_2$ over a wide range of the internuclear distance R, and compares well with the standard exact results from the Full Configuration Interaction method. Similar dimensional interpolations may be useful for complex many-body systems.

preprint2019arXiv

A quantum algorithm for evolving open quantum dynamics on quantum computing devices

Designing quantum algorithms for simulating quantum systems has seen enormous progress, yet few studies have been done to develop quantum algorithms for open quantum dynamics despite its importance in modeling the system-environment interaction found in most realistic physical models. In this work we propose and demonstrate a general quantum algorithm to evolve open quantum dynamics on quantum computing devices. The Kraus operators governing the time evolution can be converted into unitary matrices with minimal dilation guaranteed by the Sz.-Nagy theorem. This allows the evolution of the initial state through unitary quantum gates, while using significantly less resource than required by the conventional Stinespring dilation. We demonstrate the algorithm on an amplitude damping channel using the IBM Qiskit quantum simulator and the IBM Q 5 Tenerife quantum device. The proposed algorithm does not require particular models of dynamics or decomposition of the quantum channel, and thus can be easily generalized to other open quantum dynamical models.

preprint2019arXiv

Enhancement of Photovoltaic Current Generation through Dark States in Donor-Acceptor Pairs of Tungsten-based Transition Metal Di-Chalcogenides (TMDCs)

As several photovoltaic materials experimentally approach the Shockley-Queisser limit, there has been a growing interest in unconventional materials and approaches with the potential to cross this efficiency barrier. One such candidate is dark state protection induced by the dipole-dipole interaction between molecular excited states. This phenomenon has been shown to significantly reduce carrier recombination rate and enhance photon-to-current conversion, in elementary models consisting of few interacting chromophore centers. Atomically thin 2D transition metal di-chalcogenides (TMDCs) have shown great potential for use as ultra-thin photovoltaic materials in solar cells due to their favorable photon absorption and electronic transport properties. TMDC alloys exhibit tunable direct bandgaps and significant dipole moments. In this work, we introduce the dark state protection mechanism to a TMDC based photovoltaic system with pure tungsten diselenide (WSe2) as the acceptor material and the TMDC alloy tungsten sulfo-selenide (WSeS) as the donor material. Our numerical model demonstrates the first application of the dark state protection mechanism to a photovoltaic material with a photon current enhancement of up to 35% and an ideal photon-to-current efficiency exceeding the Shockley-Queisser limit.

preprint2018arXiv

Connecting bright and dark states through accidental degeneracy caused by lack of symmetry

Coupled excitonic structures are found in natural and artificial light harvesting systems where optical transitions link different excitation manifolds. In systems with symmetry, some optical transitions are allowed, while others are forbidden. Here we examine an excitonic ring structure and identify an accidental degeneracy between two categories of double-excitation eigenstates with distinct symmetries and optical transition properties. To understand the accidental degeneracy, a complete selection rule between two arbitrary excitation manifolds is derived with a physically motivated proof. Remarkably, symmetry analysis shows the lack of certain symmetry elements in the Hamiltonian is responsible for this degeneracy, which is unique to rings with size N=4l+2 (l being an integer).

preprint2018arXiv

Double-excitation manifold's effect on exciton transfer dynamics and the efficiency of coherent light harvesting

The efficiency of natural light harvesting systems is largely determined by their ability to transfer excitations from the antenna to the energy trapping center before recombination. Exciton diffusion length similarly limits organic photovoltaics and demands bulk heterojunction architectures. Dark state protection, achieved by coherent coupling between subunits within the antenna, can significantly reduce radiative recombination and enhance the efficiency of energy trapping. In this work we extend the dark state concept to the double-excitation manifold by studying the dynamical flow of excitations. We show the lowest double-excitation state carries minimal oscillator strength but relaxation to this state from higher lying double excitations can be relatively rapid such that the lowest double excitation state can act as a dynamical dark state protecting excitation from radiative recombination. This mechanism is sensitive to topology and operates differently for chain and ring structures, while becoming more pronounced in both geometries when the size of the antenna increases. When the exciton-exciton annihilation (EEA) mechanism is considered, the double-excitation population is quickly depleted and the dynamics changes dramatically. However the efficiency and output power are still significantly different from those calculated using the single-excitation manifold alone, justifying the necessity of considering the double-excitation manifold. Remarkably, in certain scenarios, the EEA can even increase the overall light harvesting efficiency by bringing population down from the double-excitation dark states to the single-excitation manifold.

preprint2017arXiv

Dark states and delocalization: competing effects of quantum coherence on the efficiency of light harvesting systems

Natural light harvesting systems exploit electronic coupling of identical chromophores to generate efficient and robust excitation transfer and conversion. Dark states created by strong coupling between chromophores in the antenna structure can significantly reduce radiative recombination and enhance energy conversion efficiency. Increasing the number of the chromophores increases the number of dark states and the associated enhanced energy conversion efficiency, yet also delocalizes excitations away from the trapping center and reduces the energy conversion rate. Therefore, a competition between dark state protection and delocalization must be considered when designing the optimal size of a light harvesting system. In this study, we explore the two competing mechanisms in a chain-structured antenna and show that dark state protection is the dominant mechanism, with an intriguing dependence on the parity of the number of chromophores. This dependence is linked to the exciton distribution among eigenstates which is strongly affected by the coupling strength between chromophores and the temperature. Combining these findings, we propose that increasing the coupling strength between the chromophores can significantly increase the power output of the light harvesting system.

preprint2016arXiv

An Ancilla Based Quantum Simulation Framework for Non-Unitary Matrices

The success probability in an ancilla based circuit generally decreases exponentially in the number of qubits consisted in the ancilla. Although the probability can be amplified through the amplitude amplification process, the input dependence of the amplitude amplification makes difficult to sequentially combine two or more ancilla based circuits. A new version of the amplitude amplification known as the oblivious amplitude amplification runs independently of the input to the system register. This allow us to to sequentially combine two or more ancilla based circuits. However, this type of the amplification only works when the considered system is unitary or non-unitary but somehow close to a unitary. In this paper, we present a general framework to simulate non-unitary matrices on ancilla based quantum circuits in which the success probability is maximized by using the oblivious amplitude amplification. In particular, we show how to extend a non-unitary matrix to an almost unitary matrix. We then simulate the extended matrix by using an ancilla based circuit design along with the oblivious amplitude amplification. Measuring the distance of the produced matrix to the closest unitary unitary matrix, a lower bound for the fidelity of the final state obtained from the oblivious amplitude amplification process is presented. Numerical simulations for random matrices of different sizes show that independent of the system size, the final amplified probabilities are generally around 0.75 and the fidelity of the final state is mostly high and around 0.95. Furthermore, we discuss the complexity analysis and show that combining two such ancilla based circuits, a matrix product can be implemented. This may lead us to efficiently implement matrix functions represented as infinite matrix products on quantum computers.

preprint2016arXiv

Confirmation of Kramers-Henneberger Atoms

In a remarkable experiment, Eichmann et al. attained unprecedented acceleration of neutral atoms by strong short-pulse IR laser fields. The driving mechanism was identified as the ponderomotive force on excited electrons bound in Rydberg orbits that survive long enough to enable the atoms to reach the detector. However, the observed velocities lie somewhat above the theoretical prediction. The systematic discrepancy was attributed to "absolute laser intensity uncertainties or a slightly non-Gaussian intensity distribution". Here, we examine the process by transforming to the Kramers-Henneberger (KH) reference frame. We find that in addition to the ponderomotive potential there exists a smaller but significant term that comes from the binding energy of the KH atom. Including this KH term brings the calculated maximum velocities to a close match with experimental results over the full range of laser pulse durations.

preprint2016arXiv

Efficient estimation of perturbative error with cellular automata

From celestial mechanics to quantum theory of atoms and molecules, perturbation theory has played a central role in natural sciences. Particularly in quantum mechanics, the amount of information needed for specifying the state of a many-body system commonly scales exponentially as the system size. This poses a fundamental difficulty in using perturbation theory at arbitrary order. As one computes the terms in the perturbation series at increasingly higher orders, it is often important to determine whether the series converges and if so, what is an accurate estimation of the total error that comes from the next order of perturbation up to infinity. Here we present a set of efficient algorithms that compute tight upper bounds to perturbation terms at arbitrary order. We argue that these tight bounds often take the form of symmetric polynomials on the parameter of the quantum system. We then use cellular automata as our basic model of computation to compute the symmetric polynomials that account for all of the virtual transitions at any given order. At any fixed order, the computational cost of our algorithm scales polynomially as a function of the system size. We present a non-trivial example which shows that our error estimation is nearly tight with respect to exact calculation.

preprint2016arXiv

Electronic Structure Calculatins and the Ising Machine

The exact solution of Schrodinger equation for atoms, molecules and extended systems continues to be a "Holy Grail" problem that the entire field has been striving to solve since its inception. Recently, breakthroughs have been made in the development of quantum annealing and coherent Ising machines capable of simulating hundreds of connected spins interacting with an Ising type Hamiltonian. One of the most vital questions pertaining to these new devices is: 'can these machine be used to perform electronic structure calculations?' Here we discuss the general procedure used by these devices and show that there is an exact mapping between the electronic structure Hamiltonian of the Hydrogen molecule and the Ising Hamiltonian

preprint2016arXiv

Enhancing Intrinsic Stability of Hybrid Perovskite Solar Cell by Strong, yet Balanced, Electronic Coupling

In the past few years, the meteoric development of hybrid organic--inorganic perovskite solar cells (PSC) astonished the community. The efficiency has already reached the level needed for commercialization; however, the instability hinders its deployment on the market. Here, we report a mechanism to chemically stabilize PSC absorbers. We propose to replace the widely used methylammonium cation (\ce{CH3NH3+}) by alternative molecular cations allowing an enhanced electronic coupling between the cation and the \ce{PbI6} octahedra while maintaining the band gap energy within the suitable range for solar cells. The mechanism exploits establishing a balance between the electronegativity of the materials' constituents and the resulting ionic electrostatic interactions. The calculations demonstrate the concept of enhancing the electronic coupling, and hence the stability, by exploring the stabilizing features of \ce{CH3PH3+}, \ce{CH3SH2+}, and \ce{SH3+} cations, among several other possible candidates. Chemical stability enhancement hence results from a strong, yet balanced, electronic coupling between the cation and the halides in the octahedron. This shall unlock the hindering instability problem for PSCs and allow them to hit the market as a serious low-cost competitor to silicon based solar cell technologies.

preprint2016arXiv

Hydrogen Bonding: A Mechanism for Tuning Electronic and Optical Properties of Hybrid Organic-Inorganic Frameworks

Here we demonstrate that significant progress in this area may be achieved by introducing structural elements that form hydrogen bonds with environment. Considering several examples of hybrid framework materials with different structural ordering containing protonated sulfonium cation H3S+ that forms strong hydrogen bonds with electronegative halogen anions (Cl-, F-), we found that hydrogen bonding increases the structural stability of the material and may be used for tuning electronic states near the bandgap. We suggest that such a behavior has a universal character and should be observed in hybrid inorganic-organic framework materials containing protonated cations. This effect may serve as a viable route for optoelectronic and photovoltaic applications.

preprint2016arXiv

Prospects for Quantum Computing with an Array of Ultracold Polar Paramagnetic Molecules

Arrays of trapped ultracold molecules represent a promising platform for implementing a universal quantum computer. DeMille has detailed a prototype design based on Stark states of polar $^1Σ$ molecules as qubits. Herein, we consider an array of polar $^2Σ$ molecules which are, in addition, inherently paramagnetic and whose Hund's case (b) free-rotor states are Bell states. We show that by subjecting the array to combinations of concurrent homogeneous and inhomogeneous electric and magnetic fields, the entanglement of the array's Stark and Zeeman states can be tuned and the qubit sites addressed. Two schemes for implementing an optically controlled CNOT gate are proposed and their feasibility discussed in the face of the broadening of spectral lines due to dipole-dipole coupling and the inhomogeneity of the electric and magnetic fields.

preprint2016arXiv

Quantum computation using arrays of N polar molecules in pendular states

We investigate several aspects of realizing quantum computation using entangled polar molecules in pendular states. Quantum algorithms typically start from a product state |00...0> and we show that up to a negligible error, the ground states of polar molecule arrays can be considered as the unentangled qubit basis state |00...0>. This state can be prepared by simply allowing the system to reach thermal equilibrium at low temperature (<1 mK). We also evaluate entanglement, characterized by the concurrence of pendular state qubits in dipole arrays as governed by the external electric field, dipole-dipole coupling and number N of molecules in the array. In the parameter regime that we consider for quantum computing, we find that qubit entanglement is modest, typically no greater than 0.0001, confirming the negligible entanglement in the ground state. We discuss methods for realizing quantum computation in the gate model, measurement based model, instantaneous quantum polynomial time circuits and the adiabatic model using polar molecules in pendular states.

preprint2016arXiv

Solving Set Cover with Pairs Problem Using Quantum Annealing

Here we consider using quantum annealing to solve Set Cover with Pairs (SCP), an NP-hard combinatorial optimization problem that play an important role in networking, computational biology, and biochemistry. We show an explicit construction of Ising Hamiltonians whose ground states encode the solution of SCP instances. We numerically simulate the time-dependent Schrödinger equation in order to test the performance of quantum annealing for random instances and compare with that of simulated annealing. We also discuss explicit embedding strategies for realizing our Hamiltonian construction on the D-wave type restricted Ising Hamiltonian based on Chimera graphs. Our embedding on the Chimera graph preserves the structure of the original SCP instance and in particular, the embedding for general complete bipartite graphs and logical disjunctions may be of broader use than that the specific problem we deal with.

preprint2015arXiv

Correction to kinetic energy density using exactly solvable model

An accurate non-gradient-expansion based correction to Thomas--Fermi is developed using solvable model. The used model is a system of $N$ non-interacting electrons moving independently in the Coulomb field of the nuclear charge. The presented correction is applicable for atoms and should be extendable beyond that. The method exploits the fact that the difference between the Thomas--Fermi approximation and the non-interacting kinetic energy is comparable to the difference between the same values inside the proposed solvable model. The numerical experiments show that by adding this correction factor, the precision of Thomas--Fermi approximation is enhanced by an order of magnitude.

preprint2015arXiv

Domain Walls Conductivity in Hybrid Organometallic Perovskites: The Key of CH3NH3PbI3 Solar Cell High Performance

The past several years has witnessed a surge of interest in organometallic trihalide perovskites, which are at the heart of the new generation of solid-state solar cells. Here, we calculated the static conductivity of charged domain walls in n- and p- doped organometallic uniaxial ferroelectric semiconductor perovskite CH3NH3PbI3 using the Landau-Ginzburg-Devonshire (LGD) theory. We find that due to the charge carrier accumulation, the static conductivity may drastically increase at the domain wall by 3-4 orders of magnitude in comparison with conductivity through the bulk of the material. Also, a two-dimensional degenerated gas of highly mobile charge carriers could be formed at the wall. The high values of conductivity at domain walls and interfaces explain high efficiency in organometallic solution-processed perovskite films which contains lots of different point and extended defects. These results could suggest new routes to enhance the performance of this promising class of novel photovoltaic materials.

preprint2015arXiv

Electron Localization and Energy Levels' Oscillations Induced by Controlled Deformation

Manipulating energy levels while controlling the electron localization is an essential step for many applications of confined systems. In this paper we demonstrate how to achieve electron localization and induce energy level oscillation in one-dimensional quantum systems by externally controlling the deformation of the system. From a practical point of view, the one-dimensional potentials can be realized using layered structures. In the analysis, we considered three different examples. The first one is a graded quantum well between confining infinite walls where the deformation is modeled by varying slightly the graded well. The second systems is a symmetric multiple quantum well between infinite walls under the effect of biasing voltage. The third system is a layered 2D hybrid perovskites where pressure is used to induce deformation. The calculations are conducted both numerically and analytically using the perturbation theory. It is shown that the obtained oscillations are associated with level avoided crossings and that the deformation results in changing the spatial localization of the electrons.

preprint2015arXiv

Generalized Remote Preparation of Arbitrary $m$-qubit Entangled States via Genuine Entanglements

Herein, we present a feasible, general protocol for quantum communication within a network via generalized remote preparation of an arbitrary $m$-qubit entangled state designed with genuine tripartite Greenberger--Horne--Zeilinger-type entangled resources. During the implementations, we construct novel collective unitary operations; these operations are tasked with performing the necessary phase transfers during remote state preparations. We have distilled our implementation methods into a five-step procedure, which can be used to faithfully recover the desired state during transfer. Compared to previous existing schemes, our methodology features a greatly increased success probability. After the consumption of auxiliary qubits and the performance of collective unitary operations, the probability of successful state transfer is increased four-fold and eight-fold for arbitrary two- and three-qubit entanglements when compared to other methods within the literature, respectively. We conclude this paper with a discussion of the presented scheme for state preparation, including: success probabilities, reducibility and generalizability.

preprint2015arXiv

Hamiltonian gadgets with reduced resource requirements

Application of the adiabatic model of quantum computation requires efficient encoding of the solution to computational problems into the lowest eigenstate of a Hamiltonian that supports universal adiabatic quantum computation. Experimental systems are typically limited to restricted forms of 2-body interactions. Therefore, universal adiabatic quantum computation requires a method for approximating quantum many-body Hamiltonians up to arbitrary spectral error using at most 2-body interactions. Hamiltonian gadgets, introduced around a decade ago, offer the only current means to address this requirement. Although the applications of Hamiltonian gadgets have steadily grown since their introduction, little progress has been made in overcoming the limitations of the gadgets themselves. In this experimentally motivated theoretical study, we introduce several gadgets which require significantly more realistic control parameters than similar gadgets in the literature. We employ analytical techniques which result in a reduction of the resource scaling as a function of spectral error for the commonly used subdivision, 3- to 2-body and $k$-body gadgets. Accordingly, our improvements reduce the resource requirements of all proofs and experimental proposals making use of these common gadgets. Next, we numerically optimize these new gadgets to illustrate the tightness of our analytical bounds. Finally, we introduce a new gadget that simulates a $YY$ interaction term using Hamiltonians containing only $\{X,Z,XX,ZZ\}$ terms. Apart from possible implications in a theoretical context, this work could also be useful for a first experimental implementation of these key building blocks by requiring less control precision without introducing extra ancillary qubits.

preprint2015arXiv

Neuroreceptor Activation by Vibration-Assisted Tunneling

G protein-coupled receptors (GPCRs) constitute a large family of receptor proteins that sense molecular signals on the exterior of a cell and activate signal transduction pathways within the cell. Modeling how an agonist activates such a receptor is fundamental for an understanding of a wide variety of physiological processes and it is of tremendous value for pharmacology and drug design. Inelastic electron tunneling spectroscopy (IETS) has been proposed as a model for the mechanism by which olfactory GPCRs are activated by a bound agonist. We apply this hypothesis to GPCRs within the mammalian nervous system using quantum chemical modeling. We found that non-endogenous agonists of the serotonin receptor share a particular IET spectral aspect both amongst each other and with the serotonin molecule: a peak whose intensity scales with the known agonist potencies. We propose an experiential validation of this model by utilizing lysergic acid dimethylamide (DAM-57), an ergot derivative, and its deuterated isotopologues; we also provide theoretical predictions for comparison to experiment. If validated our theory may provide new avenues for guided drug design and elevate methods of in silico potency/activity prediction.

preprint2015arXiv

Revealing the role of organic cations in hybrid halide perovskites CH3NH3PbI3

The hybrid halide perovskite CH$_{3}$NH$_{3}$PbI$_{3}$ has enabled solar cells to reach an efficiency of about 18\%, demonstrating a pace for improvements with no precedents in the solar energy arena. Despite such explosive progress, the microscopic origin behind the success of such material is still debated, with the role played by the organic cations in the light-harvesting process remaining unclear. Here van-der-Waals-corrected density functional theory calculations reveal that the orientation of the organic molecules plays a fundamental role in determining the material electronic properties. For instance, if CH$_{3}$NH$_{3}$ orients along a (011)-like direction, the PbI$_{6}$ octahedral cage will distort and the band gap will become indirect. Our results suggest that molecular rotations, with the consequent dynamical change of the band structure, might be at the origin of the slow carrier recombination and the superior conversion efficiency of CH$_{3}$NH$_{3}$PbI$_{3}$.

preprint2015arXiv

Singularity of the time-energy uncertainty in adiabatic perturbation and cycloids on a Bloch sphere

The adiabatic perturbation is shown to be singular from the exact solution of a spin-1/2 particle in a uniformly rotating magnetic field. Due to a non-adiabatic effect, its quantum trajectory on a Bloch sphere is a cycloid traced by a circle rolling along an adiabatic path. As the magnetic field rotates more and more slowly, the time-energy uncertainty, proportional to the distance of the quantum trajectory, calculated by the exact solution is entirely different from the one obtained by the adiabatic path traced by the instantaneous state. However, the non-adiabatic Aharonov-Anandan geometric phase, measured by the area enclosed by the exact path, approaches smoothly the adiabatic Berry phase, proportional to the area enclosed by the adiabatic path. The singular limit of the time-energy uncertainty and the regular limit of the geometric phase are associated with the arc length and arc area of the cycloid on a Bloch sphere, respectively. Prolate and curtate cycloids are also traced by different initial states outside and inside of the rolling circle, respectively. The axis trajectory of the rolling circle, parallel to the adiabatic path, is shown to be an example of transitionless driving. The non-adiabatic resonance is visualized by the number of complete cycloid arcs.

preprint2015arXiv

The Radical Pair Mechanism and the Avian Chemical Compass: Quantum Coherence and Entanglement

We review the spin radical pair mechanism which is a promising explanation of avian navigation. This mechanism is based on the dependence of product yields on (1) the hyperfine interaction involving electron spins and neighboring nuclear spins and (2) the intensity and orientation of the geomagnetic field. One surprising result is that even at ambient conditions quantum entanglement of electron spins can play an important role in avian magnetoreception. This review describes the general scheme of chemical reactions involving radical pairs generated from singlet and triplet precursors; the spin dynamics of the radical pairs; and the magnetic field dependence of product yields caused by the radical pair mechanism. The main part of the review includes a description of the chemical compass in birds. We review: the general properties of the avian compass; the basic scheme of the radical pair mechanism; the reaction kinetics in cryptochrome; quantum coherence and entanglement in the avian compass; and the effects of noise. We believe that the "quantum avian compass" can play an important role in avian navigation and can also provide the foundation for a new generation of sensitive and selective magnetic-sensing nano-devices.

preprint2014arXiv

A density-matrix renormalization group Study of one-dimensional incommensurate quantum Frenkel-Kontorova model

In this paper, the one-dimensional incommensurate quantum Frenkel-Kontorova model is investigate by a density-matrix renormalization group algorithm. Special attention is given to the entanglement and the ground state energy. The energy gap between the ground state and the first excited is also calculated. From all the numerical results, we have observed an obvious property changes from the pinned state to the sliding one as the quantum fluctuation is increased. But no expected quantum critical point can be justified by the present data.

preprint2014arXiv

Cuckoo Search Inspired Hybridization of the Nelder-Mead Simplex Algorithm Applied to Optimization of Photovoltaic Cells

A new hybridization of the Cuckoo Search (CS) is developed and applied to optimize multi-cell solar systems; namely multi-junction and split spectrum cells. The new approach consists of combining the CS with the Nelder-Mead method. More precisely, instead of using single solutions as nests for the CS, we use the concept of a simplex which is used in the Nelder-Mead algorithm. This makes it possible to use the flip operation introduces in the Nelder-Mead algorithm instead of the Levy flight which is a standard part of the CS. In this way, the hybridized algorithm becomes more robust and less sensitive to parameter tuning which exists in CS. The goal of our work was to optimize the performance of multi-cell solar systems. Although the underlying problem consists of the minimization of a function of a relatively small number of parameters, the difficulty comes from the fact that the evaluation of the function is complex and only a small number of evaluations is possible. In our test, we show that the new method has a better performance when compared to similar but more compex hybridizations of Nelder-Mead algorithm using genetic algorithms or particle swarm optimization on standard benchmark functions. Finally, we show that the new method outperforms some standard meta-heuristics for the problem of interest.

preprint2014arXiv

Delocalized Quantum States Enhance Photocell Efficiency

The high quantum efficiency of photosynthetic complexes has inspired researchers to explore new routes to utilize this process for photovoltaic devices. Quantum coherence has been demonstrated to play a crucial role within this process. Herein, we propose a three-dipole system as a model of a new photocell type which exploits the coherence among its three dipoles. We have proved that the efficiency of such a photocell is greatly enhanced by quantum coherence. We have also predicted that the photocurrents can be enhanced by about 49.5 % in such a coherent coupled dipole system compared with the uncoupled dipoles. These results suggest a promising novel design aspect of photosynthesis-mimicking photovoltaic devices.

preprint2014arXiv

Efficient method for localized functions using domain transformation and Fourier sine series

An efficient approach to handle localized states by using spectral methods (SM) in one and three dimensions is presented. The method consists of transformation of the infinite domain to the bounded domain in $(0, π)$ and using the Fourier sine series as a set of basis functions for the SM. It is shown that with an appropriate choice of transformation functions, this method manages to preserve the good properties of original SMs; more precisely, superb computational efficiency when high level of accuracy is necessary. This is made possible by analytically exploiting the properties of the transformation function and the Fourier sine series. An especially important property of this approach is the possibility of calculating the Hartree energy very efficiently. This is done by exploiting the positive properties of the sine series as a basis set and conducting an extinctive part of the calculations analytically. We illustrate the efficiency of this method and implement it to solve the Poisson's and Helmholtz equations in both one and three dimensions. The efficiency of the method is verified through a comparison to recently published results for both one and three dimensional problems.

preprint2014arXiv

Error analysis of the density-matrix renormalization group algorithm for a chain of harmonic oscillators

We investigate the application of the density-matrix renormalization group (DMRG) algorithm to a one-dimensional harmonic oscillator chain and compare the results with exact solutions, aiming to improve the algorithm efficiency. It has been demonstrated that the algorithm can give quite accurate results if the procedure is proper organized, for example, by using the optimized bases. The errors of calculated ground state energy and the energy gap between the ground state and the first excited state are analyzed, which are found to be critically dependent upon the size of the system or the energy level structure of the studied system and the number of states targeted during the DMRG procedure.

preprint2014arXiv

Infrared-dressed entanglement of cold open-shell polar molecules for universal matchgate quantum computing

Implementing a scalable quantum information processor using polar molecules in optical lattices requires precise control over the long-range dipole-dipole interaction between molecules in selected lattice sites. We present here a scheme using trapped open-shell $^2Σ$ polar molecules that allows dipolar exchange processes between nearest and next-nearest neighbors to be controlled to construct a generalized transverse Ising spin Hamiltonian with tunable $XX$, $YY$ and $XY$ couplings in the rotating frame of the driving lasers. The scheme requires a moderately strong bias magnetic field with near-infrared light to provide local tuning of the qubit energy gap, and mid-infrared pulses to perform rotational state transfer via stimulated Raman adiabatic passage. No interaction between qubits is present in the absence of the infrared driving. We analyze the fidelity of the resulting two-qubit matchgate, and demonstrate its robustness as a function of the driving parameters. We discuss a realistic application of the system for universal matchgate quantum computing in optical lattices.

preprint2014arXiv

Multiparty-controlled remote preparation of four-qubit cluster-type entangled states

We present a strategy for implementing multiparty-controlled remote state preparation (MCRSP) for a family of four-qubit cluster-type states with genuine entanglements while employing, Greenberg-Horne-Zeilinger-class states as quantum channels. In this scenario, the encoded information is transmitted from the sender to a spatially separated receiver via the control of multi-party. Predicated on the collaboration of all participants, the desired state can be entirely restored within the receiver's place with high success probability by application of appropriate local operations and necessary classical communication . Moreover, this proposal for MCRSP can be faithfully achieved with unit total success probability when the quantum channels are distilled to maximally entangled ones.

preprint2014arXiv

Multiple Network Alignment on Quantum Computers

Comparative analyses of graph structured datasets underly diverse problems. Examples of these problems include identification of conserved functional components (biochemical interactions) across species, structural similarity of large biomolecules, and recurring patterns of interactions in social networks. A large class of such analyses methods quantify the topological similarity of nodes across networks. The resulting correspondence of nodes across networks, also called node alignment, can be used to identify invariant subgraphs across the input graphs. Given $k$ graphs as input, alignment algorithms use topological information to assign a similarity score to each $k$-tuple of nodes, with elements (nodes) drawn from each of the input graphs. Nodes are considered similar if their neighbors are also similar. An alternate, equivalent view of these network alignment algorithms is to consider the Kronecker product of the input graphs, and to identify high-ranked nodes in the Kronecker product graph. Conventional methods such as PageRank and HITS (Hypertext Induced Topic Selection) can be used for this purpose. These methods typically require computation of the principal eigenvector of a suitably modified Kronecker product matrix of the input graphs. We adopt this alternate view of the problem to address the problem of multiple network alignment. Using the phase estimation algorithm, we show that the multiple network alignment problem can be efficiently solved on quantum computers. We characterize the accuracy and performance of our method, and show that it can deliver exponential speedups over conventional (non-quantum) methods.

preprint2014arXiv

Non-Adiabatic Quantum Dynamics of Grover's Adiabatic Search Algorithm

We study quantum dynamics of Grover's adiabatic search algorithm with the equivalent two-level system. Its adiabatic and non-adiabatic evolutions are visualized as trajectories of Bloch vectors on a Bloch sphere. We find the change in the non-adiabatic transition probability from exponential decay for short running time to inverse-square decay for long running time. The size dependence of the critical running time is expressed in terms of Lambert $W$ function. The transitionless driving Hamiltonian is obtained to make a quantum state follow the adiabatic path. We demonstrate that a constant Hamiltonian, approximate to the exact time-dependent driving Hamiltonian, can alter the non-adiabatic transition probability from the inverse square decay to the inverse fourth power decay with running time. This may open up a new way of reducing errors in adiabatic quantum computation.

preprint2014arXiv

Probing the Role of the Eighth Bacteriochlorophyll in holo-FMO Complex by Simulated Two-Dimensional Electronic Spectroscopy

The Fenna-Matthews-Olson (FMO) protein-pigment complex acts as a molecular wire between the outer antenna system and the reaction center (RC); it is an important model system to study the excitonic energy transfer. Recent crystallographic studies report the existence of an additional (eighth) bacteriochlorophyll a (BChl a). To understand the functionality of this eighth BChl, we simulated the two-dimensional electronic spectra of both the 7-site (apo form) and the 8-site (holo form) variant of the FMO complex from green sulfur bacteria, Prosthecochloris aestuarii. By comparing the difference between the spectrum, it was found that the eighth BChl can affect two different excitonic energy transfer pathways, these being: (1) directly involve in the first pathway 6 $\rightarrow$ 3 $\rightarrow$ 1 of the apo form model by passing the excitonic energy to exciton 6; and (2) increase the excitonic wave function overlap between excitons 4 and 5 in the second pathway (7 $\rightarrow$ 4,5 $\rightarrow$ 2 $\rightarrow$ 1) and thus increase the possible downward sampling routes across the BChls.

preprint2014arXiv

Quadratic constrained mixed discrete optimization with an adiabatic quantum optimizer

We extend the family of problems that may be implemented on an adiabatic quantum optimizer (AQO). When a quadratic optimization problem has at least one set of discrete controls and the constraints are linear, we call this a quadratic constrained mixed discrete optimization (QCMDO) problem. QCMDO problems are NP-hard, and no efficient classical algorithm for their solution is known. Included in the class of QCMDO problems are combinatorial optimization problems constrained by a linear partial differential equation (PDE) or system of linear PDEs. An essential complication commonly encountered in solving this type of problem is that the linear constraint may introduce many intermediate continuous variables into the optimization while the computational cost grows exponentially with problem size. We resolve this difficulty by developing a constructive mapping from QCMDO to quadratic unconstrained binary optimization (QUBO) such that the size of the QUBO problem depends only on the number of discrete control variables. With a suitable embedding, taking into account the physical constraints of the realizable coupling graph, the resulting QUBO problem can be implemented on an existing AQO. The mapping itself is efficient, scaling cubically with the number of continuous variables in the general case and linearly in the PDE case if an efficient preconditioner is available.

preprint2014arXiv

Quantum phase transition in one-dimensional commensurate Frenkel-Kontorova model

In this paper, we have studied the one-dimensional commensurate quantum Frenkel-Kontorova model by a density-matrix renormalization group (DMRG) algorithm. The focus has been on its properties of the entanglement, the coordinate correlation, the ground state energy and the energy gap between the ground state and the first excited one. It is demonstrated that a quantum phase transition (QPT) between the pinned and the sliding phases takes place as the quantum fluctuation, measured by an effective Planck constant $\tilde\hbar$, increases to a threshold value $\tilde \hbar_c$.

preprint2014arXiv

Quantum Random State Generation with Predefined Entanglement Constraint

Entanglement plays an important role in quantum communication, algorithms, and error correction. Schmidt coefficients are correlated to the eigenvalues of the reduced density matrix. These eigenvalues are used in Von Neumann entropy to quantify the amount of the bipartite entanglement. In this paper, we map the Schmidt basis and the associated coefficients to quantum circuits to generate random quantum states. We also show that it is possible to adjust the entanglement between subsystems by changing the quantum gates corresponding to the Schmidt coefficients. In this manner, random quantum states with predefined bipartite entanglement amounts can be generated using random Schmidt basis. This provides a technique for generating equivalent quantum states for given weighted graph states, which are very useful in the study of entanglement, quantum computing, and quantum error correction.

preprint2014arXiv

Sensitivity and Entanglement in the Avian Chemical Compass

The Radical Pair Mechanism can help to explain avian orientation and navigation. Some evidence indicates that the intensity of external magnetic fields plays an important role in avian navigation. In this paper, based on a two-stage strategy, we demonstrate that birds could reasonably detect the directions of geomagnetic fields and gradients of these fields using a yield-based chemical compass that is sensitive enough for navigation. Also, we find that the lifetime of entanglement in this proposed compass is angle-dependent and long enough to allow adequate electron transfer between molecules.

preprint2014arXiv

Theoretical Limits of Photovoltaics Efficiency and Possible Improvements by Intuitive Approaches Learned from Photosynthesis and Quantum Coherence

In this review, we present and discussed the main trends in photovoltaics with emphasize on the conversion efficiency limits. The theoretical limits of various photovoltaics device concepts are presented and analyzed using a flexible detailed balance model where more discussion emphasize is toward the losses. Also, few lessons from nature and other fields to improve the conversion efficiency in photovoltaics are presented and discussed as well. From photosynthesis, the perfect exciton transport in photosynthetic complexes can be utilized for PVs. Also, we present some lessons learned from other fields like recombination suppression by quantum coherence. For example, the coupling in photosynthetic reaction centers is used to suppress recombination in photocells.

preprint2013arXiv

Avalanches in the Raise and Peel model in the presence of a wall

We investigate a non-equilibrium one-dimensional model known as the raise and peel model describing a growing surface which grows locally and has non-local desorption. For specific values of adsorption ($u_a$) and desorption($u_d$) rates the model shows interesting features. At $u_a = u_d$, the model is described by a conformal field theory (with conformal charge $c=0$) and its stationary probability canbe mapped to the ground state of the XXZ quantum chain. Moreover, for $u_a \geq u_d$, the model shows a phase in which the the avalanche distribution is scale invariant. In this work we study the surface dynamics by looking at avalanche distributions using Finite-size Scaling formalism and explore the effect of adding a wall to the model. The model shows the same universality for the cases with and without a wall for an odd number of tiles removed, but we find a new exponent in the presence of a wall for an even number of avalanches released. We provide new conjecture for the probability distribution of avalanches with a wall obtained by using exact diagonalization of small lattices and Monte-Carlo simulations.

preprint2013arXiv

Decomposition of Unitary Matrices for Finding Quantum Circuits: Application to Molecular Hamiltonians

Constructing appropriate unitary matrix operators for new quantum algorithms and finding the minimum cost gate sequences for the implementation of these unitary operators is of fundamental importance in the field of quantum information and quantum computation. Evolution of quantum circuits faces two major challenges: complex and huge search space and the high costs of simulating quantum circuits on classical computers. Here, we use the group leaders optimization algorithm to decompose a given unitary matrix into a proper-minimum cost quantum gate sequence. We test the method on the known decompositions of Toffoli gate, the amplification step of the Grover search algorithm, the quantum Fourier transform, and the sender part of the quantum teleportation. Using this procedure, we present the circuit designs for the simulation of the unitary propagators of the Hamiltonians for the hydrogen and the water molecules. The approach is general and can be applied to generate the sequence of quantum gates for larger molecular systems.

preprint2013arXiv

Degree Distribution in Quantum Walks on Complex Networks

In this theoretical study, we analyze quantum walks on complex networks, which model network-based processes ranging from quantum computing to biology and even sociology. Specifically, we analytically relate the average long time probability distribution for the location of a unitary quantum walker to that of a corresponding classical walker. The distribution of the classical walker is proportional to the distribution of degrees, which measures the connectivity of the network nodes and underlies many methods for analyzing classical networks including website ranking. The quantum distribution becomes exactly equal to the classical distribution when the walk has zero energy and at higher energies the difference, the so-called quantumness, is bounded by the energy of the initial state. We give an example for which the quantumness equals a Renyi entropy of the normalized weighted degrees, guiding us to regimes for which the classical degree-dependent result is recovered and others for which quantum effects dominate.

preprint2013arXiv

Dynamics of Entanglement in one and two-dimensional spin systems

We present a review of dynamics of entanglement in one and two dimensional systems under the effect of external magnetic field and different degrees of anisotropy at zero and finite temperatures. Different techniques for treating the spin systems with large Hilbert space dimensions are discussed such as trace minimization algorithm, time-evolution matrix transformation and step by step projection technique.

preprint2013arXiv

Entanglement creation in cold molecular gases using strong laser pulses

While many-particle entanglement can be found in natural solids and strongly interacting atomic and molecular gases, generating highly entangled states between weakly interacting particles in a controlled and scalable way presents a significant challenge. We describe here a one-step method to generate entanglement in a dilute gas of cold polar molecules. For molecules in optical traps separated by a few micrometers, we show that maximally entangled states can be created using the strong off-resonant pulses that are routinely used in molecular alignment experiments. We show that the resulting alignment-mediated entanglement can be detected by measuring laser-induced fluorescence with single-site resolution and that signatures of this molecular entanglement also appear in the microwave absorption spectra of the molecular ensemble. We analyze the robustness of these entangled molecular states with respect to intensity fluctuations of the trapping laser and discuss possible applications of the system for quantum information processing.

preprint2013arXiv

Experimental realization of quantum algorithm for solving linear systems of equations

Quantum computers have the potential of solving certain problems exponentially faster than classical computers. Recently, Harrow, Hassidim and Lloyd proposed a quantum algorithm for solving linear systems of equations: given an $N\times{N}$ matrix $A$ and a vector $\vec b$, find the vector $\vec x$ that satisfies $A\vec x = \vec b$. It has been shown that using the algorithm one could obtain the solution encoded in a quantum state $|x$ using $O(\log{N})$ quantum operations, while classical algorithms require at least O(N) steps. If one is not interested in the solution $\vec{x}$ itself but certain statistical feature of the solution ${x}|M|x$ ($M$ is some quantum mechanical operator), the quantum algorithm will be able to achieve exponential speedup over the best classical algorithm as $N$ grows. Here we report a proof-of-concept experimental demonstration of the quantum algorithm using a 4-qubit nuclear magnetic resonance (NMR) quantum information processor. For all the three sets of experiments with different choices of $\vec b$, we obtain the solutions with over 96% fidelity. This experiment is a first implementation of the algorithm. Because solving linear systems is a common problem in nearly all fields of science and engineering, we will also discuss the implication of our results on the potential of using quantum computers for solving practical linear systems.

preprint2013arXiv

Influence of the intensity gradient upon HHG from free electrons scattered by an intense laser beam

When an electron is scattered by a tightly-focused laser beam in vacuum, the intensity gradient is a critical factor to influence the electron dynamics, for example, the electron energy exchange with the laser fields as have been explored before [P.X.Wang et al.,J. Appl. Phys. 91, 856 (2002]. In this paper, we have further investigated its influence upon the electron high-harmonic generation (HHG) by treating the spacial gradient of the laser intensity as a ponderomotive potential. Based upon perturbative QED calculations, it has been found that the main effect of the intensity gradient is the broadening of the originally line HHG spectra. A one-to-one relationship can be built between the beam width and the corresponding line width. Hence this finding may provides us a promising way to measure the beam width of intense lasers in experiments. In addition, for a laser pulse, we have also studied the different influences from transverse and longitudinal intensity gradients upon HHG.

preprint2013arXiv

Manipulation of Molecules with Electromagnetic Fields

The goal of the present article is to review the major developments that have led to the current understanding of molecule-field interactions and experimental methods for manipulating molecules with electromagnetic fields. Molecule-field interactions are at the core of several, seemingly distinct, areas of molecular physics. This is reflected in the organization of this article, which includes sections on Field control of molecular beams, External field traps for cold molecules, Control of molecular orientation and molecular alignment, Manipulation of molecules by non-conservative forces, Ultracold molecules and ultracold chemistry, Controlled many-body phenomena, Entanglement of molecules and dipole arrays, and Stability of molecular systems in high-frequency super-intense laser fields. The article contains 853 references.

preprint2013arXiv

Persistence of entanglement in thermal states of spin systems

We study and compare the persistence of bipartite entanglement (BE) and multipartite entanglement (ME) in one-dimensional and two-dimensional spin XY models in an external transverse magnetic field under the effect of thermal excitations. We compare the threshold temperature at which the entanglement vanishes in both types of entanglement. We use the entanglement of formation as a measure of the BE and the geometric measure to evaluate the ME of the system. We have found that in both dimensions in the anisotropic and partially anisotropic spin systems at zero temperatures, all types of entanglement decay as the magnetic field increases but are sustained with very small magnitudes at high field values. Also we found that for the same systems, the threshold temperatures of the nearest neighbour (nn) BEs are higher than both of the next-to-nearest neighbour BEs and MEs and the three of them increase monotonically with the magnetic field strength. Thus, as the temperature increases, the ME and the far parts BE of the system become more fragile to thermal excitations compared to the nn BE. For the isotropic system, all types of entanglement and threshold temperatures vanish at the same exact small value of the magnetic field. We emphasise the major role played by both the properties of the ground state of the system and the energy gap in controlling the characteristics of the entanglement and threshold temperatures. In addition, we have shown how an inserted magnetic impurity can be used to preserve all types of entanglement and enhance their threshold temperatures. Furthermore, we found that the quantum effects in the spin systems can be maintained at high temperatures, as the different types of entanglements in the spin lattices are sustained at high temperatures by applying sufficiently high magnetic fields.

preprint2013arXiv

Quantum Circuit Design for Solving Linear Systems of Equations

Recently, it is shown that quantum computers can be used for obtaining certain information about the solution of a linear system Ax=b exponentially faster than what is possible with classical computation. Here we first review some key aspects of the algorithm from the standpoint of finding its efficient quantum circuit implementation using only elementary quantum operations, which is important for determining the potential usefulness of the algorithm in practical settings. Then we present a small-scale quantum circuit that solves a 2x2 linear system. The quantum circuit uses only 4 qubits, implying a tempting possibility for experimental realization. Furthermore, the circuit is numerically simulated and its performance under different circuit parameter settings is demonstrated.

preprint2013arXiv

Quantum Confinement and Negative Heat Capacity

Thermodynamics dictates that the specific heat of a system is strictly non-negative. However, in finite classical systems there are well known theoretical and experimental cases where this rule is violated, in particular finite atomic clusters. Here, we show for the first time that negative heat capacity can also occur in finite quantum systems. The physical scenario on which this effect might be experimentally observed is discussed. Observing such an effect might lead to the design of new light harvesting nano devices, in particular a solar nano refrigerator.

preprint2013arXiv

Quantum criticality analysis by finite size scaling and exponential basis sets

We combine the finite size scaling method with the meshfree spectral method to calculate quantum critical parameters for a given Hamiltonian. The basic idea is to expand the exact wave function in a finite exponential basis set and extrapolate the information about system criticality from a finite basis to the infinite basis set limit. The used exponential basis set -though chosen intuitively- allows handling a very wide range of exponential decay rates and calculating multiple eigenvalues simultaneously. As a benchmark system to illustrate the combined approach, we choose the Hulthen potential. The results show that the method is very accurate and converges faster when compared with other basis functions. The approach is general and can be extended to examine near threshold phenomena for atomic and molecular systems based on even-tempered exponential and Gaussian basis functions.

preprint2013arXiv

Reducing the number of ancilla qubits and the gate count required for creating large controlled operations

In this paper we show that it is possible to adapt a qudit scheme for creating a controlled-Toffoli created by Ralph et al. [Phys. Rev. A 75 011213] to be applicable to qubits. While this scheme requires more gates than standard schemes for creating large controlled gates, we show that with simple adaptations it is directly equivalent to the standard scheme in the literature. This scheme is the most gate-efficient way of creating large controlled unitaries currently known, however it is expensive in terms of the number of ancilla qubits used. We go on to show that using a combination of these standard techniques presented by Barenco et al. [Phys. Rev. A 52 3457 (1995)] we can create an n-qubit version of the Toffoli using less gates and the same number of ancilla qubits as recent work using computer optimization. This would be useful in any architecture of quantum computing where gates are cheap but qubit initialization is expensive.

preprint2012arXiv

Designing Efficient Programmable Quantum Circuits

Constructing general programmable circuits to be able to run any given unitary operator efficiently on a quantum processor is of fundamental importance. We present a new quantum circuit design technique resulting two general programmable circuit schemes. The circuit schemes can be used to simulate any given operator by setting the angle values in the circuit. This provides a fixed circuit design whose angles are determined from the elements of the given matrix, which can be non-unitary, in an efficient way. We also give both classical and quantum complexity analysis for these circuits and show that the circuits require a few classical computations, and the quantum complexities of them are almost the same as non-general circuits.

preprint2012arXiv

Finite size scaling for quantum criticality using the finite-element method

Finite size scaling for the Schrödinger equation is a systematic approach to calculate the quantum critical parameters for a given Hamiltonian. This approach has been shown to give very accurate results for critical parameters by using a systematic expansion with global basis-type functions. Recently, the finite element method was shown to be a powerful numerical method for ab initio electronic structure calculations with a variable real-space resolution. In this work, we demonstrate how to obtain quantum critical parameters by combining the finite element method (FEM) with finite size scaling (FSS) using different ab initio approximations and exact formulations. The critical parameters could be atomic nuclear charges, internuclear distances, electron density, disorder, lattice structure, and external fields for stability of atomic, molecular systems and quantum phase transitions of extended systems. To illustrate the effectiveness of this approach we provide detailed calculations of applying FEM to approximate solutions for the two-electron atom with varying nuclear charge; these include Hartree-Fock, density functional theory under the local density approximation, and an "exact"' formulation using FEM. We then use the FSS approach to determine its critical nuclear charge for stability; here, the size of the system is related to the number of elements used in the calculations. Results prove to be in good agreement with previous Slater-basis set calculations and demonstrate that it is possible to combine finite size scaling with the finite-element method by using ab initio calculations to obtain quantum critical parameters. The combined approach provides a promising first-principles approach to describe quantum phase transitions for materials and extended systems.

preprint2012arXiv

Ground-state stability and criticality of two-electron atoms with screened Coulomb potentials using the B-splines basis set

We applied the finite-size scaling method using the B-splines basis set to construct the stability diagram for two-electron atoms with a screened Coulomb potential. The results of this method for two electron atoms are very accurate in comparison with previous calculations based on Gaussian, Hylleraas, and finite-element basis sets. The stability diagram for the screened two-electron atoms shows three distinct regions: a two-electron region, a one-electron region, and a zero-electron region, which correspond to stable, ionized and double ionized atoms. In previous studies, it was difficult to extend the finite size scaling calculations to large molecules and extended systems because of the computational cost and the lack of a simple way to increase the number of Gaussian basis elements in a systematic way. Motivated by recent studies showing how one can use B-splines to solve Hartree-Fock and Kohn-Sham equations, this combined finite size scaling using the B-splines basis set, might provide an effective systematic way to treat criticality of large molecules and extended systems. As benchmark calculations, the two-electron systems show the feasibility of this combined approach and provide an accurate reference for comparison.

preprint2012arXiv

Implementation of Quantum Logic Gates Using Polar Molecules in Pendular States

We present a systematic approach to implementation of basic quantum logic gates operating on polar molecules in pendular states as qubits for a quantum computer. A static electric field prevents quenching of the dipole moments by rotation, thereby creating the pendular states; also, the field gradient enables distinguishing among qubit sites. Multi-Target Optimal Control Theory (MTOCT) is used as a means of optimizing the initial-to-target transition probability via a laser field. We give detailed calculations for the SrO molecule, a favorite candidate for proposed quantum computers. Our simulation results indicate that NOT, Hadamard and CNOT gates can be realized with high fidelity for such pendular qubit states.

preprint2012arXiv

Multipartite Quantum Entanglement Evolution in Photosynthetic Complexes

We investigate the evolution of entanglement in the Fenna-Matthew-Olson (FMO) complex based on simulations using the scaled hierarchy equation of motion (HEOM) approach. We examine the role of multipartite entanglement in the FMO complex by direct computation of the convex roof optimization for a number of measures, including some that have not been previously evaluated. We also consider the role of monogamy of entanglement in these simulations. We utilize the fact that the monogamy bounds are saturated in the single exciton subspace. This enables us to compute more measures of entanglement exactly and also to validate the evaluation of the convex roof. We then use direct computation of the convex roof to evaluate measures that are not determined by monogamy. This approach provides a more complete account of the entanglement in these systems than has been available to date. Our results support the hypothesis that multipartite entanglement is maximum primary along the two distinct electronic energy transfer pathways.

preprint2012arXiv

Population and Coherence Dynamics in Light Harvesting Complex II (LH2)

The electronic excitation population and coherence dynamics in the chromophores of the photosynthetic light harvesting complex 2 (LH2) B850 ring from purple bacteria (Rhodopseudomonas acidophila) have been studied theoretically at both physiological and cryogenic temperatures. Similar to the well-studied Fenna-Matthews-Olson (FMO) protein, oscillations of the excitation population and coherence in the site basis are observed in LH2 by using a scaled hierarchical equation of motion (HEOM) approach. However, this oscillation time (300 fs) is much shorter compared to the FMO protein (650 fs) at cryogenic temperature. Both environment and high temperature are found to enhance the propagation speed of the exciton wave packet yet they shorten the coherence time and suppress the oscillation amplitude of coherence and the population. Our calculations show that a long-lived coherence between chromophore electronic excited states can exist in such a noisy biological environment.

preprint2012arXiv

Quantum algorithm and circuit design solving the Poisson equation

The Poisson equation occurs in many areas of science and engineering. Here we focus on its numerical solution for an equation in d dimensions. In particular we present a quantum algorithm and a scalable quantum circuit design which approximates the solution of the Poisson equation on a grid with error \varepsilon. We assume we are given a supersposition of function evaluations of the right hand side of the Poisson equation. The algorithm produces a quantum state encoding the solution. The number of quantum operations and the number of qubits used by the circuit is almost linear in d and polylog in \varepsilon^{-1}. We present quantum circuit modules together with performance guarantees which can be also used for other problems.

preprint2012arXiv

Quantum Coherence and Entanglement in the Avian Compass

The radical pair mechanism is one of two distinct mechanisms used to explain the navigation of birds in geomagnetic fields. However, little research has been done to explore the role of quantum entanglement in this mechanism. In this paper, we study the lifetime of radical pair entanglement corresponding to the magnitude and direction of magnetic fields to show that the entanglement lasts long enough in birds to be used for navigation. We also demonstrate that, due to a lack of orientational sensitivity of the entanglement in the geomagnetic field, the birds are not able to orient themselves by the mechanism based directly on radical-pair entanglement. To explore the entanglement mechanism further, we propose a model in which the hyperfine interactions are replaced by local magnetic fields of similar strength. The entanglement of the radical pair in this model lasts longer and displays an angular sensitivity in weak magnetic fields, both of these factors are not present in the previous models.

preprint2012arXiv

Tuning entanglement and ergodicity in two-dimensional spin systems using impurities and anisotropy

We consider the entanglement in a two dimensional XY model in an external magnetic field h. The model consists of a set of 7 localized spin-1/2 particles in a two dimensional triangular lattice coupled through nearest neighbor exchange interaction J. We examine the effect of single and double impurities in the system as well as the degree of anisotropy on the nearest neighbor entanglement and ergodicity of the system. We have found that the entanglement of the system at the different degrees of anisotropy mimics that of the one dimensional spin systems at the extremely small and large values of the parameter h/J. The entanglement of the Ising and partially anisotropic system show phase transition in the vicinity of h/J = 2 while the entanglement of the isotropic system suddenly vanishes there. Also we investigate the dynamic response of the system containing single and double impurities to an external exponential magnetic field at different degrees of anisotropy. We have demonstrated that the ergodicity of the system can be controlled by varying the strength and location of the impurities as well as the degree of anisotropy of the coupling.

preprint2012arXiv

Universal Programmable Quantum Circuit Schemes to Emulate an Operator

Unlike fixed designs, programmable circuit designs support an infinite number of operators. The functionality of a programmable circuit can be altered by simply changing the angle values of the rotation gates in the circuit. Here, we present a new quantum circuit design technique resulting in two general programmable circuit schemes. The circuit schemes can be used to simulate any given operator by setting the angle values in the circuit. This provides a fixed circuit design whose angles are determined from the elements of the given matrix-which can be non-unitary-in an efficient way. We also give both the classical and quantum complexity analysis for these circuits and show that the circuits require a few classical computations. They have almost the same quantum complexities as non-general circuits. Since the presented circuit designs are independent from the matrix decomposition techniques and the global optimization processes used to find quantum circuits for a given operator, high accuracy simulations can be done for the unitary propagators of molecular Hamiltonians on quantum computers. As an example, we show how to build the circuit design for the hydrogen molecule.

preprint2011arXiv

Dimensional Scaling Treatment with Relativistic Corrections for Stable Multiply Charged Atomic Ions in High-Frequency Super-Intense Laser Fields

We present a theoretical framework which describes multiply-charged atomic ions, their stability within super-intense laser fields,also lay corrections to the systems due to relativistic effects. Dimensional scaling calculations with relativistic corrections for systems: H, H$^{-}$, H$^{2-}$, He, He$^{-}$, He$^{2-}$, He$^{3-}$ within super-intense laser fields were completed. Also completed were three-dimensional self consistent field calculations to verify the dimensionally scaled quantities. With the aforementioned methods the system's ability to stably bind 'additional' electrons through the development of multiple isolated regions of high potential energy leading to nodes of high electron density is shown. These nodes are spaced far enough from each other to minimized the electronic repulsion of the electrons, while still providing adequate enough attraction so as to bind the excess elections into orbitals. We have found that even with relativistic considerations these species are stably bound within the field. It was also found that performing the dimensional scaling calculations for systems within the confines of laser fields to be a much simpler and more cost-effective method than the supporting D=3 SCF method. The dimensional scaling method is general and can be extended to include relativistic corrections to describe the stability of simple molecular systems in super-intense laser fields.

preprint2011arXiv

Dynamics of entanglement in a two-dimensional spin system

We consider the time evolution of entanglement in a finite two dimensional transverse Ising model. The model consists of a set of 7 localized spin-1/2 particles in a two dimensional triangular lattice coupled through nearest neighbor exchange interaction in presence of an external time dependent magnetic field. The magnetic field is applied in different function forms: step, exponential, hyperbolic and periodic. We found that the time evolution of the entanglement shows an ergodic behavior under the effect of the time dependent magnetic fields. Also we found that while the step magnetic field causes great disturbance to the system creating rabid oscillations, the system shows great controllability under the effect of the other magnetic fields where the entanglement profile follows closely the shape of the applied field even with the same frequency for periodic fields. This follow up trend breaks down as the strength of the field, the transition constant for exponential and hyperbolic, or frequency for periodic field increase leading to rapid oscillations. We observed that the entanglement is very sensitive to the initial value of the applied periodic field, the smaller the initial value the less distorted is the entanglement profile. Furthermore, the effect of thermal fluctuations is very devastating to the entanglement which decays very rapidly as the temperature increases. Interestingly, although large value of the magnetic field strength may yield small entanglement, it was found to be more persistent against thermal fluctuations than the small field strengths.

preprint2011arXiv

Entanglement dynamics of one-dimensional driven spin systems in time-varying magnetic fields

We study the dynamics of entanglement for a one-dimensional spin chain with a nearest neighbor time dependent Heisenberg coupling J(t) between the spins in presence of a time dependent external magnetic field h(t) at zero and finite temperatures. We consider different forms of time dependence for the coupling and magnetic field; exponential, hyperbolic and periodic. We examined the system size effect on the entanglement asymptotic value. It was found that for a small system size the entanglement starts to fluctuate within a short period of time after applying the time dependent coupling. The period of time increases as the system size increases and disappears completely as the size goes to infinity. We also found that when J(t) is periodic the entanglement shows a periodic behavior with the same period, which disappears upon applying periodic magnetic field with the same frequency. Solving the particular case where J(t) and h(t) are proportional exactly, we showed that the asymptotic value of entanglement depends only on the initial conditions regardless of the form of J(t) and h(t) applied at t > 0.

preprint2011arXiv

Entanglement of polar molecules in pendular states

In proposals for quantum computers using arrays of trapped ultracold polar molecules as qubits, a strong external field with appreciable gradient is imposed in order to prevent quenching of the dipole moments by rotation and to distinguish among the qubit sites. That field induces the molecular dipoles to undergo pendular oscillations, which markedly affect the qubit states and the dipole-dipole interaction. We evaluate entanglement of the pendular qubit states for two linear dipoles, characterized by pairwise concurrence, as a function of the molecular dipole moment and rotational constant, strengths of the external field and the dipole-dipole coupling, and ambient temperature. We also evaluate a key frequency shift, $\triangleω$, produced by the dipole-dipole interaction. Under conditions envisioned for the proposed quantum computers, both the concurrence and $\triangleω$ become very small for the ground eigenstate. In principle, such weak entanglement can be sufficient for operation of logic gates, provided the resolution is high enough to detect the $\triangleω$ shift unambiguously. In practice, however, for many candidate polar molecules it appears a challenging task to attain adequate resolution. Simple approximate formulas fitted to our numerical results are provided from which the concurrence and $\triangleω$ shift can be obtained in terms of unitless reduced variables.

preprint2011arXiv

Entanglement of polar symmetric top molecules as candidate qubits

Proposals for quantum computing using rotational states of polar molecules as qubits have previously considered only diatomic molecules. For these the Stark effect is second-order, so a sizable external electric field is required to produce the requisite dipole moments in the laboratory frame. Here we consider use of polar symmetric top molecules. These offer advantages resulting from a first-order Stark effect, which renders the effective dipole moments nearly independent of the field strength. That permits use of much lower external field strengths for addressing sites. Moreover, for a particular choice of qubits, the electric dipole interactions become isomorphous with NMR systems for which many techniques enhancing logic gate operations have been developed. Also inviting is the wider chemical scope, since many symmetric top organic molecules provide options for auxiliary storage qubits in spin and hyperfine structure or in internal rotation states.

preprint2011arXiv

Group Leaders Optimization Algorithm

We present a new global optimization algorithm in which the influence of the leaders in social groups is used as an inspiration for the evolutionary technique which is designed into a group architecture. To demonstrate the efficiency of the method, a standard suite of single and multidimensional optimization functions along with the energies and the geometric structures of Lennard-Jones clusters are given as well as the application of the algorithm on quantum circuit design problems. We show that as an improvement over previous methods, the algorithm scales as N^2.5 for the Lennard-Jones clusters of N-particles. In addition, an efficient circuit design is shown for two qubit Grover search algorithm which is a quantum algorithm providing quadratic speed-up over the classical counterpart.

preprint2011arXiv

Modified-scaled hierarchical equation of motion approach for the study of quantum coherence in photosynthetic complexes

We present a detailed theoretical study of the transfer of electronic excitation energy through the Fenna-Matthews-Olson (FMO) pigment-protein complex, using the new developed modified scaled hierarchical approach [Shi Q. et al, J Chem Phys 2009, 130, 084105]. We show that this approach is computationally more efficient than the original hierarchical approach. The modified approach reduces the truncation levels of the auxiliary density operators and the correlation function. We provide a systematic study of how the number of auxiliary density operators and the higher-order correlation functions affect the exciton dynamics. The time scales of the coherent beating are consistent with experimental observations. Furthermore, our theoretical results exhibit population beating at physiological temperature. Additionally, the method does not require a low-temperature correction to obtain the correct thermal equilibrium at long times.

preprint2011arXiv

Supersymmetric factorization yields exact solutions to the molecular Stark effect problem for "stretched" states

By invoking supersymmetry, we found a condition under which the Stark effect problem for a polar and polarizable molecule subject to nonresonant electric fields becomes exactly solvable. The exact solvability condition for the interaction parameters involved yields exact wavefunction for the "stretched" states, $|J=m,m>$, and for the $|1,1>$ state in the case of a purely induced-dipole interaction. The analytic expressions for the eigenenergy, the space-fixed dipole moment, the alignment cosine, and the expectation value of the angular momentum allow to readily reverse-engineer the problem of finding the values of the interaction parameters required for creating quantum states with preordained characteristics.

preprint2011arXiv

Supersymmetry identifies molecular Stark states whose eigenproperties can be obtained analytically

We made use of supersymmetric (SUSY) quantum mechanics to find a condition under which the Stark effect problem for a polar and polarizable closed-shell diatomic molecule subject to collinear electrostatic and nonresonant radiative fields becomes exactly solvable. The condition, $Δω= \frac{ω^2}{4 (m+1)^2}$, connects values of the dimensionless parameters $ω$ and $Δω$ that characterize the strengths of the permanent and induced dipole interactions of the molecule with the respective fields. The exact solutions are obtained for the $|\tilde{J}=m,m;ω,Δω>$ family of "stretched" states. The field-free and strong-field limits of the combined-fields problem were found to exhibit supersymmetry and shape-invariance, which is indeed the reason why they are analytically solvable. By making use of the analytic form of the $|\tilde{J}=m,m;ω,Δω>$ wavefunctions, we obtained simple formulae for the expectation values of the space-fixed electric dipole moment, the alignment cosine, the angular momentum squared, and derived a "sum rule" which combines the above expectation values into a formula for the eigenenergy. The analytic expressions for the characteristics of the strongly oriented and aligned states provide a direct access to the values of the interaction parameters required for creating such states in the laboratory.

preprint2010arXiv

Entanglement Switch for Dipole Arrays

We propose a new entanglement switch of qubits consisting of electric dipoles, oriented along or against an external electric field and coupled by the electric dipole-dipole interaction. The pairwise entanglement can be tuned and controlled by the ratio of the Rabi frequency and the dipole-dipole coupling strength. Tuning the entanglement can be achieved for one, two and three-dimensional arrangements of the qubits. The feasibility of building such an entanglement switch is also discussed.

preprint2010arXiv

Exact Calculation of Entanglement in a 19-site 2D Spin System

Using the Trace Minimization Algorithm, we carried out an exact calculation of entanglement in a 19-site two-dimensional transverse Ising model. This model consists of a set of localized spin-1/2 particles in a two dimensional triangular lattice coupled through exchange interaction J and subject to an external magnetic field of strength h. We demonstrate, for such a class of two-dimensional magnetic systems, that entanglement can be controlled and tuned by varying the parameter $λ$ = h/J in the Hamiltonian and by introducing impurities into the systems. Examining the derivative of the concurrence as a function of $λ$ shows that the system exhibits a quantum phase transition at about $λ_c$ = 3.01, a transition induced by quantum fluctuations at the absolute zero of temperature.

preprint2010arXiv

Finite Size Scaling for Criticality of the Schrödinger Equation

By solving the Schrödinger equation one obtains the whole energy spectrum, both the bound and the continuum states. If the Hamiltonian depends on a set of parameters, these could be tuned to a transition from bound to continuum states. The behavior of systems near the threshold, which separates bound-states from continuum states, is important in the study of such phenomenon as: ionization of atoms and molecules, molecule dissociation, scattering collisions and stability of matter. In general, the energy is non-analytic as a function of the Hamiltonian parameters or a bound-state does not exist at the threshold energy. The overall goal of this chapter is to show how one can predict, generate and identify new class of stable quantum systems using large-dimensional models and the finite size scaling approach. Within this approach, the finite size corresponds not to the spatial dimension but to the number of elements in a complete basis set used to expand the exact eigenfunction of a given Hamiltonian. This method is efficient and very accurate for estimating the critical parameters, $\{λ_i\}$, for stability of a given Hamiltonian, $H(λ_i)$. We present two methods of obtaining critical parameters using finite size scaling for a given quantum Hamiltonian: The finite element method and the basis set expansion method.

preprint2009arXiv

Simulated Quantum Computation of Global Minima

Finding the optimal solution to a complex optimization problem is of great importance in practically all fields of science, technology, technical design and econometrics. We demonstrate that a modified Grover's quantum algorithm can be applied to real problems of finding a global minimum using modest numbers of quantum bits. Calculations of the global minimum of simple test functions and Lennard-Jones clusters have been carried out on a quantum computer simulator using a modified Grover's algorithm. The number of function evaluations $N$ reduced from O(N) in classical simulation to $O(\sqrt{N})$ in quantum simulation. We also show how the Grover's quantum algorithm can be combined with the classical Pivot method for global optimization to treat larger systems.

preprint2008arXiv

Nuclear-induced time evolution of entanglement of two-electron spins in anisotropically coupled quantum dot

We study the time evolution of entanglement of two spins in anisotropically coupled quantum dot interacting with the unpolarized nuclear spins environment. We assume that the exchange coupling strength in the z-direction $J_z$ is different from the lateral one $J_l$. We observe that the entanglement decays as a result of the coupling to the nuclear environment and reaches a saturation value, which depends on the value of the exchange interaction difference $J=\| J_l-J_z\|$ between the two spins and the strength of the applied external magnetic field. We find that the entanglement exhibits a critical behavior controlled by the competition between the exchange interaction $J$ and the external magnetic field. The entanglement shows a quasi-symmetric behavior above and below a critical value of the exchange interaction. It becomes more symmetric as the external magnetic field increases. The entanglement reaches a large saturation value, close to unity, when the exchange interaction is far above or below its critical value and a small one as it closely approaches the critical value. Furthermore, we find that the decay rate profile of entanglement is linear when the exchange interaction is much higher or lower than the critical value but converts to a power law and finally to a Gaussian as the critical value is approached from both directions. The dynamics of entanglement is found to be independent of the exchange interaction for isotropically coupled quantum dot.

preprint2004arXiv

Scaling of entanglement in finite arrays of exchange-coupled quantum dots

We present a finite-size scaling analysis of the entanglement in a two-dimensional arrays of quantum dots modeled by the Hubbard Hamiltonian on a triangular lattice. Using multistage block renormalization group approach, we have found that there is an abrupt jump of the entanglement when a first-order quantum phase transition occurs. At the critical point, the entanglement is constant, independent of the block size.