BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

O. Melchert

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

ResearcherUnclaimed source record
physics.comp-phcond-mat.dis-nncond-mat.stat-mechBiological Physicsphysics.data-anphysics.opticsData Structures and Algorithmsnlin.CGnlin.PSquant-ph

Catalog footprint

What is connected

23works
10topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

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

23 published item(s)

preprint2022arXiv

A generalized nonlinear Schrödinger Python module implementing different models of input pulse quantum noise

We provide Python tools enabling numerical simulation and analysis of the propagation dynamics of ultrashort laser pulses in nonlinear waveguides. The modeling approach is based on the widely used generalized nonlinear Schrödinger equation for the pulse envelope. The presented software implements the effects of linear dispersion, pulse self-steepening, and the Raman effect. The focus lies on the implementation of input pulse shot noise, i.e. classical background fields that mimick quantum noise, which are often not thoroughly presented in the scientific literature. We discuss and implement commonly adopted quantum noise models based on pure spectral phase noise, as well as Gaussian noise. Coherence properties of the resulting spectra can be calculated. We demonstrate the functionality of the software by reproducing results for a supercontinuum generation process in a photonic crystal fiber, documented in the scientific literature. The presented Python tools are are open-source and released under the MIT license in a publicly available software repository.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Resonant Kushi-comb-like multi-frequency radiation of oscillating two-color soliton molecules

Nonlinear waveguides with two distinct domains of anomalous dispersion can support the formation of molecule-like two-color pulse compounds. They consist of two tightly bound subpulses with frequency loci separated by a vast frequency gap. Perturbing such a two-color pulse compound triggers periodic amplitude and width variations, reminiscent of molecular vibrations. With increasing strength of perturbation, the dynamics of the pulse compound changes from harmonic to nonlinear oscillations. The periodic amplitude variations enable coupling of the pulse compound to dispersive waves, resulting in the resonant emission of multi-frequency radiation. We demonstrate that the location of the resonances can be precisely predicted by phase-matching conditions. If the pulse compound consists of a pair of identical subpulses, inherent symmetries lead to degeneracies in the resonance spectrum. Weak perturbations lift existing degeneracies and cause a splitting of the resonance lines into multiple lines. Strong perturbations result in more complex emission spectra, characterized by well separated spectral bands caused by resonant Cherenkov radiation and additional four-wave mixing processes.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

A python package for ultrashort optical pulse propagation in terms of forward models for the analytic signal

We present a flexible, open-source Python package for the accurate simulation of the $z$-propagation dynamics of ultrashort optical pulses in nonlinear waveguides, especially valid for few-cycle pulses and their interaction. The simulation approach is based on unidirectional propagation equations for the analytic signal. The provided software allows to account for dispersion, attenuation, four-wave mixing processes including, e.g., third-harmonic generation, and features various models for the Raman response. The propagation equations are solved on a periodic temporal domain. For $z$-propagation, a selection of pseudospectral methods is available. Propagation scenarios for a custom propagation constant and initial field pulses can either be specified in terms of a HDF5 based input file format or by direct implementation using a python script. We demonstrate the functionality for a test-case for which an exact solution is available, by reproducing exemplary results documented in the scientific literature, and a complex propagation scenario involving multiple pulses. The py-fmas code, its reference manual, an extended user guide, and further usage examples are available online at https://github.com/omelchert/py-fmas.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Detection, numerical simulation and approximate inversion of optoacoustic signals generated in multi-layered PVA hydrogel based tissue phantoms

In this article we characterize optoacoustic signals generated from layered tissue phantoms via short laser pulses by experimental and numerical means. In particular, we consider the case where scattering is effectively negligible and the absorbed energy density follows Beer-Lambert's law, i.e. is characterized by an exponential decay within the layers and discontinuities at interfaces. We complement experiments on samples with multiple layers, where the material properties are known a priori, with numerical calculations for a pointlike detector, tailored to suit our experimental setup. Experimentally, we characterize the acoustic signal observed by a piezoelectric detector in the acoustic far-field in backward mode and we discuss the implication of acoustic diffraction on our measurements. We further attempt an inversion of an OA signal in the far-field approximation.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Efficient polar convolution based on the discrete Fourier-Bessel transform for application in computational biophotonics

We discuss efficient algorithms for the accurate forward and reverse evaluation of the discrete Fourier-Bessel transform (dFBT) as numerical tools to assist in the 2D polar convolution of two radially symmetric functions, relevant, e.g., to applications in computational biophotonics. In our survey of the numerical procedure we account for the circumstance that the objective function might result from a more complex measurement process and is, in the worst case, known on a finite sequence of coordinate values, only. We contrast the performance of the resulting algorithms with a procedure based on a straight forward numerical quadrature of the underlying integral transform and asses its efficienty for two benchmark Fourier-Bessel pairs. An application to the problem of finite-size beam-shape convolution in polar coordinates, relevant in the context of tissue optics and optoacoustics, is used to illustrate the versatility and computational efficiency of the numerical procedure.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Optoacoustic inversion via Volterra kernel reconstruction

In this letter we address the numeric inversion of optoacoustic signals to initial stress profiles. Therefore we put under scrutiny the optoacoustic kernel reconstruction problem in the paraxial approximation of the underlying wave-equation. We apply a Fourier-series expansion of the optoacoustic Volterra kernel and obtain the respective expansion coefficients for a given "apparative" setup by performing a gauge procedure using synthetic input data. The resulting effective kernel is subsequently used to solve the optoacoustic source reconstruction problem for general signals. We verify the validity of the proposed inversion protocol for synthetic signals and explore the feasibility of our approach to also account for the diffraction transformation of signals beyond the paraxial approximation.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Analysis of the phase transition in the $2D$ Ising ferromagnet using a Lempel-Ziv string parsing scheme and black-box data-compression utilities

In this work we consider information-theoretical observables to analyze short symbolic sequences, comprising time-series that represent the orientation of a single spin in a $2D$ Ising ferromagnet on a square lattice of size $L^2=128^2$, for different system temperatures $T$. The latter were chosen from an interval enclosing the critical point $T_{\rm c}$ of the model. At small temperatures the sequences are thus very regular, at high temperatures they are maximally random. In the vicinity of the critical point, nontrivial, long-range correlations appear. Here, we implement estimators for the entropy rate, excess entropy (i.e. "complexity") and multi-information. First, we implement a Lempel-Ziv string parsing scheme, providing seemingly elaborate entropy rate and multi-information estimates and an approximate estimator for the excess entropy. Furthermore, we apply easy-to-use black-box data compression utilities, providing approximate estimators only. For comparison and to yield results for benchmarking purposes we implement the information-theoretic observables also based on the well-established M-block Shannon entropy, which is more tedious to apply compared to the the first two "algorithmic" entropy estimation procedures. To test how well one can exploit the potential of such data compression techniques, we aim at detecting the critical point of the $2D$ Ising ferromagnet. Among the above observables, the multi-information, which is known to exhibit an isolated peak at the critical point, is very easy to replicate by means of both efficient algorithmic entropy estimation procedures. Finally, we assess how good the various algorithmic entropy estimates compare to the more conventional block entropy estimates and illustrate a simple modification that yields enhanced results.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Biased and greedy random walks on two-dimensional lattices with quenched randomness: the "greedy" ant within a disordered environment

The principle characteristics of biased greedy random walks (BGRWs) on two-dimensional lattices with real-valued quenched disorder on the lattice edges are studied. Here, the disorder allows for negative edge-weights. In previous studies, considering the negative-weight percolation (NWP) problem, this was shown to change the universality class of the existing, static percolation transition. In the presented study, four different types of BGRWs and an algorithm based on the ant colony optimization (ACO) heuristic were considered. Regarding the BGRWs, the precise configurations of the lattice walks constructed during the numerical simulations were influenced by two parameters: a disorder parameter rho that controls the amount of negative edge weights on the lattice and a bias strength B that governs the drift of the walkers along a certain lattice direction. Here, the pivotal observable is the probability that, after termination, a lattice walk exhibits a total negative weight, which is here considered as percolating. The behavior of this observable as function of rho for different bias strengths B is put under scrutiny. Upon tuning rho, the probability to find such a feasible lattice walk increases from zero to one. This is the key feature of the percolation transition in the NWP model. Here, we address the question how well the transition point rho_c, resulting from numerically exact and "static" simulations in terms of the NWP model can be resolved using simple dynamic algorithms that have only local information available, one of the basic questions in the physics of glassy systems.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Percolation thresholds on planar Euclidean relative neighborhood graphs

In the presented article, statistical properties regarding the topology and standard percolation on relative neighborhood graphs (RNGs) for planar sets of points, considering the Euclidean metric, are put under scrutiny. RNGs belong to the family of ``proximity graphs'', i.e. their edge-set encodes proximity information regarding the close neighbors for the terminal nodes of a given edge. Therefore they are, e.g., discussed in the context of the construction of backbones for wireless ad-hoc networks that guarantee connectedness of all underlying nodes. Here, by means of numerical simulations, we determine the asymptotic degree and diameter of RNGs and we estimate their bond and site percolation thresholds, which were previously conjectured to be nontrivial. We compare the results to regular 2D graphs for which the degree is close to that of the RNG. Finally, we deduce the common percolation critical exponents from the RNG data to verify that the associated universality class is that of standard 2D percolation.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Typical and large-deviation properties of minimum-energy paths on disordered hierarchical lattices

We perform numerical simulations to study the optimal path problem on disordered hierarchical graphs with effective dimension d=2.32. Therein, edge energies are drawn from a disorder distribution that allows for positive and negative energies. This induces a behavior which is fundamentally different from the case where all energies are positive, only. Upon changing the subtleties of the distribution, the scaling of the minimum energy path length exhibits a transition from self-affine to self-similar. We analyze the precise scaling of the path length and the associated ground-state energy fluctuations in the vincinity of the disorder critical point, using a decimation procedure for huge graphs. Further, using an importance sampling procedure in the disorder we compute the negative-energy tails of the ground-state energy distribution up to 12 standard deviations away from its mean. We find that the asymptotic behavior of the negative-energy tail is in agreement with a Tracy-Widom distribution. Further, the characteristic scaling of the tail can be related to the ground-state energy flucutations, similar as for the directed polymer in a random medium.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

A computational mechanics approach to estimate entropy and (approximate) complexity for the dynamics of the 2D Ising Ferromagnet

We present a numerical analysis of the entropy rate and statistical complexity related to the spin flip dynamics of the 2D Ising Ferromagnet at different temperatures T. We follow an information theoretic approach and test three different entropy estimation algorithms to asses entropy rate and statistical complexity of binary sequences. The latter are obtained by monitoring the orientation of a single spin on a square lattice of side-length L=256 at a given temperature parameter over time. The different entropy estimation procedures are based on the M-block Shannon entropy (a well established method that yields results for benchmarking purposes), non-sequential recursive pair substitution (providing an elaborate and an approximate estimator) and a convenient data compression algorithm contained in the zlib-library (providing an approximate estimator only). We propose an approximate measure of statistical complexity that emphasizes on correlations within the sequence and which is easy to implement, even by means of black-box data compression algorithms. Regarding the 2D Ising Ferromagnet simulated using Metropolis dynamics and for binary sequences of finite length, the proposed approximate complexity measure is peaked close to the critical temperature. For the approximate estimators, a finite-size scaling analysis reveals that the peak approaches the critical temperature as the sequence length increases. Results obtained using different spin-flip dynamics are briefly discussed. The suggested complexity measure can be extended to non-binary sequences in a straightforward manner.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Analysis of the loop length distribution for the negative weight percolation problem in dimensions d=2 through 6

We consider the negative weight percolation (NWP) problem on hypercubic lattice graphs with fully periodic boundary conditions in all relevant dimensions from d=2 to the upper critical dimension d=6. The problem exhibits edge weights drawn from disorder distributions that allow for weights of either sign. We are interested in in the full ensemble of loops with negative weight, i.e. non-trivial (system spanning) loops as well as topologically trivial ("small") loops. The NWP phenomenon refers to the disorder driven proliferation of system spanning loops of total negative weight. While previous studies where focused on the latter loops, we here put under scrutiny the ensemble of small loops. Our aim is to characterize -using this extensive and exhaustive numerical study- the loop length distribution of the small loops right at and below the critical point of the hypercubic setups by means of two independent critical exponents. These can further be related to the results of previous finite-size scaling analyses carried out for the system spanning loops. For the numerical simulations we employed a mapping of the NWP model to a combinatorial optimization problem that can be solved exactly by using sophisticated matching algorithms. This allowed us to study here numerically exact very large systems with high statistics.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Basic Data Analysis and More - A Guided Tour Using Python

In these lecture notes, a selection of frequently required statistical tools will be introduced and illustrated. They allow to post-process data that stem from, e.g., large-scale numerical simulations (aka sequence of random experiments). From a point of view of data analysis, the concepts and techniques introduced here are of general interest and are, at best, employed by computational aid. Consequently, an exemplary implementation of the presented techniques using the Python programming language is provided. The contents of these lecture notes is rather selective and represents a computational experimentalist's view on the subject of basic data analysis, ranging from the simple computation of moments for distributions of random variables to more involved topics such as hierarchical cluster analysis and the parallelization of Python code.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Encryption dynamics and avalanche parameter for "delayed dynamics"-based cryptosystems

The presented article attempts to characterize the encryption dynamics of delayed dynamics based block ciphers, designed for the encryption of binary data. For such encryption algorithms, the encryption process relies on a coupling dynamics with time delay between different bits in the plaintext (i.e.\ the "initial" message to be encrypted). Here, the principal dynamics of the encryption process is examined and the Hammingdistance is used to quantify the change in ciphertext (i.e.\ the plaintext after encryption) upon changing a single bit in the plaintext message or slightly perturbing the key used during encryption. More precisely, the previously proposed "encryption via delayed dynamics" (in short: EDDy) algorithm as well as its extended version (termed ExEDDy) are analyzed by means of numerical simulations. As a result it is found that while EDDy exhibits a rather poor perfomance, ExEDDy performes considerably better and hence constitutes a significant improvement over EDDy. Consequently, the results are contrasted with those obtained for a block cipher that implements the encryption/decryption dynamics by means of reversible cellular automata.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Information theoretic approach to ground-state phase transitions for two and three-dimensional frustrated spin systems

The information theoretic observables entropy (a measure of disorder), excess entropy (a measure of complexity) and multi information are used to analyze ground-state spin configurations for disordered and frustrated model systems in 2D and 3D. For both model systems, ground-state spin configurations can be obtained in polynomial time via exact combinatorial optimization algorithms, which allowed us to study large systems with high numerical accuracy. Both model systems exhibit a continuous transition from an ordered to a disordered ground state as a model parameter is varied. By using the above information theoretic observables it is possible to detect changes in the spatial structure of the ground states as the critical point is approached. It is further possible to quantify the scaling behavior of the information theoretic observables in the vicinity of the critical point. For both model systems considered, the estimates of critical properties for the ground-state phase transitions are in good agreement with existing results reported in the literature.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Is negative-weight percolation compatible with SLE?

We study numerically the geometrical properties of minimally weighted paths that appear in the negative-weight percolation (NWP) model on two-dimensional lattices assuming a combination of periodic and free boundary conditions (BCs). Each realization of the disorder consists of a random fraction 1-rho of bonds with unit strength and a fraction rho of bond strengths drawn from a Gaussian distribution with zero mean and unit width. For each such sample, the path is forced to span the lattice along the direction with the free BCs. The path and a set of negatively weighted loops form a ground state (GS). A ground state on such a lattice can be determined performing a non-trivial transformation of the original graph and applying sophisticated matching algorithms. Here we examine whether the geometrical properties of the paths are in accordance with predictions of Schramm-Loewner evolution (SLE). Measuring the fractal dimension and reviewing Schramm's left passage formula indicates that the paths cannot be described in terms of SLE.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Minimum weight spanning trees of weighted scale free networks

In this lecture we will consider the minimum weight spanning tree (MST) problem, i.e., one of the simplest and most vital combinatorial optimization problems. We will discuss a particular greedy algorithm that allows to compute a MST for undirected weighted graphs, namely Kruskal's algorithm, and we will study the structure of MSTs obtained for weighted scale free random graphs. This is meant to clarify whether the structure of MSTs is sensitive to correlations between edge weights and topology of the underlying scale free graphs.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

The universality class of the continuous phase transition in the 2D "Touch and Stop" cluster growth percolation model

We consider the "Touch and Stop" cluster growth percolation (CGP) model on the two dimensional square lattice. A key-parameter in the model is the fraction p of occupied "seed" sites that act as nucleation centers from which a particular cluster growth procedure is started. Here, we consider two growth-styles: rhombic and disk-shaped cluster growth. For intermediate values of p the final state, attained by the growth procedure, exhibits a cluster of occupied sites that spans the entire lattice. Using numerical simulations we investigate the percolation probability and the order parameter and perform a finite-size scaling analysis for lattices of side length up to L=1024 in order to carefully determine the critical exponents that govern the respective transition. In contrast to previous studies, reported in [Tsakiris et al., Phys. Rev. E 82 (2010) 041108], we find strong numerical evidence that the CGP model is in the standard percolation universality class.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Mean-field behavior of the negative-weight percolation model on random regular graphs

We investigate both analytically and numerically the ensemble of minimum-weight loops and paths in the negative-weight percolation model on random graphs with fixed connectivity and bimodal weight distribution. This allows us to study the mean-field behavior of this model. The analytical study is based on a conjectured equivalence with the problem of self-avoiding walks in a random medium. The numerical study is based on a mapping to a standard minimum-weight matching problem for which fast algorithms exist. Both approaches yield results which are in agreement, on the location of the phase transition, on the value of critical exponents, and on the absence of any sizeable indications of a glass phase. By these results, the previously conjectured upper critical dimension of d_u=6 is confirmed.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

A dedicated algorithm for calculating ground states for the triangular random bond Ising model

In the presented article we present an algorithm for the computation of ground state spin configurations for the 2d random bond Ising model on planar triangular lattice graphs. Therefore, it is explained how the respective ground state problem can be mapped to an auxiliary minimum-weight perfect matching problem, solvable in polynomial time. Consequently, the ground state properties as well as minimum-energy domain wall (MEDW) excitations for very large 2d systems, e.g. lattice graphs with up to N=384x384 spins, can be analyzed very fast. Here, we investigate the critical behavior of the corresponding T=0 ferromagnet to spin-glass transition, signaled by a breakdown of the magnetization, using finite-size scaling analyses of the magnetization and MEDW excitation energy and we contrast our numerical results with previous simulations and presumably exact results.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Configurational statistics of densely and fully packed loops in the negative-weight percolation model

By means of numerical simulations we investigate the configurational properties of densely and fully packed configurations of loops in the negative-weight percolation (NWP) model. In the presented study we consider 2d square, 2d honeycomb, 3d simple cubic and 4d hypercubic lattice graphs, where edge weights are drawn from a Gaussian distribution. For a given realization of the disorder we then compute a configuration of loops, such that the configurational energy, given by the sum of all individual loop weights, is minimized. For this purpose, we employ a mapping of the NWP model to the "minimum-weight perfect matching problem" that can be solved exactly by using sophisticated polynomial-time matching algorithms. We characterize the loops via observables similar to those used in percolation studies and perform finite-size scaling analyses, up to side length L=256 in 2d, L=48 in 3d and L=20 in 4d (for which we study only some observables), in order to estimate geometric exponents that characterize the configurations of densely and fully packed loops. One major result is that the loops behave like uncorrelated random walks from dimension d=3 on, in contrast to the previously studied behavior at the percolation threshold, where random-walk behavior is obtained for d>=6.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

The upper critical dimension of the negative-weight percolation problem

By means of numerical simulations we investigate the geometric properties of loops on hypercubic lattice graphs in dimensions d=2 through 7, where edge weights are drawn from a distribution that allows for positive and negative weights. We are interested in the appearance of system-spanning loops of total negative weight. The resulting negative-weight percolation (NWP) problem is fundamentally different from conventional percolation, as we have seen in previous studies of this model for the 2d case. Here, we characterize the transition for hypercubic systems, where the aim of the present study is to get a grip on the upper critical dimension d_u of the NWP problem. For the numerical simulations we employ a mapping of the NWP model to a combinatorial optimization problem that can be solved exactly by using sophisticated matching algorithms. We characterize the loops via observables similar to those in percolation theory and perform finite-size scaling analyses, e.g. 3d hypercubic systems with side length up to L=56 sites, in order to estimate the critical properties of the NWP phenomenon. We find our numerical results consistent with an upper critical dimension d_u=6 for the NWP problem.

0 citations0 reviewsNo code linkedOpen access