Tectonic evolution of a continental collision zone: A thermomechanical numerical model

We model evolution of a continent-continent collision and draw some parallels with the tectonic evolution of the Himalaya. We use a large-scale visco-plasto-elastic thermomechanical model that has a free upper surface, accounts for erosion and deposition and allows for all modes of lithospheric deformation. For quartz/olivine rheology and 60 mm/yr convergence rate, the continental subduction is stable, and the model predicts three distinct phases. During the phase 1 (120 km or 6% of shortening), deformation is characterized by back thrusting around the suture zone. Some amount of delaminated lower crust accumulates at depth. During phase 2 (120 km–420 km or 6%–22% of shortening), this crustal root is exhumed (medium- to high-grade rocks) along a newly formed major thrust fault. This stage bears similarities with the period of coeval activity of the Main Central thrust and of the South Tibetan Detachment between 20–16 Myr ago. During phase 3 (>420 km or 22% of shortening), the crust is scraped off from the mantle lithosphere and is incorporated into large crustal wedge. Deformation is localized around frontal thrust faults. This kinematics should produce only low- to medium-grade exhumation. This stage might be compared with the tectonics that has prevailed in the Himalaya over the last 15 Myr allowing for the formation of the Lesser Himalaya. The experiment is conducted at constant convergence rate, which implies increasing compressive force. Considering that this force is constant in nature, this result may be equivalent to a slowing down of the convergence rate as was observed during the India-Asia collision

The strange quark condensate in the nucleon in 2+1 flavor QCD

We calculate the "strange quark content of the nucleon", , which is important for interpreting the results of some dark matter detection experiments. The method is to evaluate quark-line disconnected correlations on the MILC lattice ensembles, which include the effects of dynamical strange quarks. After continuum and chiral extrapolations, the result is <N |s s_bar |N> = 0.69 +- 0.07(statistical) +- 0.09(systematic), in the modified minimal subtraction scheme (2 GeV), or for the renormalization scheme invariant form, m_s partial{M_N}/partial{m_s} = 59(6)(8) MeV.Comment: Added figures and references, especially for fit range choice. Other changes for clarity. Version to appear in publicatio

Granular Rayleigh-Taylor Instability: Experiments and Simulations

A granular instability driven by gravity is studied experimentally and numerically. The instability arises as grains fall in a closed Hele-Shaw cell where a layer of dense granular material is positioned above a layer of air. The initially flat front defined by the grains subsequently develops into a pattern of falling granular fingers separated by rising bubbles of air. A transient coarsening of the front is observed right from the start by a finger merging process. The coarsening is later stabilized by new fingers growing from the center of the rising bubbles. The structures are quantified by means of Fourier analysis and quantitative agreement between experiment and computation is shown. This analysis also reveals scale invariance of the flow structures under overall change of spatial scale.Comment: 4 pages, 11 figure

Interplay of seismic and aseismic deformations during earthquake swarms: An experimental approach

Observations of earthquake swarms and slow propagating ruptures on related faults suggest a close relation between the two phenomena. Earthquakes are the signature of fast unstable ruptures initiated on localized asperities while slow aseismic deformations are experienced on large stable segments of the fault plane. The spatial proximity and the temporal coincidence of both fault mechanical responses highlight the variability of fault rheology. However, the mechanism relating earthquakes and aseismic processes is still elusive due to the difficulty of imaging these phenomena of large spatiotemporal variability at depth. Here we present laboratory experiments that explore, in great detail, the deformation processes of heterogeneous interfaces in the brittle-creep regime. We track the evolution of an interfacial crack over 7 orders of magnitude in time and 5 orders of magnitude in space using optical and acoustic sensors. We explore the response of the system to slow transient loads and show that slow deformation episodes are systematically accompanied by acoustic emissions due to local fracture energy disorder. Features of acoustic emission activities and deformation rate distributions of our experimental system are similar to those in natural faults. On the basis of an activation energy model, we link our results to the Rate and State friction model and suggest an active role of local creep deformation in driving the seismic activity of earthquake swarms

Influence of pore-scale disorder on viscous fingering during drainage

We study viscous fingering during drainage experiments in linear Hele-Shaw cells filled with a random porous medium. The central zone of the cell is found to be statistically more occupied than the average, and to have a lateral width of 40% of the system width, irrespectively of the capillary number $Ca$. A crossover length $w_f \propto Ca^{-1}$ separates lower scales where the invader's fractal dimension $D\simeq1.83$ is identical to capillary fingering, and larger scales where the dimension is found to be $D\simeq1.53$. The lateral width and the large scale dimension are lower than the results for Diffusion Limited Aggregation, but can be explained in terms of Dielectric Breakdown Model. Indeed, we show that when averaging over the quenched disorder in capillary thresholds, an effective law $v\propto (\nabla P)^2$ relates the average interface growth rate and the local pressure gradient.Comment: 4 pages, 4 figures, submitted to Phys Rev Letter

Integrated music and math projects in secondary education

The introduction of projects involving music and mathematics in Secondary Education should allow the integration of these disciplines by nonspecialists. In the present work we describe an experience carried out with future mathematics teachers with solid scientific-technical training, but little musical training. This pilot contributes with concrete orientations and results para the creation and development of STEAM activities, which can be found in [5].This research work has been carried out within the framework of the project EDU2017-84979-R, of the Spanish State Program of R&D and Innovation Oriented to the Challenges of the Society

Persistence in One-dimensional Ising Models with Parallel Dynamics

We study persistence in one-dimensional ferromagnetic and anti-ferromagnetic nearest-neighbor Ising models with parallel dynamics. The probability P(t) that a given spin has not flipped up to time t, when the system evolves from an initial random configuration, decays as P(t) \sim 1/t^theta_p with theta_p \simeq 0.75 numerically. A mapping to the dynamics of two decoupled A+A \to 0 models yields theta_p = 3/4 exactly. A finite size scaling analysis clarifies the nature of dynamical scaling in the distribution of persistent sites obtained under this dynamics.Comment: 5 pages Latex file, 3 postscript figures, to appear in Phys Rev.

A Novel Approach for Ellipsoidal Outer-Approximation of the Intersection Region of Ellipses in the Plane

In this paper, a novel technique for tight outer-approximation of the intersection region of a finite number of ellipses in 2-dimensional (2D) space is proposed. First, the vertices of a tight polygon that contains the convex intersection of the ellipses are found in an efficient manner. To do so, the intersection points of the ellipses that fall on the boundary of the intersection region are determined, and a set of points is generated on the elliptic arcs connecting every two neighbouring intersection points. By finding the tangent lines to the ellipses at the extended set of points, a set of half-planes is obtained, whose intersection forms a polygon. To find the polygon more efficiently, the points are given an order and the intersection of the half-planes corresponding to every two neighbouring points is calculated. If the polygon is convex and bounded, these calculated points together with the initially obtained intersection points will form its vertices. If the polygon is non-convex or unbounded, we can detect this situation and then generate additional discrete points only on the elliptical arc segment causing the issue, and restart the algorithm to obtain a bounded and convex polygon. Finally, the smallest area ellipse that contains the vertices of the polygon is obtained by solving a convex optimization problem. Through numerical experiments, it is illustrated that the proposed technique returns a tighter outer-approximation of the intersection of multiple ellipses, compared to conventional techniques, with only slightly higher computational cost