Source author record

M. A. Novotny

M. A. Novotny appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

22works
14topics
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

22 published item(s)

preprint2025arXiv

Universal Scaling of Electron Transmission for Nearly Ballistic and Quantum Dragon Nanodevices

We predict two different universal scaling regimes for the quantum transmission of metallic nanodevices following the addition of a small amount of uncorrelated disorder. A nanodevice is connected to two thin semi-infinite uniform leads, and the Non-Equilibrium Green's Function (NEGF) methodology yields the electron transmission ${\cal T}(E)$ as a function of the injected electron energy $E$. Ballistic nanodevices have no disorder and have ${\cal T}(E)=1$ for all $E$ that allow electron propagation in the leads. Quantum dragon nanodevices can have extremely strong properly correlated disorder, and still have ${\cal T}(E)=1$ for all $E$. Additional uncorrelated site disorder leads to Fano resonances in ${\cal T}(E)$. Averaging over the uncorrelated disorder we predict using perturbation theory two universal scaling regimes for ${\cal T}_{\rm ave}(E)$. The functional form of both universal scaling regimes depend on the device length and width, energy, and variance of the uncorrelated disorder. The second scaling regime, valid for small but somewhat larger uncorrelated disorder than the first scaling regime, also has the form dependent on the density of states of the system. These two scaling regimes are demonstrated to be valid via large scale computer calculations.

preprint2022arXiv

Benchmarking quantum annealers using symmetries in embedded subgraphs

We investigate an efficient, generic method for evaluating the performance of quantum annealing devices that does not require the prior knowledge of the true ground state of the benchmark problem. This approach exploits symmetry properties inherent to the ground states of a composite Hamiltonian comprising the benchmark problem Hamiltonian and its symmetric counterpart. Using this method, we compare the performance of two generations of D-Wave machines. Although we do not observe a noticeable difference in the probability of finding solutions with the required symmetry, our results suggest that the current generation of D-Wave machines notably outperforms its predecessor in finding states closer to those with the required symmetry.

preprint2021arXiv

Example Exact Solutions of the Time-independent Gross-Pitaevskii and Schrödinger Equations

A prescription is given to obtain some exact results for certain external potentials $V\left({\vec r}\right)$ of the time-independent Gross-Pitaevskii and Schrödinger equations. The study motivation is the ability to program $V\left({\vec r}\right)$ experimentally in Bose-Einstein condensates. Rather than derive wavefunctions that are solutions for a given $V\left({\vec r}\right)$, we ask which $V\left({\vec r}\right)$ will have a given pdf (probability density function) $P\left({\vec r}\right)$. Several examples in 1D and 2D are presented for well-known pdfs and for the hydrogen atom in momentum space.

preprint2016arXiv

Site and bond percolation thresholds in $K_{n,n}$-based lattices: Vulnerability of quantum annealers to random qubit and coupler failures on chimera topologies

We estimate the critical thresholds of bond and site percolation on nonplanar, effectively two-dimensional graphs with chimera like topology. The building blocks of these graphs are complete and symmetric bipartite subgraphs of size $2n$, referred to as $K_{n,n}$ graphs. For the numerical simulations we use an efficient union-find based algorithm and employ a finite-size scaling analysis to obtain the critical properties for both bond and site percolation. We report the respective percolation thresholds for different sizes of the bipartite subgraph and verify that the associated universality class is that of standard two-dimensional percolation. For the canonical chimera graph used in the D-Wave Systems Inc.~quantum annealer ($n = 4$), we discuss device failure in terms of network vulnerability, i.e., we determine the critical fraction of qubits and couplers that can be absent due to random failures prior to losing large-scale connectivity throughout the device.

preprint2015arXiv

How ubiquitous are dragon segments in quantum transmission?

Quantum dragon segments are nanodevices that have energy-independent total transmission of electrons. At the level of the single-band tight-binding model a nanodevice is viewed as a weighted undirected graph, with a vertex weight given by the on-site energy and the edge weight given by the tight-binding hopping parameter. A quantum dragon is a weighted undirected graph which when connected to idealized semi-infinite input and output leads, has the electron transmission probability ${\cal T}(E)$$=$$1$ for all electron energies $E$. The probability ${\cal T}(E)$ is obtained from the solution of the time-independent Schrödinger equation. A graph must have finely tuned tight-binding parameters in order to have ${\cal T}(E)$$=$$1$. This paper addresses classes of weighted graphs which can be tuned, by adjusting a small fraction of the total weights, to be a quantum dragon. We prove that with proper tuning any nanodevice can be a quantum dragon. Three prescriptions are presented to tune a weighted graph into a quantum dragon nanodevice. The implications of the prescriptions for physical nanodevices is discussed.

preprint2015arXiv

Nonuniversal effects in mixing correlated-growth processes with randomness: Interplay between bulk morphology and surface roughening

To construct continuum stochastic growth equations for competitive nonequilibrium surface-growth processes of the type RD+X that mixes random deposition (RD) with a correlated-growth process X, we use a simplex decomposition of the height field. A distinction between growth processes X that do and do not create voids in the bulk leads to the definition of the {\it effective probability} $p_{\mathrm{eff}}$ of the process X that is a measurable property of the bulk morphology and depends on the {\it activation probability} $p$ of X in the competitive process RD+X. The bulk morphology is reflected in the surface roughening via {\it nonuniversal} prefactors in the universal scaling of the surface width that scales in $p_{\mathrm{eff}}$. The equation and the resulting scaling are derived for X in either a Kardar-Parisi-Zhang or Edwards-Wilkinson universality class in $(1+1)$ dimensions, and illustrated by an example of X being a ballistic deposition. We obtain full data collapse on its corresponding universal scaling function for all $p \in (0;1]$. We outline the generalizations to $(1+n)$ dimensions and to many-component competitive growth processes.

preprint2015arXiv

Quantum Decoherence at Finite Temperatures

We study measures of decoherence and thermalization of a quantum system $S$ in the presence of a quantum environment (bath) $E$. The whole system is prepared in a canonical thermal state at a finite temperature. Applying perturbation theory with respect to the system-environment coupling strength, we find that under common Hamiltonian symmetries, up to first order in the coupling strength it is sufficient to consider the uncoupled system to predict decoherence and thermalization measures of $S$. This decoupling allows closed form expressions for perturbative expansions for the measures of decoherence and thermalization in terms of the free energies of $S$ and of $E$. Numerical results for both coupled and decoupled systems with up to 40 quantum spins validate these findings.

preprint2010arXiv

A new battery-charging method suggested by molecular dynamics simulations

Based on large-scale molecular dynamics simulations, we propose a new charging method that should be capable of charging a Lithium-ion battery in a fraction of the time needed when using traditional methods. This charging method uses an additional applied oscillatory electric field. Our simulation results show that this charging method offers a great reduction in the average intercalation time for Li+ ions, which dominates the charging time. The oscillating field not only increases the diffusion rate of Li+ ions in the electrolyte but, more importantly, also enhances intercalation by lowering the corresponding overall energy barrier.

preprint2010arXiv

A New Charging Method for Li-ion Batteries: Dependence of the charging time on the Direction of an Additional Oscillating Field

We have recently proposed a new method for charging Li-ion batteries based on large-scale molecular dynamics studies (I. Abou Hamad et al, Phys. Chem. Chem. Phys., 12, 2740 (2010)). Applying an additional oscillating electric field in the direction perpendicular to the graphite sheets of the anode showed an exponential decrease in charging time with increasing amplitude of the applied oscillating field. Here we present new results exploring the effect on the charging time of changing the orientation of the oscillating field. Results for oscillating fields in three orthogonal directions are compared.

preprint2010arXiv

Mapping the dynamics of complex multi-dimensional systems onto a discrete set of states conserving mean first passage times: a Projective Dynamics approach

We consider any dynamical system that starts from a given ensemble of configurations and evolves in time until the system reaches a certain fixed stopping criterion, with the mean first-passage time the quantity of interest. We present a general method, Projective Dynamics, which maps the multi-dimensional dynamics of the system onto an arbitrary discrete set of states ${ζ_k}$, subject only to the constraint that the dynamics is restricted to transitions not further than the neighboring states $ζ_{k\pm 1}$. We prove that with this imposed condition there exists a master equation with nearest-neighbor coupling with the same mean first-passage time as the original dynamical system. We show applications of the method for Brownian motion of particles in one and two dimensional potential energy landscapes and the folding process of small bio-polymers. We compare results for the mean first passage time and the mean folding time obtained with the Projective Dynamics method with those obtained by a direct measurement, and where possible with a semi-analytical solution.

preprint2010arXiv

Quantum Transport through Hierarchical Structures

The transport of quantum electrons through hierarchical lattices is of interest because such lattices have some properties of both regular lattices and random systems. We calculate the electron transmission as a function of energy in the tight binding approximation for two related Hanoi networks. HN3 is a Hanoi network with every site having three bonds. HN5 has additional bonds added to HN3 to make the average number of bonds per site equal to five. We present a renormalization group approach to solve the matrix equation involved in this quantum transport calculation. We observe band gaps in HN3, while no such band gaps are observed in linear networks or in HN5.

preprint2009arXiv

Comment on "Dynamic properties in a family of competitive growing models"

The article [Phys. Rev. E {\bf 73}, 031111 (2006)] by Horowitz and Albano reports on simulations of competitive surface-growth models RD+X that combine random deposition (RD) with another deposition X that occurs with probability $p$. The claim is made that at saturation the surface width $w(p)$ obeys a power-law scaling $w(p) \propto 1/p^δ$, where $δ$ is only either $δ=1/2$ or $δ=1$, which is illustrated by the models where X is ballistic deposition and where X is RD with surface relaxation. Another claim is that in the limit $p \to 0^+$, for any lattice size $L$, the time evolution of $w(t)$ generally obeys the scaling $w(p,t) \propto (L^α/p^δ) F(p^{2δ}t/L^z)$, where $F$ is Family-Vicsek universal scaling function. We show that these claims are incorrect.

preprint2007arXiv

Two Modes of Magnetization Switching in a Simulated Iron Nanopillar in an Obliquely Oriented Field

Finite-temperature micromagnetics simulations are employed to study the magnetization-switching dynamics driven by a field applied at an angle to the long axis of an iron nanopillar. A bi-modal distribution in the switching times is observed, and evidence for two competing modes of magnetization-switching dynamics is presented. For the conditions studied here, temperature $T = 20$ K and the reversal field 3160 Oe at an angle of 75$^\circ$ to the long axis, approximately 70% of the switches involve unstable decay (no free-energy barrier) and 30% involve metastable decay (a free-energy barrier is crossed). The latter are indistinguishable from switches which are constrained to start at a metastable free-energy minimum. Competition between unstable and metastable decay could greatly complicate applications involving magnetization switches near the coercive field.

preprint2006arXiv

A rejection-free Monte Carlo method for the hard-disk system

We construct a rejection-free Monte Carlo method for the hard-disk system. Rejection-free Monte Carlo methods preserve the time-evolution behavior of the standard Monte Carlo method, and this relationship is confirmed for our method by observing nonequilibrium relaxation of a bond-orientational order parameter. The rejection-free method gives a greater computational efficiency than the standard method at high densities. The rejection free method is implemented in a shrewd manner using optimization methods to calculate a rejection probability and to update the system. This method should allow an efficient study of the dynamics of two-dimensional solids at high density.

preprint2003arXiv

On the Possibility of Quasi Small-World Nanomaterials

The possibility of materials that are governed by a fixed point related to small world networks is discussed. In particular, large-scale Monte Carlo simulations are performed on Ising ferromagnetic models on two different small-world networks generated from a one-dimensional spin chain. One has the small-world bond strengths independent of the length, and exhibits a finite-temperature phase transition. The other has small-world bonds built from atoms, and although there is no finite-temperature phase transition the system shows a slow power-law change of the effective critical temperature of a finite system as a function of the system size. An outline of a possible synthesis route for quasi small-world nanomaterials is presented.

preprint2001arXiv

A Tutorial on Advanced Dynamic Monte Carlo Methods for Systems with Discrete State Spaces

Advanced algorithms are necessary to obtain faster-than-real-time dynamic simulations in a number of different physical problems that are characterized by widely disparate time scales. Recent advanced dynamic Monte Carlo algorithms that preserve the dynamics of the model are described. These include the $n$-fold way algorithm, the Monte Carlo with Absorbing Markov Chains (MCAMC) algorithm, and the Projective Dynamics (PD) algorithm. To demonstrate the use of these algorithms, they are applied to some simplified models of dynamic physical systems. The models studied include a model for ion motion through a pore such as a biological ion channel and the metastable decay of the ferromagnetic Ising model. Non-trivial parallelization issues for these dynamic algorithms, which are in the class of parallel discrete event simulations, are discussed. Efforts are made to keep the article at an elementary level by concentrating on a simple model in each case that illustrates the use of the advanced dynamic Monte Carlo algorithm.

preprint1999arXiv

Advanced Dynamic Algorithms for the Decay of Metastable Phases in Discrete Spin Models: Bridging Disparate Time Scales

An overview of advanced dynamical algorithms capable of spanning the widely disparate time scales that govern the decay of metastable phases in discrete spin models is presented. The algorithms discussed include constrained transfer-matrix, Monte Carlo with Absorbing Markov Chains (MCAMC), and projective dynamics (PD) methods. The strengths and weaknesses of each of these algorithms are discussed, with particular emphasis on identifying the parameter regimes (system size, temperature, and field) in which each algorithm works best.

preprint1999arXiv

Application of the Projected Dynamics Method to an Anisotropic Heisenberg Model

The Projected Dynamics method was originally developed to study metastable decay in ferromagnetic discrete spin models. Here, we apply it to a classical, continuous Heisenberg model with anisotropic ferromagnetic interactions, which evolves under a Monte Carlo dynamic. The anisotropy is sufficiently large to allow comparison with the Ising model. We describe the Projected Dynamics method and how to apply it to this continuous-spin system. We also discuss how to extract metastable lifetimes and how to extrapolate from small systems to larger systems.

preprint1998arXiv

A projection method for statics and dynamics of lattice spin systems

A method based on Monte Carlo sampling of the probability flows projected onto the subspace of one or more slow variables is proposed for investigation of dynamic and static properties of lattice spin systems. We illustrate the method by applying it, with projection onto the order-parameter subspace, to the three-dimensional 3-state Potts model in equilibrium and to metastable decay in a three-dimensional 3-state kinetic Potts model.

preprint1998arXiv

Test of the Kolmogorov-Johnson-Mehl-Avrami picture of metastable decay in a model with microscopic dynamics

The Kolmogorov-Johnson-Mehl-Avrami (KJMA) theory for the time evolution of the order parameter in systems undergoing first-order phase transformations has been extended by Sekimoto to the level of two-point correlation functions. Here, this extended KJMA theory is applied to a kinetic Ising lattice-gas model, in which the elementary kinetic processes act on microscopic length and time scales. The theoretical framework is used to analyze data from extensive Monte Carlo simulations. The theory is inherently a mesoscopic continuum picture, and in principle it requires a large separation between the microscopic scales and the mesoscopic scales characteristic of the evolving two-phase structure. Nevertheless, we find excellent quantitative agreement with the simulations in a large parameter regime, extending remarkably far towards strong fields (large supersaturations) and correspondingly small nucleation barriers. The original KJMA theory permits direct measurement of the order parameter in the metastable phase, and using the extension to correlation functions one can also perform separate measurements of the nucleation rate and the average velocity of the convoluted interface between the metastable and stable phase regions. The values obtained for all three quantities are verified by other theoretical and computational methods. As these quantities are often difficult to measure directly during a process of phase transformation, data analysis using the extended KJMA theory may provide a useful experimental alternative.

preprint1995arXiv

A Method to Study Relaxation of Metastable Phases: Macroscopic Mean-Field Dynamics

We propose two different macroscopic dynamics to describe the decay of metastable phases in many-particle systems with local interactions. These dynamics depend on the macroscopic order parameter $m$ through the restricted free energy $F(m)$ and are designed to give the correct equilibrium distribution for $m$. The connection between macroscopic dynamics and the underlying microscopic dynamic are considered in the context of a projection- operator formalism. Application to the square-lattice nearest-neighbor Ising ferromagnet gives good agreement with droplet theory and Monte Carlo simulations of the underlying microscopic dynamic. This includes quantitative agreement for the exponential dependence of the lifetime on the inverse of the applied field $H$, and the observation of distinct field regions in which the derivative of the lifetime with respect to $1/H$ depends differently on $H$. In addition, at very low temperatures we observe oscillatory behavior of this derivative with respect to $H$, due to the discreteness of the lattice and in agreement with rigorous results. Similarities and differences between this work and earlier works on finite Ising models in the fixed-magnetization ensemble are discussed.

preprint1994arXiv

Application of a constrained-transfer-matrix method to metastability in the $d$$=$$2$ Ising ferromagnet

Applying a numerical transfer-matrix formalism, we obtain complex-valued constrained free energies for the two-dimensional square-lattice nearest-neighbor Ising ferromagnet below its critical temperature and in an external magnetic field. In particular, we study the imaginary part of the constrained free-energy branch that corresponds to the metastable phase. Although droplets are not introduced explicitly, the metastable free energy is obtained in excellent agreement with field-theoretical droplet-model predictions. The finite-size scaling properties are different in the weak-field and intermediate-field regimes, and we identify the corresponding different critical-droplet shapes. For intermediate fields, we show that the surface free energy of the critical droplet is given by a Wulff construction with the equilibrium surface tension. We also find a prefactor exponent in complete agreement with the field-theoretical droplet model. Our results extend the region of validity for known results of this field-theoretical droplet model, and they indicate that this transfer-matrix approach provides a nonperturbative numerical continuation of the equilibrium free energy into the metastable phase.