Inexact Newton Dogleg Methods

The dogleg method is a classical trust-region technique for globalizing Newton\u27s method. While it is widely used in optimization, including large-scale optimization via truncated-Newton approaches, its implementation in general inexact Newton methods for systems of nonlinear equations can be problematic. In this paper, we first outline a very general dogleg method suitable for the general inexact Newton context and provide a global convergence analysis for it. We then discuss certain issues that may arise with the standard dogleg implementational strategy and propose modified strategies that address them. Newton-Krylov methods have provided important motivation for this work, and we conclude with a report on numerical experiments involving a Newton-GMRES dogleg method applied to benchmark CFD problems

Construction and Performance of a Micro-Pattern Stereo Detector with Two Gas Electron Multipliers

The construction of a micro-pattern gas detector of dimensions 40x10 cm**2 is described. Two gas electron multiplier foils (GEM) provide the internal amplification stages. A two-layer readout structure was used, manufactured in the same technology as the GEM foils. The strips of each layer cross at an effective crossing angle of 6.7 degrees and have a 406 um pitch. The performance of the detector has been evaluated in a muon beam at CERN using a silicon telescope as reference system. The position resolutions of two orthogonal coordinates are measured to be 50 um and 1 mm, respectively. The muon detection efficiency for two-dimensional space points reaches 96%.Comment: 21 pages, 17 figure

Cones, pringles, and grain boundary landscapes in graphene topology

A polycrystalline graphene consists of perfect domains tilted at angle {\alpha} to each other and separated by the grain boundaries (GB). These nearly one-dimensional regions consist in turn of elementary topological defects, 5-pentagons and 7-heptagons, often paired up into 5-7 dislocations. Energy G({\alpha}) of GB computed for all range 0<={\alpha}<=Pi/3, shows a slightly asymmetric behavior, reaching ~5 eV/nm in the middle, where the 5's and 7's qualitatively reorganize in transition from nearly armchair to zigzag interfaces. Analysis shows that 2-dimensional nature permits the off-plane relaxation, unavailable in 3-dimensional materials, qualitatively reducing the energy of defects on one hand while forming stable 3D-landsapes on the other. Interestingly, while the GB display small off-plane elevation, the random distributions of 5's and 7's create roughness which scales inversely with defect concentration, h ~ n^(-1/2)Comment: 9 pages, 4 figure

Theoretical analysis of the role of chromatin interactions in long-range action of enhancers and insulators

Long-distance regulatory interactions between enhancers and their target genes are commonplace in higher eukaryotes. Interposed boundaries or insulators are able to block these long distance regulatory interactions. The mechanistic basis for insulator activity and how it relates to enhancer action-at-a-distance remains unclear. Here we explore the idea that topological loops could simultaneously account for regulatory interactions of distal enhancers and the insulating activity of boundary elements. We show that while loop formation is not in itself sufficient to explain action at a distance, incorporating transient non-specific and moderate attractive interactions between the chromatin fibers strongly enhances long-distance regulatory interactions and is sufficient to generate a euchromatin-like state. Under these same conditions, the subdivision of the loop into two topologically independent loops by insulators inhibits inter-domain interactions. The underlying cause of this effect is a suppression of crossings in the contact map at intermediate distances. Thus our model simultaneously accounts for regulatory interactions at a distance and the insulator activity of boundary elements. This unified model of the regulatory roles of chromatin loops makes several testable predictions that could be confronted with \emph{in vitro} experiments, as well as genomic chromatin conformation capture and fluorescent microscopic approaches.Comment: 10 pages, originally submitted to an (undisclosed) journal in May 201

Optimal General Matchings

Given a graph $G=(V,E)$ and for each vertex $v \in V$ a subset $B(v)$ of the set $\{0,1,\ldots, d_G(v)\}$, where $d_G(v)$ denotes the degree of vertex $v$ in the graph $G$, a $B$-factor of $G$ is any set $F \subseteq E$ such that $d_F(v) \in B(v)$ for each vertex $v$, where $d_F(v)$ denotes the number of edges of $F$ incident to $v$. The general factor problem asks the existence of a $B$-factor in a given graph. A set $B(v)$ is said to have a {\em gap of length} $p$ if there exists a natural number $k \in B(v)$ such that $k+1, \ldots, k+p \notin B(v)$ and $k+p+1 \in B(v)$. Without any restrictions the general factor problem is NP-complete. However, if no set $B(v)$ contains a gap of length greater than $1$, then the problem can be solved in polynomial time and Cornuejols \cite{Cor} presented an algorithm for finding a $B$-factor, if it exists. In this paper we consider a weighted version of the general factor problem, in which each edge has a nonnegative weight and we are interested in finding a $B$-factor of maximum (or minimum) weight. In particular, this version comprises the minimum/maximum cardinality variant of the general factor problem, where we want to find a $B$-factor having a minimum/maximum number of edges. We present an algorithm for the maximum/minimum weight $B$-factor for the case when no set $B(v)$ contains a gap of length greater than $1$. This also yields the first polynomial time algorithm for the maximum/minimum cardinality $B$-factor for this case

Electronic transport in polycrystalline graphene

Most materials in available macroscopic quantities are polycrystalline. Graphene, a recently discovered two-dimensional form of carbon with strong potential for replacing silicon in future electronics, is no exception. There is growing evidence of the polycrystalline nature of graphene samples obtained using various techniques. Grain boundaries, intrinsic topological defects of polycrystalline materials, are expected to dramatically alter the electronic transport in graphene. Here, we develop a theory of charge carrier transmission through grain boundaries composed of a periodic array of dislocations in graphene based on the momentum conservation principle. Depending on the grain boundary structure we find two distinct transport behaviours - either high transparency, or perfect reflection of charge carriers over remarkably large energy ranges. First-principles quantum transport calculations are used to verify and further investigate this striking behaviour. Our study sheds light on the transport properties of large-area graphene samples. Furthermore, purposeful engineering of periodic grain boundaries with tunable transport gaps would allow for controlling charge currents without the need of introducing bulk band gaps in otherwise semimetallic graphene. The proposed approach can be regarded as a means towards building practical graphene electronics.Comment: accepted in Nature Material

Effect of the disorder in graphene grain boundaries: A wave packet dynamics study

Chemical vapor deposition (CVD) on Cu foil is one of the most promising methods to produce graphene samples despite of introducing numerous grain boundaries into the perfect graphene lattice. A rich variety of GB structures can be realized experimentally by controlling the parameters in the CVD method. Grain boundaries contain non-hexagonal carbon rings (4, 5, 7, 8 membered rings) and vacancies in various ratios and arrangements. Using wave packet dynamic (WPD) simulations and tight-binding electronic structure calculations, we have studied the effect of the structure of GBs on the transport properties. Three model GBs with increasing disorder were created in the computer: a periodic 5-7 GB, a "serpentine" GB, and a disordered GB containing 4, 8 membered rings and vacancies. It was found that for small energies (E = EF ± 1 eV) the transmission decreases with increasing disorder. Four membered rings and vacancies are identified as the principal scattering centers. Revealing the connection between the properties of GBs and the CVD growth method may open new opportunities in the graphene based nanoelectronics. © 2013 Elsevier B.V. All rights reserved

Enrichment analysis of Alu elements with different spatial chromatin proximity in the human genome

Transposable elements (TEs) have no longer been totally considered as “junk DNA” for quite a time since the continual discoveries of their multifunctional roles in eukaryote genomes. As one of the most important and abundant TEs that still active in human genome, Alu, a SINE family, has demonstrated its indispensable regulatory functions at sequence level, but its spatial roles are still unclear. Technologies based on 3C(chromosomeconformation capture) have revealed the mysterious three-dimensional structure of chromatin, and make it possible to study the distal chromatin interaction in the genome. To find the role TE playing in distal regulation in human genome, we compiled the new released Hi-C data, TE annotation, histone marker annotations, and the genome-wide methylation data to operate correlation analysis, and found that the density of Alu elements showed a strong positive correlation with the level of chromatin interactions (hESC: r=0.9, P<2.2×1016; IMR90 fibroblasts: r = 0.94, P < 2.2 × 1016) and also have a significant positive correlation withsomeremote functional DNA elements like enhancers and promoters (Enhancer: hESC: r=0.997, P=2.3×10−4; IMR90: r=0.934, P=2×10−2; Promoter: hESC: r = 0.995, P = 3.8 × 10−4; IMR90: r = 0.996, P = 3.2 × 10−4). Further investigation involving GC content and methylation status showed the GC content of Alu covered sequences shared a similar pattern with that of the overall sequence, suggesting that Alu elements also function as the GC nucleotide and CpG site provider. In all, our results suggest that the Alu elements may act as an alternative parameter to evaluate the Hi-C data, which is confirmed by the correlation analysis of Alu elements and histone markers. Moreover, the GC-rich Alu sequence can bring high GC content and methylation flexibility to the regions with more distal chromatin contact, regulating the transcription of tissue-specific genes