Showing 1–22 of 22 results for author: Knudsen, J

Engineering two-dimensional materials from single-layer NbS$_2$

Authors: Timo Knispel, Daniela Mohrenstecher, Carsten Speckmann, Affan Safeer, Camiel van Efferen, Virgínia Boix, Alexander Grüneis, Wouter Jolie, Alexei Preobrajenski, Jan Knudsen, Nicolae Atodiresei, Thomas Michely, Jeison Fischer

Abstract: Starting from a single layer of NbS$_2$ grown on graphene by molecular beam epitaxy, the single unit cell thick 2D materials Nb$_{5/3}$S$_3$-2D and Nb$_2$S$_3$-2D are created using two different pathways. Either annealing under sulfur-deficient conditions at progressively higher temperatures or deposition of increasing amounts of Nb at elevated temperature result in phase-pure Nb$_{5/3}$S$_3$-2D f… ▽ More Starting from a single layer of NbS$_2$ grown on graphene by molecular beam epitaxy, the single unit cell thick 2D materials Nb$_{5/3}$S$_3$-2D and Nb$_2$S$_3$-2D are created using two different pathways. Either annealing under sulfur-deficient conditions at progressively higher temperatures or deposition of increasing amounts of Nb at elevated temperature result in phase-pure Nb$_{5/3}$S$_3$-2D followed by Nb$_2$S$_3$-2D. The materials are characterized by scanning tunneling microscopy, scanning tunneling spectroscopy and X-ray photoemission spectroscopy. The experimental assessment combined with systematic density functional theory calculations reveals their structure. The 2D materials are covalently bound without any van der Waals gap. Their stacking sequence and structure are at variance with expectations based on corresponding bulk materials highlighting the importance of surface and interface effects in structure formation. △ Less

Submitted 24 July, 2024; originally announced July 2024.

Novel 2D vanadium sulphides: synthesis, atomic structure engineering and charge density waves

Authors: Camiel van Efferen, Joshua Hall, Nicolae Atodiresei, Virgínia Boix, Affan Safeer, Tobias Wekking, Nikolay Vinogradov, Alexei Preobrajensk, Jan Knudsen, Jeison Fischer, Wouter Jolie, Thomas Michely

Abstract: Two new ultimately thin vanadium rich 2D materials based on VS2 are created via molecular beam epitaxy and investigated using scanning tunneling microscopy, X-ray photoemission spectroscopy and density-functional theory calculations. The controlled synthesis of stoichiometric single-layer VS2 or either of the two vanadium-rich materials is achieved by varying the sample coverage and the sulphur pr… ▽ More Two new ultimately thin vanadium rich 2D materials based on VS2 are created via molecular beam epitaxy and investigated using scanning tunneling microscopy, X-ray photoemission spectroscopy and density-functional theory calculations. The controlled synthesis of stoichiometric single-layer VS2 or either of the two vanadium-rich materials is achieved by varying the sample coverage and the sulphur pressure during annealing. Through annealing of small stoichiometric single-layer VS2 islands without S pressure, S-vacancies spontaneously order in 1D arrays, giving rise to patterned adsorption. Via the comparison of density-functional theory calculations with scanning tunneling microscopy data, the atomic structure of the S-depleted phase, with a stoichiometry of V4S7, is determined. By depositing larger amounts of vanadium and sulphur, which are subsequently annealed in a S-rich atmosphere, self-intercalated ultimately thin V5S8-derived layers are obtained, which host 2x2 V-layers between sheets of VS2. We provide atomic models for the thinnest V5S8-derived structures. Finally, we use scanning tunneling spectroscopy to investigate the charge density wave observed in the 2D V5S8-derived islands. △ Less

Submitted 13 February, 2024; v1 submitted 29 June, 2023; originally announced June 2023.

The Diffuse Auroral Eraser

Authors: R. N. Troyer, A. N. Jaynes, S. L. Jones, D. J. Knudsen, T. S. Trondsen

Abstract: The source of diffuse aurora has been widely studied and linked to electron cyclotron harmonic (ECH) and upper-band chorus (UBC) waves. It is known that these waves scatter 100s of eV to 10s of keV electrons from the plasma sheet, but the relative contribution of each wave type is still an open question. In this paper, we report on an interesting and unusual auroral feature observed on March 15, 2… ▽ More The source of diffuse aurora has been widely studied and linked to electron cyclotron harmonic (ECH) and upper-band chorus (UBC) waves. It is known that these waves scatter 100s of eV to 10s of keV electrons from the plasma sheet, but the relative contribution of each wave type is still an open question. In this paper, we report on an interesting and unusual auroral feature observed on March 15, 2002. We believe that these observations could help further our understanding of waves associated with diffuse aurora. This diffuse auroral feature is characterized by four phases: (1) the initial phase exhibiting regular diffuse aurora, (2) the brightening phase, where an east-west auroral stripe rapidly brightens, (3) the eraser phase, where the stripe dims to below its initial state, and (4) the recovery phase, where the diffuse aurora returns to its original brightness. Using a superposed epoch analysis of 22 events, we calculate the average recovery phase time to be 20 seconds, although this varies widely between events. We hypothesize that the process responsible for these auroral eraser events could be the result of chorus waves modulating diffuse aurora. △ Less

Submitted 11 May, 2021; originally announced May 2021.

Comments: 19 pages, 8 figures

Journal ref: Journal of Geophysical Research: Space Physics, 126, e2020JA028805

Load Balancing with Dynamic Set of Balls and Bins

Authors: Anders Aamand, Jakob Bæk Tejs Knudsen, Mikkel Thorup

Abstract: In dynamic load balancing, we wish to distribute balls into bins in an environment where both balls and bins can be added and removed. We want to minimize the maximum load of any bin but we also want to minimize the number of balls and bins affected when adding or removing a ball or a bin. We want a hashing-style solution where we given the ID of a ball can find its bin efficiently. We are given… ▽ More In dynamic load balancing, we wish to distribute balls into bins in an environment where both balls and bins can be added and removed. We want to minimize the maximum load of any bin but we also want to minimize the number of balls and bins affected when adding or removing a ball or a bin. We want a hashing-style solution where we given the ID of a ball can find its bin efficiently. We are given a balancing parameter $c=1+ε$, where $ε\in (0,1)$. With $n$ and $m$ the current numbers of balls and bins, we want no bin with load above $C=\lceil c n/m\rceil$, referred to as the capacity of the bins. We present a scheme where we can locate a ball checking $1+O(\log 1/ε)$ bins in expectation. When inserting or deleting a ball, we expect to move $O(1/ε)$ balls, and when inserting or deleting a bin, we expect to move $O(C/ε)$ balls. Previous bounds were off by a factor $1/ε$. These bounds are best possible when $C=O(1)$ but for larger $C$, we can do much better: Let $f=εC$ if $C\leq \log 1/ε$, $f=ε\sqrt{C}\cdot \sqrt{\log(1/(ε\sqrt{C}))}$ if $\log 1/ε\leq C<\tfrac{1}{2ε^2}$, and $C=1$ if $C\geq \tfrac{1}{2ε^2}$. We show that we expect to move $O(1/f)$ balls when inserting or deleting a ball, and $O(C/f)$ balls when inserting or deleting a bin. For the bounds with larger $C$, we first have to resolve a much simpler probabilistic problem. Place $n$ balls in $m$ bins of capacity $C$, one ball at the time. Each ball picks a uniformly random non-full bin. We show that in expectation and with high probability, the fraction of non-full bins is $Θ(f)$. Then the expected number of bins that a new ball would have to visit to find one that is not full is $Θ(1/f)$. As it turns out, we obtain the same complexity in our more complicated scheme where both balls and bins can be added and removed. △ Less

Submitted 11 April, 2021; originally announced April 2021.

Comments: Accepted at STOC'21

On Sums of Monotone Random Integer Variables

Authors: Anders Aamand, Noga Alon, Jakob Bæk Tejs Knudsen, Mikkel Thorup

Abstract: We say that a random integer variable $X$ is monotone if the modulus of the characteristic function of $X$ is decreasing on $[0,π]$. This is the case for many commonly encountered variables, e.g., Bernoulli, Poisson and geometric random variables. In this note, we provide estimates for the probability that the sum of independent monotone integer variables attains precisely a specific value. We do… ▽ More We say that a random integer variable $X$ is monotone if the modulus of the characteristic function of $X$ is decreasing on $[0,π]$. This is the case for many commonly encountered variables, e.g., Bernoulli, Poisson and geometric random variables. In this note, we provide estimates for the probability that the sum of independent monotone integer variables attains precisely a specific value. We do not assume that the variables are identically distributed. Our estimates are sharp when the specific value is close to the mean, but they are not useful further out in the tail. By combining with the trick of \emph{exponential tilting}, we obtain sharp estimates for the point probabilities in the tail under a slightly stronger assumption on the random integer variables which we call strong monotonicity. △ Less

Submitted 13 April, 2021; v1 submitted 8 April, 2021; originally announced April 2021.

Comments: 8 pages

The Power of Hashing with Mersenne Primes

Authors: Thomas Dybdahl Ahle, Jakob Tejs Bæk Knudsen, Mikkel Thorup

Abstract: The classic way of computing a $k$-universal hash function is to use a random degree-$(k-1)$ polynomial over a prime field $\mathbb Z_p$. For a fast computation of the polynomial, the prime $p$ is often chosen as a Mersenne prime $p=2^b-1$. In this paper, we show that there are other nice advantages to using Mersenne primes. Our view is that the hash function's output is a $b$-bit integer that i… ▽ More The classic way of computing a $k$-universal hash function is to use a random degree-$(k-1)$ polynomial over a prime field $\mathbb Z_p$. For a fast computation of the polynomial, the prime $p$ is often chosen as a Mersenne prime $p=2^b-1$. In this paper, we show that there are other nice advantages to using Mersenne primes. Our view is that the hash function's output is a $b$-bit integer that is uniformly distributed in $\{0, \dots, 2^b-1\}$, except that $p$ (the all \texttt1s value in binary) is missing. Uniform bit strings have many nice properties, such as splitting into substrings which gives us two or more hash functions for the cost of one, while preserving strong theoretical qualities. We call this trick "Two for one" hashing, and we demonstrate it on 4-universal hashing in the classic Count Sketch algorithm for second-moment estimation. We also provide a new fast branch-free code for division and modulus with Mersenne primes. Contrasting our analytic work, this code generalizes to any Pseudo-Mersenne primes $p=2^b-c$ for small $c$. △ Less

Submitted 6 May, 2021; v1 submitted 19 August, 2020; originally announced August 2020.

No Repetition: Fast Streaming with Highly Concentrated Hashing

Authors: Anders Aamand, Debarati Das, Evangelos Kipouridis, Jakob B. T. Knudsen, Peter M. R. Rasmussen, Mikkel Thorup

Abstract: To get estimators that work within a certain error bound with high probability, a common strategy is to design one that works with constant probability, and then boost the probability using independent repetitions. Important examples of this approach are small space algorithms for estimating the number of distinct elements in a stream, or estimating the set similarity between large sets. Using sta… ▽ More To get estimators that work within a certain error bound with high probability, a common strategy is to design one that works with constant probability, and then boost the probability using independent repetitions. Important examples of this approach are small space algorithms for estimating the number of distinct elements in a stream, or estimating the set similarity between large sets. Using standard strongly universal hashing to process each element, we get a sketch based estimator where the probability of a too large error is, say, 1/4. By performing $r$ independent repetitions and taking the median of the estimators, the error probability falls exponentially in $r$. However, running $r$ independent experiments increases the processing time by a factor $r$. Here we make the point that if we have a hash function with strong concentration bounds, then we get the same high probability bounds without any need for repetitions. Instead of $r$ independent sketches, we have a single sketch that is $r$ times bigger, so the total space is the same. However, we only apply a single hash function, so we save a factor $r$ in time, and the overall algorithms just get simpler. Fast practical hash functions with strong concentration bounds were recently proposed by Aamand em et al. (to appear in STOC 2020). Using their hashing schemes, the algorithms thus become very fast and practical, suitable for online processing of high volume data streams. △ Less

Submitted 2 April, 2020; originally announced April 2020.

Comments: 10 pages

Almost Optimal Tensor Sketch

Authors: Thomas D. Ahle, Jakob B. T. Knudsen

Abstract: We construct a matrix $M\in R^{m\otimes d^c}$ with just $m=O(c\,λ\,\varepsilon^{-2}\text{poly}\log1/\varepsilonδ)$ rows, which preserves the norm $\|Mx\|_2=(1\pm\varepsilon)\|x\|_2$ of all $x$ in any given $λ$ dimensional subspace of $ R^d$ with probability at least $1-δ$. This matrix can be applied to tensors $x^{(1)}\otimes\dots\otimes x^{(c)}\in R^{d^c}$ in $O(c\, m \min\{d,m\})$ time -- hence… ▽ More We construct a matrix $M\in R^{m\otimes d^c}$ with just $m=O(c\,λ\,\varepsilon^{-2}\text{poly}\log1/\varepsilonδ)$ rows, which preserves the norm $\|Mx\|_2=(1\pm\varepsilon)\|x\|_2$ of all $x$ in any given $λ$ dimensional subspace of $ R^d$ with probability at least $1-δ$. This matrix can be applied to tensors $x^{(1)}\otimes\dots\otimes x^{(c)}\in R^{d^c}$ in $O(c\, m \min\{d,m\})$ time -- hence the name "Tensor Sketch". (Here $x\otimes y = \text{asvec}(xy^T) = [x_1y_1, x_1y_2,\dots,x_1y_m,x_2y_1,\dots,x_ny_m]\in R^{nm}$.) This improves upon earlier Tensor Sketch constructions by Pagh and Pham~[TOCT 2013, SIGKDD 2013] and Avron et al.~[NIPS 2014] which require $m=Ω(3^cλ^2δ^{-1})$ rows for the same guarantees. The factors of $λ$, $\varepsilon^{-2}$ and $\log1/δ$ can all be shown to be necessary making our sketch optimal up to log factors. With another construction we get $λ$ times more rows $m=\tilde O(c\,λ^2\,\varepsilon^{-2}(\log1/δ)^3)$, but the matrix can be applied to any vector $x^{(1)}\otimes\dots\otimes x^{(c)}\in R^{d^c}$ in just $\tilde O(c\, (d+m))$ time. This matches the application time of Tensor Sketch while still improving the exponential dependencies in $c$ and $\log1/δ$. Technically, we show two main lemmas: (1) For many Johnson Lindenstrauss (JL) constructions, if $Q,Q'\in R^{m\times d}$ are independent JL matrices, the element-wise product $Qx \circ Q'y$ equals $M(x\otimes y)$ for some $M\in R^{m\times d^2}$ which is itself a JL matrix. (2) If $M^{(i)}\in R^{m\times md}$ are independent JL matrices, then $M^{(1)}(x \otimes (M^{(2)}y \otimes \dots)) = M(x\otimes y\otimes \dots)$ for some $M\in R^{m\times d^c}$ which is itself a JL matrix. Combining these two results give an efficient sketch for tensors of any size. △ Less

Submitted 3 September, 2019; originally announced September 2019.

Oblivious Sketching of High-Degree Polynomial Kernels

Authors: Thomas D. Ahle, Michael Kapralov, Jakob B. T. Knudsen, Rasmus Pagh, Ameya Velingker, David Woodruff, Amir Zandieh

Abstract: Kernel methods are fundamental tools in machine learning that allow detection of non-linear dependencies between data without explicitly constructing feature vectors in high dimensional spaces. A major disadvantage of kernel methods is their poor scalability: primitives such as kernel PCA or kernel ridge regression generally take prohibitively large quadratic space and (at least) quadratic time, a… ▽ More Kernel methods are fundamental tools in machine learning that allow detection of non-linear dependencies between data without explicitly constructing feature vectors in high dimensional spaces. A major disadvantage of kernel methods is their poor scalability: primitives such as kernel PCA or kernel ridge regression generally take prohibitively large quadratic space and (at least) quadratic time, as kernel matrices are usually dense. Some methods for speeding up kernel linear algebra are known, but they all invariably take time exponential in either the dimension of the input point set (e.g., fast multipole methods suffer from the curse of dimensionality) or in the degree of the kernel function. Oblivious sketching has emerged as a powerful approach to speeding up numerical linear algebra over the past decade, but our understanding of oblivious sketching solutions for kernel matrices has remained quite limited, suffering from the aforementioned exponential dependence on input parameters. Our main contribution is a general method for applying sketching solutions developed in numerical linear algebra over the past decade to a tensoring of data points without forming the tensoring explicitly. This leads to the first oblivious sketch for the polynomial kernel with a target dimension that is only polynomially dependent on the degree of the kernel function, as well as the first oblivious sketch for the Gaussian kernel on bounded datasets that does not suffer from an exponential dependence on the dimensionality of input data points. △ Less

Submitted 22 December, 2020; v1 submitted 3 September, 2019; originally announced September 2019.

Fast hashing with Strong Concentration Bounds

Authors: Anders Aamand, Jakob B. T. Knudsen, Mathias B. T. Knudsen, Peter M. R. Rasmussen, Mikkel Thorup

Abstract: Previous work on tabulation hashing by Patrascu and Thorup from STOC'11 on simple tabulation and from SODA'13 on twisted tabulation offered Chernoff-style concentration bounds on hash based sums, e.g., the number of balls/keys hashing to a given bin, but under some quite severe restrictions on the expected values of these sums. The basic idea in tabulation hashing is to view a key as consisting of… ▽ More Previous work on tabulation hashing by Patrascu and Thorup from STOC'11 on simple tabulation and from SODA'13 on twisted tabulation offered Chernoff-style concentration bounds on hash based sums, e.g., the number of balls/keys hashing to a given bin, but under some quite severe restrictions on the expected values of these sums. The basic idea in tabulation hashing is to view a key as consisting of $c=O(1)$ characters, e.g., a 64-bit key as $c=8$ characters of 8-bits. The character domain $Σ$ should be small enough that character tables of size $|Σ|$ fit in fast cache. The schemes then use $O(1)$ tables of this size, so the space of tabulation hashing is $O(|Σ|)$. However, the concentration bounds by Patrascu and Thorup only apply if the expected sums are $\ll |Σ|$. To see the problem, consider the very simple case where we use tabulation hashing to throw $n$ balls into $m$ bins and want to analyse the number of balls in a given bin. With their concentration bounds, we are fine if $n=m$, for then the expected value is $1$. However, if $m=2$, as when tossing $n$ unbiased coins, the expected value $n/2$ is $\gg |Σ|$ for large data sets, e.g., data sets that do not fit in fast cache. To handle expectations that go beyond the limits of our small space, we need a much more advanced analysis of simple tabulation, plus a new tabulation technique that we call \emph{tabulation-permutation} hashing which is at most twice as slow as simple tabulation. No other hashing scheme of comparable speed offers similar Chernoff-style concentration bounds. △ Less

Submitted 10 August, 2020; v1 submitted 1 May, 2019; originally announced May 2019.

Comments: 54 pages, 3 figures. An extended abstract appeared at the 52nd Annual ACM Symposium on Theory of Computing (STOC20)

Subsets and Supermajorities: Optimal Hashing-based Set Similarity Search

Authors: Thomas Dybdahl Ahle, Jakob Bæk Tejs Knudsen

Abstract: We formulate and optimally solve a new generalized Set Similarity Search problem, which assumes the size of the database and query sets are known in advance. By creating polylog copies of our data-structure, we optimally solve any symmetric Approximate Set Similarity Search problem, including approximate versions of Subset Search, Maximum Inner Product Search (MIPS), Jaccard Similarity Search and… ▽ More We formulate and optimally solve a new generalized Set Similarity Search problem, which assumes the size of the database and query sets are known in advance. By creating polylog copies of our data-structure, we optimally solve any symmetric Approximate Set Similarity Search problem, including approximate versions of Subset Search, Maximum Inner Product Search (MIPS), Jaccard Similarity Search and Partial Match. Our algorithm can be seen as a natural generalization of previous work on Set as well as Euclidean Similarity Search, but conceptually it differs by optimally exploiting the information present in the sets as well as their complements, and doing so asymmetrically between queries and stored sets. Doing so we improve upon the best previous work: MinHash [J. Discrete Algorithms 1998], SimHash [STOC 2002], Spherical LSF [SODA 2016, 2017] and Chosen Path [STOC 2017] by as much as a factor $n^{0.14}$ in both time and space; or in the near-constant time regime, in space, by an arbitrarily large polynomial factor. Turning the geometric concept, based on Boolean supermajority functions, into a practical algorithm requires ideas from branching random walks on $\mathbb Z^2$, for which we give the first non-asymptotic near tight analysis. Our lower bounds follow from new hypercontractive arguments, which can be seen as characterizing the exact family of similarity search problems for which supermajorities are optimal. The optimality holds for among all hashing based data structures in the random setting, and by reductions, for 1 cell and 2 cell probe data structures. As a side effect, we obtain new hypercontractive bounds on the directed noise operator $T^{p_1 \to p_2}_ρ$. △ Less

Submitted 20 April, 2020; v1 submitted 8 April, 2019; originally announced April 2019.

MSC Class: 68W01

Classifying Convex Bodies by their Contact and Intersection Graphs

Authors: Anders Aamand, Mikkel Abrahamsen, Jakob Bæk Tejs Knudsen, Peter Michael Reichstein Rasmussen

Abstract: Suppose that $A$ is a convex body in the plane and that $A_1,\dots,A_n$ are translates of $A$. Such translates give rise to an intersection graph of $A$, $G=(V,E)$, with vertices $V=\{1,\dots,n\}$ and edges $E=\{uv\mid A_u\cap A_v\neq \emptyset\}$. The subgraph $G'=(V, E')$ satisfying that $E'\subset E$ is the set of edges $uv$ for which the interiors of $A_u$ and $A_v$ are disjoint is a unit dist… ▽ More Suppose that $A$ is a convex body in the plane and that $A_1,\dots,A_n$ are translates of $A$. Such translates give rise to an intersection graph of $A$, $G=(V,E)$, with vertices $V=\{1,\dots,n\}$ and edges $E=\{uv\mid A_u\cap A_v\neq \emptyset\}$. The subgraph $G'=(V, E')$ satisfying that $E'\subset E$ is the set of edges $uv$ for which the interiors of $A_u$ and $A_v$ are disjoint is a unit distance graph of $A$. If furthermore $G'=G$, i.e., if the interiors of $A_u$ and $A_v$ are disjoint whenever $u\neq v$, then $G$ is a contact graph of $A$. In this paper we study which pairs of convex bodies have the same contact, unit distance, or intersection graphs. We say that two convex bodies $A$ and $B$ are equivalent if there exists a linear transformation $B'$ of $B$ such that for any slope, the longest line segments with that slope contained in $A$ and $B'$, respectively, are equally long. For a broad class of convex bodies, including all strictly convex bodies and linear transformations of regular polygons, we show that the contact graphs of $A$ and $B$ are the same if and only if $A$ and $B$ are equivalent. We prove the same statement for unit distance and intersection graphs. △ Less

Submitted 5 February, 2019; originally announced February 2019.

Comments: 19 pages, 7 figures

MSC Class: 52C05

Absorption and luminescence spectroscopy of mass-selected Flavin Adenine Dinucleotide mono-anions

Authors: L. Giacomozzi, C. Kjær, J. Langeland Knudsen, L. H. Andersen, S. Brøndsted Nielsen, M. H. Stockett

Abstract: We report the absorption profile of isolated Flavin Adenine Dinucleotide (FAD) mono-anions recorded using Photo-Induced Dissociation action spectroscopy. In this charge state, one of the phosphoric acid groups is deprotonated and the chromophore itself is in its neutral oxidized state. These measurements cover the first four optical transitions of FAD with excitation energies from 2.3 to 6.0~eV (2… ▽ More We report the absorption profile of isolated Flavin Adenine Dinucleotide (FAD) mono-anions recorded using Photo-Induced Dissociation action spectroscopy. In this charge state, one of the phosphoric acid groups is deprotonated and the chromophore itself is in its neutral oxidized state. These measurements cover the first four optical transitions of FAD with excitation energies from 2.3 to 6.0~eV (210--550~nm). The $S_0\rightarrow S_2$ transition is strongly blue-shifted relative to aqueous solution, supporting the view that this transition has significant charge-transfer character. The remaining bands are close to their solution-phase positions. This confirms that the large discrepancy between quantum chemical calculations of vertical transition energies and solution-phase band maxima can not be explained by solvent effects. We also report the luminescence spectrum of FAD mono-anions \textit{in vacuo}. The gas-phase Stokes shift for $S_1$ is 3000~cm$^{-1}$, which is considerably larger than any previously reported for other molecular ions and consistent with a significant displacement of the ground and excited state potential energy surfaces. Consideration of vibronic structure is thus essential for simulating the absorption and luminescence spectra of flavins. △ Less

Submitted 14 May, 2018; originally announced May 2018.

Adsorption of CO on the Fe3O4(001) Surface

Authors: Jan Hulva, Zdeněk Jakub, Zbynek Novotny, Niclas Johansson, Jan Knudsen, Joachim Schnadt, Michael Schmid, Ulrike Diebold, Gareth S. Parkinson

Abstract: The interaction of CO with the Fe3O4(001)-(rt2xrt2)R45° surface was studied using temperature programmed desorption (TPD), scanning tunneling microscopy (STM) and x-ray photoelectron spectroscopy (XPS), the latter both under ultrahigh vacuum (UHV) conditions and in CO pressures up to 1 mbar. In general, the CO-Fe3O4 interaction is found to be weak. The strongest adsorption occurs at surface defect… ▽ More The interaction of CO with the Fe3O4(001)-(rt2xrt2)R45° surface was studied using temperature programmed desorption (TPD), scanning tunneling microscopy (STM) and x-ray photoelectron spectroscopy (XPS), the latter both under ultrahigh vacuum (UHV) conditions and in CO pressures up to 1 mbar. In general, the CO-Fe3O4 interaction is found to be weak. The strongest adsorption occurs at surface defects, leading to small TPD peaks at 115 K, 130 K and 190 K. Desorption from the regular surface occurs in two distinct regimes. For coverages up to 2 CO molecules per (rt2xrt2)R45° unit cell, the desorption maximum shows a large shift with increasing coverage, from initially 105 K to 70 K. For coverages between 2 and 4 molecules per (rt2xrt2)R45° unit cell, a much sharper desorption feature emerges at 50 K. Thermodynamic analysis of the TPD data suggests a phase transition from a dilute 2D gas into an ordered overlayer with CO molecules bound to surface Fe3+ sites. XPS data acquired at 45 K in UHV are consistent with physisorption. Some carbon-containing species are observed in the near-ambient-pressure XPS experiments at room temperature, but are attributed to contamination and/or reaction with CO with water from the residual gas. No evidence was found for surface reduction or carburization by CO molecules. △ Less

Submitted 13 April, 2018; originally announced April 2018.

Journal ref: J. Phys. Chem. B 122 (2018) 721-729

Low-altitude ion heating, downflowing ions, and BBELF waves in the return current region

Authors: Yangyang Shen, David J. Knudsen, Johnathan K. Burchill, Andrew D. Howarth, Andrew W. Yau, David M. Miles, H. Gordon James, Gareth W. Perry, Leroy Cogger

Abstract: Heavy (O+) ion energization and field-aligned motion in and near the ionosphere are still not well understood. Based on observations from the CASSIOPE Enhanced Polar Outflow Probe (e-POP) at altitudes between 325 km and 730 km over one year, we present a statistical study (24 events) of ion heating and its relation to field-aligned ion bulk flow velocity, low-frequency waves and field-aligned curr… ▽ More Heavy (O+) ion energization and field-aligned motion in and near the ionosphere are still not well understood. Based on observations from the CASSIOPE Enhanced Polar Outflow Probe (e-POP) at altitudes between 325 km and 730 km over one year, we present a statistical study (24 events) of ion heating and its relation to field-aligned ion bulk flow velocity, low-frequency waves and field-aligned currents (FACs). The ion temperature and field-aligned bulk flow velocity are derived from 2-D ion velocity distribution functions measured by the suprathermal electron imager (SEI) instrument. Consistent ion heating and flow velocity characteristics are observed from both the SEI and the rapid-scanning ion mass spectrometer (IRM) instruments. We find that transverse O+ ion heating in the ionosphere can be intense (up to 4.5 eV), confined to very narrow regions (~ 2 km across B), is more likely to occur in the downward current region, and is associated with broadband extremely low frequency (BBELF) waves. These waves are interpreted as linearly polarized perpendicular to the magnetic field. The amount of ion heating cannot be explained by frictional heating, and the correlation of ion heating with BBELF waves suggest that significant wave-ion heating is occurring and even dominating at altitudes as low as 350 km, a boundary that is lower than previously reported. Surprisingly, the majority of these heating events (17 out 24) are associated with core ion downflows rather than upflows. This may be explained by a downward-pointing electric field in the low-altitude return current region. △ Less

Submitted 10 March, 2018; originally announced March 2018.

Comments: 17 pages, including 12 figures, 1 table, accepted by JGR-Space Physics

Efficient Modal Decomposition of Vortex Beams via Holographically Reconstructed Phase

Authors: Jasmine M. Knudsen, Samuel N. Alperin, Andrew A. Voitiv, William G. Holtzmann, Juliet T. Gopinath, Mark E. Siemens

Abstract: We use phase-shifting digital holography to measure the amplitude and phase of twisted light. In our experiment, a spatial light modulator generates the studied vortex beams in addition to a co-propagating reference beam with a controllable relative phase. We show complex field measurements for single and superposition Laguerre Gaussian (LG) modes, demonstrate full modal decompositions into LG and… ▽ More We use phase-shifting digital holography to measure the amplitude and phase of twisted light. In our experiment, a spatial light modulator generates the studied vortex beams in addition to a co-propagating reference beam with a controllable relative phase. We show complex field measurements for single and superposition Laguerre Gaussian (LG) modes, demonstrate full modal decompositions into LG and orbital angular momentum (OAM) power spectral bases, provide error analysis and demonstrate high insensitivity to detector misalignment to show the robustness of the technique. This enables rapid determination of OAM spectra with low uncertainty, allowing us to report the first vortex beams with 99.9% purity. △ Less

Submitted 8 October, 2017; originally announced October 2017.

Swarm in situ observations of F region polar cap patches created by cusp precipitation

Authors: L. V. Goodwin, B. Iserhienrhien, D. M. Miles, S. Patra, C. van der Meeren, S. C. Buchert, J. K. Burchill, L. B. N. Clausen, D. J. Knudsen, K. A. McWilliams, J. Moen

Abstract: High-resolution in situ measurements from the three Swarm spacecraft, in a string-of-pearls configuration, provide new insights about the combined role of flow channel events and particle impact ionization in creating $\textit{F}$ region electron density structures in the northern Scandinavian dayside cusp. We present a case of polar cap patch formation where a reconnection-driven low-density rela… ▽ More High-resolution in situ measurements from the three Swarm spacecraft, in a string-of-pearls configuration, provide new insights about the combined role of flow channel events and particle impact ionization in creating $\textit{F}$ region electron density structures in the northern Scandinavian dayside cusp. We present a case of polar cap patch formation where a reconnection-driven low-density relative westward flow channel is eroding the dayside solar-ionized plasma but where particle impact ionization in the cusp dominates the initial plasma structuring. In the cusp, density features are observed which are twice as dense as the solar-ionized background. These features then follow the polar cap convection and become less structured and lower in amplitude. These are the first in situ observations tracking polar cap patch evolution from creation by plasma transport and enhancement by cusp precipitation, through entrainment in the polar cap flow and relaxation into smooth patches as they approach the nightside auroral oval. △ Less

Submitted 7 June, 2016; originally announced June 2016.

Xe Irradiation of Graphene on Ir(111): From Trapping to Blistering

Authors: Charlotte Herbig, E. Harriet Åhlgren, Ulrike A. Schröder, Antonio J. Martínez-Galera, Mohammad A. Arman, Jani Kotakoski, Jan Knudsen, Arkady V. Krasheninnikov, Thomas Michely

Abstract: Using X-ray photoelectron spectroscopy, thermal desorption spectroscopy, and scanning tunneling microscopy we show that upon keV Xe + irradiation of graphene on Ir(111), Xe atoms are trapped under the graphene. Upon annealing, aggregation of Xe leads to graphene bulges and blisters. The efficient trapping is an unexpected and remarkable phenomenon, given the absence of chemical binding of Xe to Ir… ▽ More Using X-ray photoelectron spectroscopy, thermal desorption spectroscopy, and scanning tunneling microscopy we show that upon keV Xe + irradiation of graphene on Ir(111), Xe atoms are trapped under the graphene. Upon annealing, aggregation of Xe leads to graphene bulges and blisters. The efficient trapping is an unexpected and remarkable phenomenon, given the absence of chemical binding of Xe to Ir and to graphene, the weak interaction of a perfect graphene layer with Ir(111), as well as the substantial damage to graphene due to irradiation. By combining molecular dynamics simulations and density functional theory calculations with our experiments, we uncover the mechanism of trapping. We describe ways to avoid blister formation during graphene growth, and also demonstrate how ion implantation can be used to intentionally create blisters without introducing damage to the graphene layer. Our approach may provide a pathway to synthesize new materials at a substrate - 2D material interface or to enable confined reactions at high pressures and temperatures. △ Less

Submitted 12 August, 2015; originally announced August 2015.

Journal ref: Phys. Rev. B 92, 085429 (2015)

On Graphs and Codes Preserved by Edge Local Complementation

Authors: Lars Eirik Danielsen, Matthew G. Parker, Constanza Riera, Joakim Grahl Knudsen

Abstract: Orbits of graphs under local complementation (LC) and edge local complementation (ELC) have been studied in several different contexts. For instance, there are connections between orbits of graphs and error-correcting codes. We define a new graph class, ELC-preserved graphs, comprising all graphs that have an ELC orbit of size one. Through an exhaustive search, we find all ELC-preserved graphs of… ▽ More Orbits of graphs under local complementation (LC) and edge local complementation (ELC) have been studied in several different contexts. For instance, there are connections between orbits of graphs and error-correcting codes. We define a new graph class, ELC-preserved graphs, comprising all graphs that have an ELC orbit of size one. Through an exhaustive search, we find all ELC-preserved graphs of order up to 12 and all ELC-preserved bipartite graphs of order up to 16. We provide general recursive constructions for infinite families of ELC-preserved graphs, and show that all known ELC-preserved graphs arise from these constructions or can be obtained from Hamming codes. We also prove that certain pairs of ELC-preserved graphs are LC equivalent. We define ELC-preserved codes as binary linear codes corresponding to bipartite ELC-preserved graphs, and study the parameters of such codes. △ Less

Submitted 8 August, 2013; v1 submitted 30 June, 2010; originally announced June 2010.

Comments: 23 pages, 16 figures

MSC Class: 94B05 (Primary) 05C30; 94A15; 05C75 (Secondary)

Iterative Decoding on Multiple Tanner Graphs Using Random Edge Local Complementation

Authors: Joakim Grahl Knudsen, Constanza Riera, Lars Eirik Danielsen, Matthew G. Parker, Eirik Rosnes

Abstract: In this paper, we propose to enhance the performance of the sum-product algorithm (SPA) by interleaving SPA iterations with a random local graph update rule. This rule is known as edge local complementation (ELC), and has the effect of modifying the Tanner graph while preserving the code. We have previously shown how the ELC operation can be used to implement an iterative permutation group decod… ▽ More In this paper, we propose to enhance the performance of the sum-product algorithm (SPA) by interleaving SPA iterations with a random local graph update rule. This rule is known as edge local complementation (ELC), and has the effect of modifying the Tanner graph while preserving the code. We have previously shown how the ELC operation can be used to implement an iterative permutation group decoder (SPA-PD)--one of the most successful iterative soft-decision decoding strategies at small blocklengths. In this work, we exploit the fact that ELC can also give structurally distinct parity-check matrices for the same code. Our aim is to describe a simple iterative decoder, running SPA-PD on distinct structures, based entirely on random usage of the ELC operation. This is called SPA-ELC, and we focus on small blocklength codes with strong algebraic structure. In particular, we look at the extended Golay code and two extended quadratic residue codes. Both error rate performance and average decoding complexity, measured by the average total number of messages required in the decoding, significantly outperform those of the standard SPA, and compares well with SPA-PD. However, in contrast to SPA-PD, which requires a global action on the Tanner graph, we obtain a performance improvement via local action alone. Such localized algorithms are of mathematical interest in their own right, but are also suited to parallel/distributed realizations. △ Less

Submitted 26 May, 2009; originally announced May 2009.

Comments: 5 pages, to appear in proc. IEEE ISIT, June 2009

Experimental and theoretical study of oxygen adsorption structures on Ag(111)

Authors: Joachim Schnadt, Jan Knudsen, Xiao Liang Hu, Angelos Michaelides, Ronnie T. Vang, Karsten Reuter, Zheshen Li, Erik Lægsgaard, Matthias Scheffler, Flemming Besenbacher

Abstract: The oxidized Ag(111) surface has been studied by a combination of experimental and theoretical methods, scanning tunneling microscopy (STM), x-ray photoelectron spectroscopy (XPS), and density functional theory (DFT). A large variety of different surface structures is found, depending on the detailed preparation conditions. The observed structures fall into four classes: (a) individually chemiso… ▽ More The oxidized Ag(111) surface has been studied by a combination of experimental and theoretical methods, scanning tunneling microscopy (STM), x-ray photoelectron spectroscopy (XPS), and density functional theory (DFT). A large variety of different surface structures is found, depending on the detailed preparation conditions. The observed structures fall into four classes: (a) individually chemisorbed atomic oxygen atoms, (b) three different oxygen overlayer structures, including the well-known p(4 x 4) phase, formed from the same Ag$_6$ and Ag$_{10}$ building blocks, (c) a c(4 x 8) structure not previously observed, and (d) at higher oxygen coverages structures characterized by stripes along the high-symmetry directions of the Ag(111) substrate. Our analysis provides a detailed explanation of the atomic-scale geometry of the Ag$_6$/Ag$_{10}$ building block structures, and the c(4 x 8) and stripe structures are discussed in detail. The observation of many different and co-existing structures implies that the O/Ag(111) system is characterized by a significantly larger degree of complexity than previously anticipated, and this will impact our understanding of oxidation catalysis processes on Ag catalysts. △ Less

Submitted 20 August, 2009; v1 submitted 23 April, 2009; originally announced April 2009.

Comments: 12 pages, 9 figures; peer-reviewed and published version

Journal ref: Phys. Rev. B 80 (2009) 075424

Breakdown of disordered media by surface loads

Authors: Jakob Knudsen, A. R. Massih

Abstract: We model an interface layer connecting two parts of a solid body by N parallel elastic springs connecting two rigid blocks. We load the system by a shear force acting on the top side. The springs have equal stiffness but are ruptured randomly when the load reaches a critical value. For the considered system, we calculate the shear modulus, G, as a function of the order parameter, φ, describing t… ▽ More We model an interface layer connecting two parts of a solid body by N parallel elastic springs connecting two rigid blocks. We load the system by a shear force acting on the top side. The springs have equal stiffness but are ruptured randomly when the load reaches a critical value. For the considered system, we calculate the shear modulus, G, as a function of the order parameter, φ, describing the state of damage, and also the ``spalled'' material (burst) size distribution. In particular, we evaluate the relation between the damage parameter and the applied force and explore the behaviour in the vicinity of material breakdown. Using this simple model for material breakdown, we show that damage, caused by applied shear forces, is analogous to a first-order phase transition. The scaling behaviour of G with φis explored analytically and numerically, close to φ=0 and φ=1 and in the vicinity of φ_c, when the shear load is close but below the threshold force that causes material breakdown. Our model calculation represents a first approximation of a system subject to wear induced loads. △ Less

Submitted 12 October, 2004; originally announced October 2004.

Comments: 15 pages, 7 figures