• ### Assessing Visual Quality of Omnidirectional Videos(1709.06342)

July 14, 2019 eess.IV
In contrast with traditional video, omnidirectional video enables spherical viewing direction with support for head-mounted displays, providing an interactive and immersive experience. Unfortunately, to the best of our knowledge, there are few visual quality assessment (VQA) methods, either subjective or objective, for omnidirectional video coding. This paper proposes both subjective and objective methods for assessing quality loss in encoding omnidirectional video. Specifically, we first present a new database, which includes the viewing direction data from several subjects watching omnidirectional video sequences. Then, from our database, we find a high consistency in viewing directions across different subjects. The viewing directions are normally distributed in the center of the front regions, but they sometimes fall into other regions, related to video content. Given this finding, we present a subjective VQA method for measuring difference mean opinion score (DMOS) of the whole and regional omnidirectional video, in terms of overall DMOS (O-DMOS) and vectorized DMOS (V-DMOS), respectively. Moreover, we propose two objective VQA methods for encoded omnidirectional video, in light of human perception characteristics of omnidirectional video. One method weighs the distortion of pixels with regard to their distances to the center of front regions, which considers human preference in a panorama. The other method predicts viewing directions according to video content, and then the predicted viewing directions are leveraged to allocate weights to the distortion of each pixel in our objective VQA method. Finally, our experimental results verify that both the subjective and objective methods proposed in this paper advance state-of-the-art VQA for omnidirectional video.
• ### Convolutional Sequence to Sequence Model for Human Dynamics(1805.00655)

May 2, 2018 cs.CV
Human motion modeling is a classic problem in computer vision and graphics. Challenges in modeling human motion include high dimensional prediction as well as extremely complicated dynamics.We present a novel approach to human motion modeling based on convolutional neural networks (CNN). The hierarchical structure of CNN makes it capable of capturing both spatial and temporal correlations effectively. In our proposed approach,a convolutional long-term encoder is used to encode the whole given motion sequence into a long-term hidden variable, which is used with a decoder to predict the remainder of the sequence. The decoder itself also has an encoder-decoder structure, in which the short-term encoder encodes a shorter sequence to a short-term hidden variable, and the spatial decoder maps the long and short-term hidden variable to motion predictions. By using such a model, we are able to capture both invariant and dynamic information of human motion, which results in more accurate predictions. Experiments show that our algorithm outperforms the state-of-the-art methods on the Human3.6M and CMU Motion Capture datasets. Our code is available at the project website.
• ### Sparse solutions in optimal control of PDEs with uncertain parameters: the linear case(1804.05678)

April 16, 2018 math.NA, stat.CO, math.OC
We study sparse solutions of optimal control problems governed by PDEs with uncertain coefficients. We propose two formulations, one where the solution is a deterministic control optimizing the mean objective, and a formulation aiming at stochastic controls that share the same sparsity structure. In both formulations, regions where the controls do not vanish can be interpreted as optimal locations for placing control devices. In this paper, we focus on linear PDEs with linearly entering uncertain parameters. Under these assumptions, the deterministic formulation reduces to a problem with known structure, and thus we mainly focus on the stochastic control formulation. Here, shared sparsity is achieved by incorporating the $L^1$-norm of the mean of the pointwise squared controls in the objective. We reformulate the problem using a norm reweighting function that is defined over physical space only and thus helps to avoid approximation of the random space. We show that a fixed point algorithm applied to the norm reweighting formulation leads to the well-studied iterative reweighted least squares (IRLS) algorithm, and we propose a novel preconditioned Newton-conjugate gradient method to speed up the IRLS algorithm. We combine our algorithms with low-rank operator approximations, for which we provide estimates of the truncation error. We carefully examine the computational complexity of the resulting algorithms. The sparsity structure of the optimal controls and the performance of the solution algorithms are studied numerically using control problems governed by the Laplace and Helmholtz equations with uncertain Neumann boundary data. In these experiments the Newton variant clearly outperforms the IRLS method.
• ### Dephasing and Relaxation of Bosons in 1D: Newton's Cradle Revisited(1804.01969)

We study the relaxation of ultra-cold bosons in one dimension. In a setting similar to the Newton's cradle experiment by Kinoshita et al.[Nature (London) 440, 900 (2006)], we excite the atoms to oscillate and collide in a 2D optical lattice of 1D tubes and observe the oscillations for up to 4.8 seconds (400 oscillations) with very little heating and loss. Tuning the initial distribution, we can probe the onset of relaxation in the cross over between 1D and 3D. Investigating the population of the transverse excited state by band mapping, we show that even for samples where the energy put into the system is initially not sufficient to excite transverse state, relaxation starts through a minimal (1%) seed of atoms in transverse excited states.
• ### Learning the Kernel for Classification and Regression(1712.08597)

Dec. 25, 2017 cs.LG
We investigate a series of learning kernel problems with polynomial combinations of base kernels, which will help us solve regression and classification problems. We also perform some numerical experiments of polynomial kernels with regression and classification tasks on different datasets.
• ### Transport properties of ultrathin BaFe1.84Co0.16As2 superconducting nanowires(1711.07637)

Nov. 21, 2017 cond-mat.supr-con
Superconducting nanowire single-photon detectors (SNSPDs) have an absolute advantage over other types of single photon detectors except the low operating temperature. Therefore, many efforts have been devoted to find high-temperature superconducting materials that are suitable for preparing SNSPDs. Copper-based and MgB2 ultra-thin superconducting nanowires have been already reported. However, the transport properties of iron-based ultra-thin superconducting nanowires have not been studied. In this work, a 10 nm thick 200 nm wide 30 {\mu}m long high quality superconducting nanowire was fabricated from ultrathin BaFe1.84Co0.16As2 films by a lift-off process. The precursor BaFe1.84Co0.16As2 film with a thickness of 10 nm and root-mean-square roughness of 1 nm was grown on CaF2 substrates by pulsed laser deposition. The nanowire shows a high superconducting critical temperature Tzero c=20 K with a narrow transition width of Delta T=2.5 K and exhibits a high critical current density Jc of 1.8E7 A/cm2 at 10 K. These results of ultrathin BaFe1.84Co0.16As2 nanowire will attract interest in electronic applications, including SNSPDs
• ### Density functional theory of electron transfer beyond the Born-Oppenheimer approximation: Case study of LiF(1710.11377)

Nov. 1, 2017 physics.chem-ph
We perform model calculations for a stretched LiF molecule, demonstrating that nonadiabatic charge transfer effects can be accurately and seamlessly described within a density functional framework. In alkali halides like LiF, there is an abrupt change in the ground state electronic distribution due to an electron transfer at a critical bond length $R=R_c$, where a barely avoided crossing of the lowest adiabatic potential energy surfaces calls the validity of the Born-Oppenheimer approximation into doubt. Modeling the $R$-dependent electronic structure of LiF within a two-site Hubbard model, we find that nonadiabatic electron-nuclear coupling produces a sizable elongation of the critical $R_c$ by 0.5 Bohr. This effect is very accurately captured by a simple and rigorously-derived correction, with an $M^{-1}$ prefactor, to the exchange-correlation potential in density functional theory; $M=$ reduced nuclear mass. Since this nonadiabatic term depends on gradients of the nuclear wavefunction and conditional electronic density, $\nabla_R \chi(R)$ and $\nabla_R n(\mathbf{r},R)$, it couples the Kohn-Sham equations at neighboring $R$ points. Motivated by an observed localization of nonadiabatic effects in nuclear configuration space, we propose a local conditional density approximation -- an approximation that reduces the search for nonadiabatic density functionals to the search for a single function $y(n)$.
• ### Chaos-based Wireless Communication Resisting Multipath Effects(1611.02325)

Sept. 26, 2017 cs.IT, math.IT
In additive white gaussian noise (AWGN) channel, chaos has been proved to be the optimal coherent communication waveform in the sense of using very simple matched filter to maximize the signal-to-noise ratio (SNR). Recently, Lyapunov exponent spectrum of the chaotic signals after being transmitted through a wireless channel has been shown to be unaltered, paving the way for wireless communication using chaos. In wireless communication systems, inter-symbol interference (ISI) caused by multipath propagation is one of the main obstacles to achieve high bit transmission rate and low bit error rate (BER). How to resist multipath effect is a fundamental problem in a chaos-based wireless communication system (CWCS). In this paper, implementation of a CWCS is presented. It is built to transmit chaotic signals generated by a hybrid dynamical system and then to filter the received signals by using the corresponding matched filter to decrease the noise effect and to detect the binary information. We find that the multipath effect can be effectively resisted by regrouping the return map of the received signal and by setting the corresponding threshold based on the available information. We show that the optimal threshold is a function of the channel parameters and of the transmitted information symbols. Practically, the channel parameters are time-variant, and the future information symbols are unavailable. In this case, a suboptimal threshold (SOT) is proposed, and the BER using the SOT is derived analytically. Simulation results show that the CWCS achieves a remarkable competitive performance even under inaccurate channel parameters.
• The fourth generation of the Sloan Digital Sky Survey (SDSS-IV) began observations in July 2014. It pursues three core programs: APOGEE-2, MaNGA, and eBOSS. In addition, eBOSS contains two major subprograms: TDSS and SPIDERS. This paper describes the first data release from SDSS-IV, Data Release 13 (DR13), which contains new data, reanalysis of existing data sets and, like all SDSS data releases, is inclusive of previously released data. DR13 makes publicly available 1390 spatially resolved integral field unit observations of nearby galaxies from MaNGA, the first data released from this survey. It includes new observations from eBOSS, completing SEQUELS. In addition to targeting galaxies and quasars, SEQUELS also targeted variability-selected objects from TDSS and X-ray selected objects from SPIDERS. DR13 includes new reductions of the SDSS-III BOSS data, improving the spectrophotometric calibration and redshift classification. DR13 releases new reductions of the APOGEE-1 data from SDSS-III, with abundances of elements not previously included and improved stellar parameters for dwarf stars and cooler stars. For the SDSS imaging data, DR13 provides new, more robust and precise photometric calibrations. Several value-added catalogs are being released in tandem with DR13, in particular target catalogs relevant for eBOSS, TDSS, and SPIDERS, and an updated red-clump catalog for APOGEE. This paper describes the location and format of the data now publicly available, as well as providing references to the important technical papers that describe the targeting, observing, and data reduction. The SDSS website, http://www.sdss.org, provides links to the data, tutorials and examples of data access, and extensive documentation of the reduction and analysis procedures. DR13 is the first of a scheduled set that will contain new data and analyses from the planned ~6-year operations of SDSS-IV.
• ### Localized Orbital Scaling Correction for Systematic Elimination of Delocalization Error in Density Functional Approximations(1707.00856)

The delocalization error of popular density functional approximations (DFAs) leads to diversified problems in present-day density functional theory calculations. For achieving a universal elimination of delocalization error, we develop a localized orbital scaling correction (LOSC) framework, which unifies our previously proposed global and local scaling approaches. The LOSC framework accurately characterizes the distributions of global and local fractional electrons, and is thus capable of correcting system energy, energy derivative and electron density in a self-consistent and size-consistent manner. The LOSC-DFAs lead to systematically improved results, including the dissociation of cationic species, the band gaps of molecules and polymer chains, the energy and density changes upon electron addition and removal, and photoemission spectra.
• ### Familia: An Open-Source Toolkit for Industrial Topic Modeling(1707.09823)

July 31, 2017 cs.CL, cs.IR
Familia is an open-source toolkit for pragmatic topic modeling in industry. Familia abstracts the utilities of topic modeling in industry as two paradigms: semantic representation and semantic matching. Efficient implementations of the two paradigms are made publicly available for the first time. Furthermore, we provide off-the-shelf topic models trained on large-scale industrial corpora, including Latent Dirichlet Allocation (LDA), SentenceLDA and Topical Word Embedding (TWE). We further describe typical applications which are successfully powered by topic modeling, in order to ease the confusions and difficulties of software engineers during topic model selection and utilization.
• ### Compressed sensing with local structure: uniform recovery guarantees for the sparsity in levels class(1601.01988)

June 2, 2017 math.NA, cs.IT, math.IT
In compressed sensing, it is often desirable to consider signals possessing additional structure beyond sparsity. One such structured signal model - which forms the focus of this paper - is the local sparsity in levels class. This class has recently found applications in problems such as compressive imaging, multi-sensor acquisition systems and sparse regularization in inverse problems. In this paper we present uniform recovery guarantees for this class when the measurement matrix corresponds to a subsampled isometry. We do this by establishing a variant of the standard restricted isometry property for sparse in levels vectors, known as the restricted isometry property in levels. Interestingly, besides the usual log factors, our uniform recovery guarantees are simpler and less stringent than existing nonuniform recovery guarantees. For the particular case of discrete Fourier sampling with Haar wavelet sparsity, a corollary of our main theorem yields a new recovery guarantee which improves over the current state-of-the-art.
• ### AstroCloud: A Distributed Cloud Computing and Application Platform for Astronomy(1701.05641)

Jan. 19, 2017 astro-ph.IM
Virtual Observatory (VO) is a data-intensively online astronomical research and education environment, which takes advantages of advanced information technologies to achieve seamless and global access to astronomical information. AstroCloud is a cyber-infrastructure for astronomy research initiated by Chinese Virtual Observatory (China-VO) project, and also a kind of physical distributed platform which integrates lots of tasks such as telescope access proposal management, data archiving, data quality control, data release and open access, cloud based data processing and analysis. It consists of five application channels, i.e. observation, data, tools, cloud and public and is acting as a full lifecycle management system and gateway for astronomical data and telescopes. Physically, the platform is hosted in six cities currently, i.e. Beijing, Nanjing, Shanghai, Kunming, Lijiang and Urumqi, and serving more than 17 thousand users. Achievements from international Virtual Observatories and Cloud Computing are adopted heavily. In the paper, backgrounds of the project, architecture, Cloud Computing environment, key features of the system, current status and future plans are introduced.
• ### A Joint Framework for Argumentative Text Analysis Incorporating Domain Knowledge(1701.05343)

Jan. 19, 2017 cs.CL
For argumentation mining, there are several sub-tasks such as argumentation component type classification, relation classification. Existing research tends to solve such sub-tasks separately, but ignore the close relation between them. In this paper, we present a joint framework incorporating logical relation between sub-tasks to improve the performance of argumentation structure generation. We design an objective function to combine the predictions from individual models for each sub-task and solve the problem with some constraints constructed from background knowledge. We evaluate our proposed model on two public corpora and the experiment results show that our model can outperform the baseline that uses a separate model significantly for each sub-task. Our model also shows advantages on component-related sub-tasks compared to a state-of-the-art joint model based on the evidence graph.
• ### PandaX-III: Searching for Neutrinoless Double Beta Decay with High Pressure $^{136}$Xe Gas Time Projection Chambers(1610.08883)

Oct. 28, 2016 hep-ex, nucl-ex, physics.ins-det
Searching for the Neutrinoless Double Beta Decay (NLDBD) is now regarded as the topmost promising technique to explore the nature of neutrinos after the discovery of neutrino masses in oscillation experiments. PandaX-III (Particle And Astrophysical Xenon Experiment III) will search for the NLDBD of $^{136}$Xe at the China Jin Ping underground Laboratory (CJPL). In the first phase of the experiment, a high pressure gas Time Projection Chamber (TPC) will contain 200 kg, 90% $^{136}$Xe enriched gas operated at 10 bar. Fine pitch micro-pattern gas detector (Microbulk Micromegas) will be used at both ends of the TPC for the charge readout with a cathode in the middle. Charge signals can be used to reconstruct tracks of NLDBD events and provide good energy and spatial resolution. The detector will be immersed in a large water tank to ensure $\sim$5 m of water shielding in all directions. The second phase, a ton-scale experiment, will consist of five TPCs in the same water tank, with improved energy resolution and better control over backgrounds.
• ### Ferroelectric domain inversion and its stability in lithium niobate thin film on insulator with different thicknesses(1604.05444)

April 19, 2016 cond-mat.mtrl-sci
Ferroelectric domain inversion and its effect on the stability of lithium niobate thin films on insulator (LNOI) are experimentally characterized. Two sets of specimens with different thicknesses varying from submicron to microns are selected. For micron thick samples (~28 um), domain structures are achieved by pulsed electric field poling with electrodes patterned via photolithography. No domain structure deterioration has been observed for a month as inspected using polarizing optical microscopy and etching. As for submicron (540 nm) films, large-area domain inversion is realized by scanning a biased conductive tip in a piezoelectric force microscope. A graphic processing method is taken to evaluate the domain retention. A domain life time of 25.0 h is obtained and possible mechanisms are discussed. Our study gives a direct reference for domain structure-related applications of LNOI, including guiding wave nonlinear frequency conversion, nonlinear wavefront tailoring, electro-optic modulation, and piezoelectric devices.
• ### Sparsity and Parallel Acquisition: Optimal Uniform and Nonuniform Recovery Guarantees(1603.08050)

March 25, 2016 cs.IT, math.IT, math.FA
The problem of multiple sensors simultaneously acquiring measurements of a single object can be found in many applications. In this paper, we present the optimal recovery guarantees for the recovery of compressible signals from multi-sensor measurements using compressed sensing. In the first half of the paper, we present both uniform and nonuniform recovery guarantees for the conventional sparse signal model in a so-called distinct sensing scenario. In the second half, using the so-called sparse and distributed signal model, we present nonuniform recovery guarantees which effectively broaden the class of sensing scenarios for which optimal recovery is possible, including to the so-called identical sampling scenario. To verify our recovery guarantees we provide several numerical results including phase transition curves and numerically-computed bounds.
• ### A review on locomotion robophysics: the study of movement at the intersection of robotics, soft matter and dynamical systems(1602.04712)

In this review we argue for the creation of a physics of moving systems -- a locomotion "robophysics" -- which we define as the pursuit of the discovery of principles of self generated motion. Robophysics can provide an important intellectual complement to the discipline of robotics, largely the domain of researchers from engineering and computer science. The essential idea is that we must complement study of complex robots in complex situations with systematic study of simplified robophysical devices in controlled laboratory settings and simplified theoretical models. We must thus use the methods of physics to examine successful and failed locomotion in simplified (abstracted) devices using parameter space exploration, systematic control, and techniques from dynamical systems. Using examples from our and other's research, we will discuss how such robophysical studies have begun to aid engineers in the creation of devices that begin to achieve life-like locomotor abilities on and within complex environments, have inspired interesting physics questions in low dimensional dynamical systems, geometric mechanics and soft matter physics, and have been useful to develop models for biological locomotion in complex terrain. The rapidly decreasing cost of constructing sophisticated robot models with easy access to significant computational power bodes well for scientists and engineers to engage in a discipline which can readily integrate experiment, theory and computation.
• ### The Improvement of Joule Balance NIM-1 and the Design of New Joule Balance NIM-2(1412.6895)

Dec. 22, 2014 physics.ins-det
The development of the joule balance method to measure the Planck constant, in support of the redefinition of the kilogram, has been going on at the National Institute of Metrology of China (NIM) since 2007. The first prototype has been built to check the feasibility of the principle. In 2011, the relative uncertainty of the Planck constant measurement at NIM is $7.7\times10^{-5}$. Self-heating of the coils, swing of the coil, are the main uncertainty contributions. Since 2012, some improvements have been made to reduce these uncertainties. The relative uncertainty of the joule balance is reduced to $7.2\times10^{-6}$ at present. The Planck constant measured with the joule balance is $h=6.6261041(470)\times10^{-34}$Js. The relative difference between the determined h and the CODATA2010 recommendation value is $5\times10^{-6}$. Further improvements are still being carried out on the NIM-1 apparatus. At the same time, the design and construction of a brand new and compact joule balance NIM-2 are also in progress and presented here.
• ### AsterixDB: A Scalable, Open Source BDMS(1407.0454)

July 2, 2014 cs.DB
AsterixDB is a new, full-function BDMS (Big Data Management System) with a feature set that distinguishes it from other platforms in today's open source Big Data ecosystem. Its features make it well-suited to applications like web data warehousing, social data storage and analysis, and other use cases related to Big Data. AsterixDB has a flexible NoSQL style data model; a query language that supports a wide range of queries; a scalable runtime; partitioned, LSM-based data storage and indexing (including B+-tree, R-tree, and text indexes); support for external as well as natively stored data; a rich set of built-in types; support for fuzzy, spatial, and temporal types and queries; a built-in notion of data feeds for ingestion of data; and transaction support akin to that of a NoSQL store. Development of AsterixDB began in 2009 and led to a mid-2013 initial open source release. This paper is the first complete description of the resulting open source AsterixDB system. Covered herein are the system's data model, its query language, and its software architecture. Also included are a summary of the current status of the project and a first glimpse into how AsterixDB performs when compared to alternative technologies, including a parallel relational DBMS, a popular NoSQL store, and a popular Hadoop-based SQL data analytics platform, for things that both technologies can do. Also included is a brief description of some initial trials that the system has undergone and the lessons learned (and plans laid) based on those early "customer" engagements.
• ### Digital Logarithmic Airborne Gamma Ray Spectrometer(1403.3147)

March 13, 2014 physics.ins-det
A new digital logarithmic airborne gamma ray spectrometer is designed in this study. The spectrometer adopts a high-speed and high-accuracy logarithmic amplifier (LOG114) to amplify the pulse signal logarithmically and to improve the utilization of the ADC dynamic range, because the low-energy pulse signal has a larger gain than the high-energy pulse signal. The spectrometer can clearly distinguish the photopeaks at 239, 352, 583, and 609keV in the low-energy spectral sections after the energy calibration. The photopeak energy resolution of 137Cs improves to 6.75% from the original 7.8%. Furthermore, the energy resolution of three photopeaks, namely, K, U, and Th, is maintained, and the overall stability of the energy spectrum is increased through potassium peak spectrum stabilization. Thus, effectively measuring energy from 20keV to 10MeV is possible.
• ### GreenMail: Reducing Email Service's Carbon Emission with Minimum Cost(1401.5546)

Jan. 22, 2014 cs.OH
Internet services contribute a large fraction of worldwide carbon emission nowadays, in a context of increasing number of companies tending to provide and more and more developers use Internet services. Noticeably, a trend is those service providers are trying to reduce their carbon emissions by utilizing on-site or off-site renewable energy in their datacenters in order to attract more customers. With such efforts have been paid, there are still some users who are aggressively calling for even cleaner Internet services. For example, over 500,000 Facebook users petitioned the social networking site to use renewable energy to power its datacenter. However, it seems impossible for such demand to be satisfied merely from the inside of those production datacenters, considering the transition cost and stability. Outside the existing Internet services, on the other hand, may easily set up a proxy service to attract those renewable-energy-sensitive users, by 1) using carbon neutral or even over-offsetting cloud instances to bridge the end user and traditional Internet services; and 2) estimating and offsetting the carbon emissions from the traditional Internet services. In our paper, we proposed GreenMail, which is a general IMAP proxy caching system that connects email users and traditional email services. GreenMail runs on green web hosts to cache users' emails on green cloud instances. Besides, it offsets the carbon emitted by traditional backend email services. With GreenMail, users could set a carbon emission constraint and use traditional email service without breaking any code modification of user side and email server side.
• ### The effect of limb kinematics on the speed of a legged robot on granular media(1303.7276)

Achieving effective locomotion on diverse terrestrial substrates can require subtle changes of limb kinematics. Biologically inspired legged robots (physical models of organisms) have shown impressive mobility on hard ground but suffer performance loss on unconsolidated granular materials like sand. Because comprehensive limb-ground interaction models are lacking, optimal gaits on complex yielding terrain have been determined empirically. To develop predictive models for legged devices and to provide hypotheses for biological locomotors, we systematically study the performance of SandBot, a small legged robot, on granular media as a function of gait parameters. High performance occurs only in a small region of parameter space. A previously introduced kinematic model of the robot combined with a new anisotropic granular penetration force law predicts the speed. Performance on granular media is maximized when gait parameters minimize body acceleration and limb interference, and utilize solidification features of granular media.
• ### Multi-functional foot use during running in the zebra-tailed lizard (Callisaurus draconoides)(1303.7292)

A diversity of animals that run on solid, level, flat, non-slip surfaces appear to bounce on their legs; elastic elements in the limbs can store and return energy during each step. The mechanics and energetics of running in natural terrain, particularly on surfaces that can yield and flow under stress, is less understood. The zebra-tailed lizard (Callisaurus draconoides), a small desert generalist with a large, elongate, tendinous hind foot, runs rapidly across a variety of natural substrates. We use high speed video to obtain detailed three-dimensional running kinematics on solid and granular surfaces to reveal how leg, foot, and substrate mechanics contribute to its high locomotor performance. Running at ~10 body length/s (~1 m/s), the center of mass oscillates like a spring-mass system on both substrates, with only 15% reduction in stride length on the granular surface. On the solid surface, a strut-spring model of the hind limb reveals that the hind foot saves about 40% of the mechanical work needed per step, significant for the lizard's small size. On the granular surface, a penetration force model and hypothesized subsurface foot rotation indicates that the hind foot paddles through fluidized granular medium, and that the energy lost per step during irreversible deformation of the substrate does not differ from the reduction in the mechanical energy of the center of mass. The upper hind leg muscles must perform three times as much mechanical work on the granular surface as on the solid surface to compensate for the greater energy lost within the foot and to the substrate.
• ### A Terradynamics of Legged Locomotion on Granular Media(1303.7065)

The theories of aero- and hydrodynamics predict animal movement and device design in air and water through the computation of lift, drag, and thrust forces. Although models of terrestrial legged locomotion have focused on interactions with solid ground, many animals move on substrates that flow in response to intrusion. However, locomotor-ground interaction models on such flowable ground are often unavailable. We developed a force model for arbitrarily-shaped legs and bodies moving freely in granular media, and used this "terradynamics" to predict a small legged robot's locomotion on granular media using various leg shapes and stride frequencies. Our study reveals a complex but generic dependence of stresses in granular media on intruder depth, orientation, and movement direction and gives insight into the effects of leg morphology and kinematics on movement.