abstract: The capability to sense and respond to external mechanical stimuli at various timescales is essential to many physiological aspects in plants, including self-protection, intake of nutrients, and reproduction. Remarkably, some plants have evolved the ability to react to mechanical stimuli within a few seconds despite a lack of muscles and nerves. The fast movements of plants in response to mechanical stimuli have long captured the curiosity of scientists and engineers, but the mechanisms behind these rapid thigmonastic movements still are not understood completely. In this article, we provide an overview of such thigmonastic movements in several representative plants, including *Dionaea*, *Utricularia*, *Aldrovanda*, *Drosera*, and *Mimosa*. In addition, we review a series of studies that present biomimetic structures inspired by fast moving plants. We hope that this article will shed light on the current status of research on the fast movements of plants and bioinspired structures and also promote interdisciplinary studies on both the fundamental mechanisms of plants' fast movements and biomimetic structures for engineering applications, such as artificial muscles, multi-stable structures, and bioinspired robots.
author: Q. Guo$^{1,2}$, E. Dai$^{3}$, X. Han$^{4}$, S. Xie$^{5}$, E. Chao$^{3}$, Z. Chen$^{4}$
date: 2024-10-03
title: Fast nastic motion of plants and bio-inspired structures

# Introduction

When one thinks of fast movement, plants do not usually come to mind. The movements of most plants typically involve a tropism, the growth or turning movement of a biological organism in response to a direction-dependent stimulus, and are normally only noticeable over the course of hours or days. Classic examples of tropisms include phototropic movements towards light, directional growth or reorientation in response to gravity (gravitropism), and touch sensitive growth of the roots (an example of thigmotropism) . In addition to tropisms, some carnivorous plants trap prey using nastic movements, or direction-independent, reversible orientational changes in response to direction-independent stimuli. Among these, of particular interest is the Venus flytrap (*Dionaea*) which, upon consecutive triggering of the sensitive hairs on its namesake trap, can close a trap in less than a fraction of a second to capture its arthropod prey. This unusual plant, living in bogs with nutrient-deficient soil, was first discovered in 1760 by the North Carolina colonial governor Arthur Dobbs and was considered as "the great wonder of the \[plant\] kingdom". Later, Charles Darwin systematically investigated the mechanisms of the Venus flytrap's fast movement and called the plant "one of the most wonderful in the world."

The movements of live plants are typcally classified into two categories, tropistic and nastic movements . Tropistic movements refer to directional movements in response to external stimuli that have certain directionality (e.g., gravity, light, and mechanical touches). For example, the development of tendrils may depend on where and how they touch an object that acts as a support, so this type of movement is tropistic (termed thigmotropic movement or thigmotropism, where the prefix 'thigmo' denotes touch responsive ). In contrast, nastic movements refer to direction-independent mechanical movements in response to external stimuli . For instance, the way that a *Dionaea* trap closes usually does not depend on where and in what direction the sensitive hairs are touched, so it is an example of a nastic (or more specifically, thigmonastic) movement. On the other hand, researchers have also distinguished between active and passive systems among plant movements . Action potentials are typically involved in active movements, whereby the live tissues or cells can actively deform in response to the stimuli by the transport of ions and fluids or the change in membrane permeability, often assisted by mechanical principles such as buckling, swelling, or cavitation . In contrast, passive movements refer to configurational changes of dead tissues that occur in response to environmental changes . Because of the vast amount of literature on plants' movement, we choose to focus primarily on fast, nastic movements and by "fast movement" we refer to any movement that takes place in few seconds or less.

Although it has been known for centuries that certain plants are capable of fast movement, the exact biomechanical, biochemical, and electrophysiological mechanisms which drive their fast movement have only begun to be comprehensively understood in recent years. For a more detailed overview of existing knowledge of these mechanisms, the readers may refer to some recent review articles on fast moving plants . It is also worth noting that plant biomechanics and mechanobiology have garnered increasing attention from the scientific community, as they represent "convergent paths to flourishing interdisciplinary research" . In the past decades, the role of mechanical forces, both internal and external, has been considered critical in the shaping of biological shapes, such as in stem cell differentiation , embryonic morphogenesis , cancer cell migration , plant morphogenesis , and touch-sensitive responses of plants . In this review, nevertheless, our focus is on the mechanically stimulated movement of plants.

The mechanical sensing, actuation, and movement in plants have since become sources of inspiration in biomimetic design, which have a wide range of engineering applications, including structural mechanics, biomedical engineering, and chemical engineering. In this article, we review the characteristics and mechanisms of several representative fast-moving plants and some recent development of engineering structures inspired by these plants. It is foreseeable that interdisciplinary collaborations between biologists, physicists, chemists, mathematicians, and engineers will bring forth new insights that ultimately lead to a unified picture of how these non-muscular systems operate. We hereby hope that our work can promote further interdisciplinary studies on the movements of plants as well as the advancement of bio-inspired and biomimetic technologies.

# Fast Movement of the Venus flytrap

In this section, we examine the important structures and features of the Venus flytrap and review the biomechanical, chemical, and electrophysiological mechanisms behind its fast movement. Further study of these mechanisms will contribute to more comprehensive understanding of the nastic movements of plants.

## The structural features of the Venus flytrap

The most spectacular feature of the Venus flytrap, is the "open mouth", a pair of leaves that traps arthropods to nourish the plant. As shown in Figure (a), the modified leaves of the Venus flytrap have two parts, the upper and lower leaf. The lower leaf, also known as the petiolus, has an expanded leaf-like structure . The upper leaf has three main features: a shell-like "open mouth" composed of a pair of symmetrical lobes, three to five (typically three) sensitive hairs (see Figure (b)) on the inner side of each lobe, and cilia, the numerous interlocked "teeth" located on the border of each lobe. When an external force is applied to the trigger hairs, the trap swiftly closes. Usually, two consecutive touches are needed to trigger this rapid closure, a mechanism that helps ensure the plant does not waste energy on capturing non-prey. This movement can also be triggered by other external stimuli, such as electricity or heat . In addition, upon mechanical stimulation, touch receptors, the small trichomes located on the outer leaves and stems of the Venus flytrap, generate action potentials which can sensitize the trigger hairs for contact .

## Biomechanics of fast movements in the Venus flytrap

The physical mechanisms behind the trap of the Venus flytrap are highly complex and remain incompletely understood. Through systematic observation and experimentation, Darwin discovered that the inner layers of the Venus flytrap's lobes contract during closure, while the entire leaf's shape changes from convex to concave. Following these observations, researchers have proposed that the trap is mechanically bistable, with two stable states, i.e., the open and closed configurations, and an intermediate state .

Following Darwin, several additional explanations have been proposed for the underlying mechanisms behind this trap closure. Brown proposed that trap closure is caused by the expansion of the outside of the lobes, and furthermore, that the opposite expansion of the inner surface of the lobes drives the reopening process of the Venus flytrap. Inspired by Darwin's hypotheses, Forterre *et al.* used modern experimental technology to examine the biomechanics of leaf closure . Using elasticity theory, they proposed that a dimensionless parameter, $\alpha = W^4 \kappa^2/h^2$ (where $W$ denotes the size of the lobe, $h$ is the thickness, and $\kappa$ is the curvature), controls bistability (see Figure (a)-(c)). To test this hypothesis, Forterre and coworkers drew arrays of sub-millimetric ultraviolet-fluorescent dots on the surface of the leaves and irradiated the plant to fluoresce the dots, with the movements of these dots recorded by high speed camera during trap closure. From these experiments, the researchers inferred that the main source of the fast and dramatic closure of the trap is the bistable, doubly curved structure of the leaves, which snap-buckle to reverse their Gaussian curvature upon closing. Moreover, the sudden change of intrinsic curvature along the direction perpendicular to the midrib was determined to be the main source of the driving force behind trap closure, verified by measuring the strain field upon closure. The experiments also showed that the speed of the lobe during closure increases as a function of $\alpha$ (Figure (d)).

While Forterre and coworkers demonstrated that the fast motion of the Venus flytrap involves a snap-buckling instability of the shell-like geometry of the lobes, the cell and tissue scale sources of this active movement remain poorly understood. To address these concerns, Colombani and Forterre studied the cell scale *in vivo* measurements of the pressure and other poroelastic properties using a microfluidic pressure probe. Their initial results shed light on the time course of cellular pressure and volume change during plant movement, as well as membrane permeability and Young's modulus. The determined poroelastic time of 20-150s, much too large to account for rapid closure of the trap, suggests that some other cellular property (e.g., pre-stress of mesophyll cells) rather than active, osmotically driven water transport between cells could be the the cell scale source of rapid trap closure. Additional studies are necessary to fully understand the exact nature of this source.

In addition, Markin *et al.* proposed a hydroelastic curvature model to interpret the mechanism of the active movement of the Venus Flytrap. Here, based on the assumption that the trap possesses curvature elasticity and consists of outer and inner hydraulic layers where different hydrostatic pressures can build up, the natural principal curvatures are determined by the hydraulic states of the two layers of the plant, which are in turn defined by their differing hydrostatic pressures. Similar to a bimetallic couple, this so-called bilayer couple can change shape rapidly in response to triggering of the sensitive hairs of the plant. This change is thought to occur through transportation of water via aquaporins between the two layers. Yang *et al.* further developed a mathematical model to simulate the closing and opening dynamics of the Venus flytrap from the nonlinear dynamics and control perspective. In a separate study, Joyeux presented an elasticity model of the Venus flytrap by accounting for the anisotropic nature of the strain field of the trap and assuming that the open geometry is the minimum energy state, which embodies improvement over the previous model developed by Joyeux and coworkers. More recently, Pandolfi *et al*. showed that this movement is intimately dependent on gravity, which has not been considered before; however, the detailed biophysical mechanism remains inconclusively explained and requires further examination.

## Different ways of triggering the Venus flytrap

The Venus flytrap can be triggered in a number of ways. In nature, the trigger mechanism is mechanical stimulation, in which two consecutive touches on the sensitive hairs within around 30 seconds lead to closure. These touches can be on the same hair or on separate hairs. Each touch generates an electrical signal, which transmits between the lobes and the midrib of the trap. In addition to the typical two touches, a single sustained touch to one trigger hair can also close the trap, if two electrical signals are generated about two seconds apart in the same manner . Furthermore, there exists a temperature dependence of mechanically stimulated closure, with evidence suggesting that higher temperatures result in a greater sensitivity to mechanical stimulation .

External electrical stimulation, involving the transmission of electrical charge between a lobe and the midrib by means of electrodes and a charged capacitor can also induce closure. At room temperature, a total charge of about 14 $\mu$C , applied at once or in smaller increments, also causes the Venus flytrap to close, though the exact amount of charge required increases with trap size and decreases at higher temperatures. Closure by repetitive application of smaller charges demonstrates electrical memory in the Venus flytrap (see the section *Electrophysiology behind the Venus flytrap's fast movement*), subject to a time interval similar to that associated with mechanical stimulation .

The Venus flytrap also closes in response to optical stimulation in the mid-infrared region. The snapping mechanism can be triggered with 7.35 $\mu$m radiation from a continuous-wave laser directed at the interior of a Venus flytrap. Furthermore, electrical signals are found to propagate in the Venus flytrap concurrently with each optical stimulation, in the same manner observed in mechanical stimulation . Also noteworthy are recent experiments involving chemical stimuli of trap closure. Volkov *et al.* showed that application of 10 $\mu$L of ether or chloroform on the midrib induces slow (10s) closure of the trap, and similarly, that application of 10 $\mu$L of chloroform directly on the sensitive hairs can induce closure as fast as mechanical stimulation of the hairs.

## Biochemistry behind the Venus flytrap's fast movement

A number of plant hormones have been identified and linked to several key properties of the Venus flytrap and its prey-catching and digesting mechanisms. In plants, lipid-based hormone signals known as jasmonates regulate defense mechanisms in response to abiotic and biotic stresses. In contrast, these jasmonates play a role in nutrion retrieval for the Venus flytrap: the phytohormone jasmonic acid and coronatine (the molecular mimic of the isoleucine of jasmonic acid) are jasmonates involved in the secretion of digestive enzymes. Perez *et al.* demonstrated that spray-application of these jasmonates on just one leaf can trigger secretion of digestive enzymes throughout a Venus flytrap without any mechanical stimulation, indicating that these jasmonates are systemic in nature . In addition, the researchers showed that jasmonates are involved in the trap closing mechanism. They also showed that isolated traps from *D. muscipula* can be induced to close when aqueous-phase extracts containing 12-oxo-phytodienoic acid are adsorbed on the leaves. In these experiments, the total quantity of bioactive substance, as opposed to the concentration, influences closure. One hypothesis suggests that this bioactive trap-closing chemical may be secreted stepwise in response to each action potential, inducing trap-leaf closure past a specific threshold . This hypothesis would thereby link the observed memory effect associated with the accumulation of stimuli to the accumulation of this chemical.

Abscisic acid (ABA) is also involved in the trap closing mechanism. ABA is a plant hormone that activates anion channels to close stomata, reducing transpiration. ABA is produced in plants in response to low water availability. In *D. muscipula*, ABA has been shown to lower trap sensitivity to mechanical stimuli. Compared to well-watered Venus flytraps, both ABA-sprayed and water-deficient traps require additional mechanical stimulation to trigger fast trap closure .

## Electrophysiology behind the Venus flytrap's fast movement

The Venus flytrap's rapid closure is typically triggered by consecutive mechanical stimuli within a window of 30 seconds on the one or multiple sensitive hairs. It is known that such mechanical stimulation activates certain ion channels to generate both a receptor potential and an action potential, which may be responsible for the opening of fluid channels between different hydraulic layers. Nevertheless, two consecutive stimuli only put the trap into a semi-closed state. Without any further stimulations, the trap will gradually reopen in a few hours. In the case where a trapped insect keeps struggling and stimulating the sensitive hairs, the trap will close tighter and tighter until it is fully sealed, upon which the digestive fluids are secreted to help digest the nutrients.

As mentioned before, an alternative way to trigger closure is by applying an electrical pulse above a threshold value between the upper layer of one lobe and the midrib. In such an experiment, the polarity is critical because an electrical pulse with inverted polarity does not trigger trap closure . Whether triggered mechanically or electrically, the rapid closure of the Venus flytrap always follows the generation of receptor and action potentials . There are usually three phases during closure: a mechanically silent phase (no noticeable movement), an accelerating movement phase, and a fast movement phase . Volkov and co-workers employed a high-speed data acquisition device to record the duration and amplitude of the action potentials and showed that the duration of an action potential in the Venus flytrap is typically about 1.5 milliseconds. Furthermore, the researchers used ion channel uncouplers and blockers to explore the mechanisms of different closure phases. They found that both ion channel uncouplers and blockers can increase the time delay of trap closure dramatically while decreasing the speed of closure.

It is also interesting to note that the Venus flytrap exhibits short-term electrical memory , i.e., separate electrical pulses, each of which is less than the threshold and within a time window of no more than 50 seconds, can initiate closure as long as the sum of electrical stimulus exceeds the threshold value. Volkov and coworkers demonstrated that this stimuli threshold is temperature dependent, with a charge threshold between 8 to 9 $\mu C$ (with a larger electrical threshold corresponding to larger traps) at room temperature compared to only 4.1 $\mu C$ at 28 to 36 degree Celsius . This electrical memory was also demonstrated to be time dependent through application of below-threshold electric shocks separated by a 15 second interval, with decreased number of below-threshold electric shocks required for closure of the trap with increasing current of each shock. Researchers have also begun to discover the biomolecular components underlying this electrical memory. Ueda and co-workers demonstrated that there exists a threshold of accumulated bio-metabolite for stimulating trap closure. Using bioassays to separate *Dionaea* extracts, a bioactive polysaccharide was identified that is capable of triggering trap closure in absence of any external mechanical or electrical stimuli. The molecular electronics and different electrical characteristics of the complete hunting cycle (the open state, closed state, locked state, constriction and digestion state, and semi-open state) were further examined . Putting the results from these studies together, it is promising that a more comprehensive understanding of the biological electronic system and the associated "memory" and actuation mechanisms in plants such as *Dionaea* might be achieved. Furthermore, while there are clear parallels to be drawn between the electrical signals in plants and animals, additional studies are required to elucidate the unique physiological origins and functions of these signals in plants .

# Comparative biomechanics on active plant traps

The mechanics of movements in *D. muscipula* differ from those in the underwater carnivorous traps *Aldrovanda vesiculosa* and *Utricularia*. Poppinga and Joyeux compared the biomechanical mechanisms of *Aldrovanda* and *Dionaea* by modeling the *Aldrovandas* trap as a pair of thin elastic shells hinged to the midrib (see Figure (c)-(g)). The midrib connecting the two lobes can bend inward in the closed configuration. This study indicates that the *Aldrovanda*'s trap closure is due to swelling and shrinking of tissues in or around the midrib, rather than mechanical buckling that takes place in *Dionaea muscipula*. Moreover, the model suggests that this mechanism is the source of the rapid speed of trap closure in *Aldrovanda*, i.e., a large opening or closing movement of the trap can result from a relatively small bending deformation of the midrib, in contrast with the snap-buckling of convex to concave curvature seen in the lobes of *Dionaea*. More recently, Joyeux proposed a precise mechanism for the action of motor cells surrounding the midrib of a *Aldrovanda* trap that could lead to closing and opening of the trap. This new model confirms that the reversible movements of *Aldrovanda* trap can be controlled by the swelling and shrinking of the motor cells near the midrib, in contrast to the irreversible mechanism operating in *Dionaea* traps where mechanical buckling plays a dominant role.

Bladderwort (*Utricularia*) traps also employ a mechanical buckling mechanism which differs from that of *Dionaea* . *Utricularia*'s suction traps are among the most complicated of all active traps in the plant kingdom . The trap is composed of a hollow trap body, a "trap door" (which grants entrance to the trap body's interior), and four outward facing trigger hairs located on the trap door . During the trap-setting phase, active cells pump water out of the trap body through two-armed glands, termed bifids, creating a negative pressure differential up to -16 kPa between the interior and exterior of the trap body. As a result, the thin walls of *Utricularia* become concave (modeled as spherical) with stored elastic energy in the doubly-curved shell of the trap body. Joyeux *et al.* modeled this buckling scenario like a trap door. The researchers also showed that this buckling can also be initiated by a local deformation of the shell from mechanical stimulation of the trigger hair. Figure further suggests that the swinging motion of the panel starts with local curvature inversion, and that the "door" acts like a flexible valve that flips under external pressure (and changes curvature from concave to convex). Additionally, *Utricularia* traps fire not only upon mechanical stimulation of the trigger hairs, but also spontaneously to capture small organisms and detritus (e.g., phytoplankton or bacteria) which are too small to mechanically trigger the trap . This spontaneous action occurs through mechanical oscillations of regular pumping and mechanical buckling of the trap door. The underlying mechanisms of the spontaneous firing could inspire biomimetic design of engineering structures that are capable of autonomous elastic deformation.
Furthermore, Llorens and co-workers developed a model to simulate the dynamical process of *Utricularia*'s trap closure by considering hydrodynamics, buckling of the door, and deformation of the trap body. This comprehensive treatment can predict the main mechanical features of the *Utricularia*'s trap closure, namely triggered closure and spontaneous firing. Firing and resetting of traps in *Utricularia* are associated with water flow and trap volume changes, which mechanism is studied experimentally by Lubomir Adamec. The experimental data suggest that only a direct recording of a physiological activity in trigger hairs of a stimulated *Utricularia* trap can determine whether or not the mechanical stimulus results in an action potential (like in *Dionaea muscipula* and *Aldrovanda vesiculosa*), which could serve as an electrical signal for the trap-door buckling . Bartosz J. Plachno *et al.* studied the relationship between trap anatomy and its firing-resetting efficiency in nine Australian *Utricularia* species and found that the trap firing-resetting rate, while conceivably dependent on trap thickness, did not seem to depend on the number of cell layers in the trap wall, as quantitative functional trap characteristics of the 3-4 cell layer thick traps of *U. volubilis* were similar to those of the 2-3 cell layer thick traps of *U. dichotoma* .

# Catapult-flypaper traps

Sundews (*Drosera*) are mostly known for capturing prey with sticky traps. However, one recently discovered sundew (*Drosera glanduligera*, *Droseraceae*) demonstrates fast motion for the capture of prey. Upon sensing prey with its touch-sensitive snap tentacle heads, *Drosera glanduligera* catapults prey using outstretched snap-tentacles (as shown in Figure (a)) onto adjacent sticky glue-tentacles (Figure (b) and (c)) in as little as 75 milliseconds . These sticky glue-tentacles then slowly fold up and surround the prey over several minutes to several hours while digesting the prey with secreted digestive enzymes. The rapid catapulting movement of the snap-tentacle is considered to be primarily due to change in hydraulic pressures within the tentacle. Two possible mechanisms of this hydraulic pressure change have been proposed : rapid transport of water from cells in the adaxial half of the tentacle to cells in the abaxial half, resulting in adaxial contraction and abaxial extension and the fast catapulting movement; or, a sudden loss of turgor pressure in adaxial cells, followed by tentacle bending due to prestress of the abaxial surface to release the elastic energy stored in the epidermal cells on the adaxial side. However, further investigations are needed to test these hypotheses and elucidate the operating mechanisms of the fast catapulting motion of the snap-tentacles. Notably, unlike the snap-through movement of the Venus flytrap, the catapulting movement of the snap-tentacles is irreversible, possibly because non-recoverable fracture of the epidermal cells occurs in the hinge zone of the snap-tentacle upon initiation of movement. Some other flypaper traps, such as *Pinguicula* (), less mobile in nature, are nevertheless not reviewed in detail in the current article.

# Other fast movements in plants

Besides these carnivorous plants with active traps, there are other plants that can move at a comparably dramatic speed. Since the main focus of this review article is on the fast, nastic movements mostly discovered in carnivorous plants, here we only mention a few examples of other types of fast movements in plants. Interested readers can refer to some excellent review articles . As yet another well-known moving plant, *Mimosa pudica* exhibits a rapid, defensive response to external stimuli, such as the closing of its leaves (Figure (a)-(c)) and bending of its pulvinus, the joint-like thickening of the plant near the bottom of a leaf or leaflet where the plants "motor cells" reside. (Figure (d)). These motor cells, divided into flexor and extensor cells on the ventral and dorsal side of the leaf, respectively, regulate the volume and shape according to their relative turgor pressure. Upon mechanical stimulation of the plant, an action potential is initiated which propagates a signal from the site of the stimulus to the pulvinus. This results in another action potential which initiates the differential transport of K$^{+}$ and Cl$^{-}$, leading to a change in turgor pressure in the cells (Figure (e)). Consequently, the water potential in the extensor cells increases, resulting in loss of water and shrinkage of the cells, with the opposite occurring in the flexor cells, which swell to a turgid state. As a result, the leaflet folds up, ensnaring the trapped prey within. Calcium ions are shown to play an important role in this nastic movement as well . More recently, Song *et al.* , using advanced bio-imaging techniques including 2D X-ray micro-imaging and 3D X-ray tomography, non-destructively monitored the pulvinis at a high resolution to analyze the morphological characteristics of the pulvinis in real time.

In contrast with the plant initiated, hydraulically driven movements of *Mimosa pudica*, some plants use purely mechanical mechanisms to achieve rapid movements. Among such plants, a variety of mechanisms exist to achieve the fast dispersion of small solid objects, such as spores, pollen, or seeds . For example, the stamen of the bunchberry dogwood (*Cornus canadensis*) has evolved a trebuchet-like method for "catapulting" their pollen skywards after the flower opens , and the sphagnum moss can shoot its spores out like an airgun using built-up internal pressure . Yet another example of projectile motion is the release of fern spores , where through a cavitation catapult mechanism, generation of bubbles causes rapid cell expansion, inducing dramatic bending of the annulus cells and subsequently converting stored elastic energy into kinetic energy .

Last but not least, a number of plants also exhibit thigmonastic movements in their reproductive parts such as the floral organs and fruits . For example, when touched gently, the stamens of the moss rose, *Portulaca grandiflora* can bend within less than a second and return to the original positions in about five minutes , similar to the thigmonastic movement of *Mimosa* leaves in terms of the response times. Orchids have also evolved remarkable actuating structures and morphologies, such as sexually dimorphic flowers, to induce insects to cross-pollinate for them . Such touch sensitive movements of the reproductive organs represent an active adaption to enhance cross-pollination by transporting the pollen grains to arthropods (such as bees) and birds who serve as pollinators. Interesting questions then naturally arise about how these plants and their pollinators have co-evolved over millions of years. Addressing the molecular and biophysical mechanisms in these thigmonastic movements may help elucidate the co-evolutionary relationships between the species of flowering plants and the pollinators as well as provide information about the evolution paths of these versatile floral morphologies .

# Engineered actuating structures inspired by fast-moving plants

Fast-moving plants, as examples of dynamically morphing biological structures, are good model systems in the field of biomimetics, a relatively new discipline that takes inspirations from biology for problems in engineering. For example, the leaves of the *Mimosa* plant, which contract upon certain external stimuli (physical, heat, chemical etc.), have inspired research in structures that can change on demand between flexible and stiff states through changes in hydraulic pressure . In particular, *Dionaea muscipula*, the Venus flytrap, has drawn strong interest from engineers and plant biologists alike. The quick sensing and closing properties of *Dionaea* have inspired a number of imitative designs with potential applications in drug delivery systems, artificial muscles, and shape changing structures .

Holmes and Crosby used mechanical buckling of plates (in a two-dimensional analog of Euler buckling) and bi-axial compressive loads to fabricate an array of lens-like bistable shells, shown in Figure . A layer of biaxially prestretched polydimethylsiloxane (PDMS), which is patterned with periodic arrays of holes, is topped via spin coating by a thin film of uncured PDMS. The two layers of PDMS are then bonded via crosslinking. The micro-lenses formed by this process can possess either a concave or convex shape; a characteristic snap-through movement, similar to that of the Venus flytrap, occurs upon transition between the shapes. The focal point of each micro-lens is either above the structure surface (when the micro-lens is convex) or below the structure surface (when the micro-lens is concave). The researchers demonstrated finely tunable mechanical properties of the micro-lens structure through tuning of the geometric properties of the structure, e.g., lens size and spacing between adjacent lenses.

Chen and coworkers studied the geometric and mechanical properties of a bistable strip (Figures. (a) to (e)). They proposed a theoretical model for the large deformation of shell structures by the conformation of a strip on the doubly-curved surface of a torus. The theory captures the main feature of the competition between bending and stretching energy of an elastic plate when the Gaussian curvature changes. In addition, they fabricated a physical model of a bistable structure by uniaxially pre-stretching two rubber sheets in perpendicular directions, "sandwiching" a much thicker rubber sheet between the two pre-stretched sheets and bonding the three sheets together. Using the theoretical model, the researchers identified two dimensionless parameters (related to the mechanical forces and geometry) that control the bistable behavior of the system. In particular, the geometric parameter, $\eta \equiv W \sqrt{\kappa/H}$, also discussed by Armon *et al.* in an independent study, is equivalent to $\alpha^{1/4}$ (where $\alpha$ is the geometrical parameter that controls the elastic energy barrier in the Venus flytrap's closure ). Their work quantified the conditions for bistability, providing a mathematical framework for the study of bistable and multi-stable structures, with application in artificial muscles, bio-inspired robots, and shape-changing structures. Subsequent finite element simulations and experiments further support the validity of the theory . Based on this mechanical principle, bistable helical ribbons have been designed which can morph from one helical shape to the other upon certain external stimulus .

Hydrogel polymers, which swell and contract in response to a wide range of input stimuli (chemical, temperature, mechanical), can lead to a variety of sensing and actuating systems with a wide range of engineering applications . Inspired by *Dionaea*, Lee *et al.* fabricated a jumping microgel device as illustrated in Figure (h), utilizing the sensing-actuating properties of hydrogel to engineer a microgel jumping device that generates a snapping motion within 12 milliseconds. The researchers created a doubly curved shape from hydrogel and employed elastic instability to create the snap-buckling mechanism responsible for the quick actuation of the device. Upon swelling, the microgel legs of the device snap-buckle, resulting in a rapid jump. The power density of this device, at 34.2 mW g$^{-1}$, approaches that of human muscle. Also inspired by *Dionaea*, Gracias and coworkers have designed and fabricated smart photo-patterned bilayer hydrogel actuators. These hydrogel actuators exhibit external stimuli dependent swelling behavior, which, employed in a bilayer design, can act as a hinge for a Venus flytrap-like actuator. Santulli and coworkers , combining hydrogel polymer with braided fibers, developed and characterized actuators with tunable stiffness and force generation capacity.

In order to understand the behavior of these hydrogel devices, the underlying behaviors and properties of the hydrogel must be characterized and quantified. Stoychev *et al.* explored the folding behavior of rectangular hydrogel bilayers on a substrate, characterizing the effects of both non-homogenous swelling and adhesion of the hydrogel to the substrate. The capillary forces which drive liquid absorption based hydrogel actuation have been reviewed by Roman and Bico . Nevertheless, many properties and behaviors of hydrogel structures remain to be explored.

Hydrogel composites allow for the design of complex, dynamically morphing pre-programmed shapes. Erb and coworkers , inspired by the shape changes in natural seed dispersal units and carnivorous plants, designed hydrogel composites with embedded structural reinforcements which can deform into complex, bioinspired shapes, such as a twisting seed pod. Hydrogels and other polymers may also be utilized with origami design principals, combining the sensing-actuating properties of polymers with the pre-programmed, dynamically tunable properties of origami designs. For instance, Wei *et al.* designed hybrid hydrogels which can dynamically morph into specific shapes, including helical sheets with reversible chirality, while Ryu and coworkers develped photo-origami designs which, in reaction to certain wavelengths of light, fold into precise 3-D structures from 2-D sheets . Origami designs incorporating stiff plate elements and biomimetic leaf folding patterns might also be used towards the development of novel actuators.

In addition, ionic polymer-metal composites (IPMCs) have garnered much interest from the scientific community for both their biomimetic sensor and actuator properties . IPMCs exhibit characteristic sensing and mechanical bending in an electrical field, actions which strongly resemble the sensing and bending capabilities of the leaves of the Venus flytrap. IPMCs also have potential medical and industrial applications; current study of IPMCs focuses on further modelling of IPMC actuators and development of application prototypes . Because of the similarities between the two, studies of the bending and sensing mechanics behind IPMC actuators and Venus flytrap mechanics go hand in hand. As illustrated in Figure (f) and (g), a novel flytrap-inspired robot was fabricated by Shahinpoor using IPMC artificial muscles as distributed sensors and actuators . A conductive copper spine served as a "midrib", which transmitted the signals sensed by the IPMC trigger hair to the solid-state relay system and to activate a small dynamic voltage generator capable of actuating the robotic trap. When certain stimuli were applied to the IPMC trigger hairs, the functioning of the voltage generator can cause the flytrap robot to close. Such biomimetic design of smart systems/structures, integrated with functional sensors and actuators, can be potentially useful in a number of engineering applications .

The hydraulic redistribution driven movements of plants have also influenced the development of osmotically actuating, shape changing structures. Sinibaldi **et al.** , inspired by osmotic-driven plant movements, modelled the dynamic behavior of two implementations of the osmotic actuator concept. Both models were formulated on the basis of two chambers, a reservoir chamber (RC) filled with solute and actuation chamber (AC) filled with solute and solvent, which are separated by a solute-impermeable osmotic membrane (OM). The two models are different in their implementation of the actuator: the first has actuation work stored through the elastic deformation of a spring by a piston which displaces through solvent flux; the second has actuation work stored through the elastic deformation of a non-permeable bulging membrane. These models represent preliminary steps towards understanding towards the development of biomimetic, biorobotic solutions inspired by plant movement.

There are additional examples of morphing structures inspired by the hydraulic redistribution driven plant movement. Inspired by the hydraulic movements of *D. muscipula* and *M. pudica*, Pagitz *et al.* developed a model for morphing structures driven by pressure-actuated cellular structures . The researchers further refined this model by introducing a shell-like structure that can, through its characteristic "snap-through" morphology, change the sign of its Gaussian curvature . These shape changing structures, driven by multi-layer cellular pressure actuation, have potential applications in aerospace engineering, automobile engineering, architectural design, etc. Moreover, Freeman designed biomimetic flexible morphing membranes, so called "engineered nastic membranes", for on-demand change of aircraft wing shape during flight. Upon introduction of adenine triphosphate (ATP) to these nastic membranes, proton pumps are activated which move select ions across the membrane into a cylindrical inclusion (analogous to protons being pumped into the intracellular space). The resulting electrochemical gradient activates co-transporters and exchangers located in the nastic membrane, which drive the ions back out of the inclusion and down their electrochemical gradient while concurrently using this ionic flux to move fluid into the inclusion. This inward flow of fluid results in expansion of the cylindrical inclusion until fluid pressures are equalized and hydrostatic equilibrium reached and is responsible for the dynamic morphology of the structure. Detailed analysis of the mechanical, morphological, and electrochemical properties of this bio-inspired actuator in ABAQUS for the surrounding polymer and UFLUID for the driving transport processes revealed that various properties including slow reaction time, difficulty in upwards scaling, and limited displacement prevented this design from application in its original purpose of morphing aircraft wings; however, the design showed promising applications in enhanced osmotic actuation, water purification through protein-driven selective transport and highly targeted drug or vaccine delivery devices, which burst and release their contents upon contact with the proper stimulus in the body. More recent developments to this design by Freeman and Leiland incorporated the use of voltage-gated ion channels for the control of activation and a novel dual activation chamber design capable of simultaneous swelling and shrinking, with potential applications in biomorphic structures.

The mechanical principles demonstrated by moving plants also have potential application in systems engineering. Recently, Lenau and Hesselberg studied the self-organization principles demonstrated in stimulated motion plants. They divided the plant motions into three categories (according to the purpose): dispersal of spores, pollen, or seeds (e.g., sphagnum mosses); obtainment of nutrients (e.g., carnivorous plants such as *Dionaea*); and protection of vital organs or obtainment of optimal light and temperature conditions (e.g. *Mimosa*). They outlined the self-organization principles behind sun-tracking plants, i.e., detect sunlight, communicate this signal to an actuator, and perform a motion to optimize internal and external conditions based on the amount of sunlight, and showed that these design principles could be applied to the design of window blinds that open and close based on the amount of sunlight or solar energy devices that track the sun's position through the sky in order to maximize exposure. These designs for sun tracking devices have potential applications in energy conservation on an industrial or consumer scale. Also inspired by the sun-tracking movements of plants, Dicker and coworkers developed a biomimetic analogue of the sun-tracking leaves of the Cornish Mallow (*Lavatera cretica*), which reversibly move the plane of the leaf to be perpendicular to available sunlight, termed diaheliotropic movement. As shown in Figure , the designed biomimetic photo-actuator consists of a reversible photobase as the sensor and pH sensitive hydrogels as the actuator. The photobase 4,4-bis(dimethylamino)triphenlmethane leucohydroxide (malachite green carbinol base) is stored in a pair of artificial leaf veins which limit light exposure via shades oriented 45 degrees to the tube surfaces and changes pH when exposed to UV light. Upon pH change, the pH sensitive hydrogel, an epoxy hydrogel described by Yoshioka and Calvert , reversibly swells within a pair of novel flexible matrix composite (FMC) tubes in a twisted double T joint configuration. Due to the twisted double T joint configuration of the tubes and limited light exposure of the artificial leaf veins, the biomimetic photo-actuator continues to bend until the light exposure of the photobase is maximized, analogous to the sun-tracking movements of plants. This biomimetic sun-tracking device acts as a proof of concept for chemical sensing, control and actuation with potential applications in solar tracking for solar power generation.

# Conclusion

Plants with fast thigmonastic movements have evolved the capability to quickly and dramatically alter their shapes in response to certain external stimuli. Although these mechanically-induced active movements may appear to serve different purposes, from a protective mechanism (e.g., the folding of $Mimosa$ leaflets) to a hunting function (for example, the rapid closure of the traps of $\emph{Dionaea}$, $\emph{Utricularia}$, and $\emph{Aldrovanda}$), some common ground exists. For instance, recent investigations have started to show that most thigmonastic movements result from differential changes of turgor pressure in the tissues, and the opening of certain key ion channels are often involved . While some biophysical and biochemical mechanisms behind the fast motions of plants have been extensively studied, many questions remain unanswered. For example, how the action potentials are generated through the touch receptors in $Dionaea$ remains to be investigated. Moreover, the source and role of action potentials in fast-moving plants, especially in carnivores other than $Dionaea$, are not completely understood. The coupling between these action potentials and molecular-level activities (which drive macroscopic movement) also remains not well understood. In $\emph{Dionaea}$, recent studies suggest that it is the speed of ion transport, instead of water transport, that constrains the speed of osmotic flows; therefore, alternative mechanisms of perturbing the water potential in cells should be searched for . Last but not least, the roles of actin cytoskeleton, microtublues, and auxin during these active movements require elucidation.

Understanding the molecular, biomechanical, biochemical, and electrophysical mechanisms among the diverse types of movements in plants may be the key to reveal the universal mechanisms governing plants' movements and how they evolve as a result of adaption to the ever-changing living environments . However, it remains a grand challenge how to develop more realistic, multi-scale biomechanical models that incorporate various factors at the molecular, cellular, tissue, and organism level. The difficulty in bridging the physical and biological aspects of the problem also arises from the fact that there is a lack of effective non-invasive experimental techniques for measuring the stresses inside cells and monitoring the transport of fluids (not to mention ions) inside the plant cells during these rapid movements . The recent *in vivo* measurement of pressures and poroelastic properties in cells of $Dionaea$ and the measurement of mechanical stresses in developing embryos seem to have paved way for further breakthroughs in investigating the underlying mechanisms in fast-moving plants. In addition, the comparisons between the behaviors and operating mechanisms in the moving plants, $Mimosa$, $Dionaea$, $Aldrovanda$, $Utricularia$, and $Drosera$ also bring forth interesting questions from the evolutionary biology and comparative biomechanics perspectives that remain to be addressed. The study of fast-moving plants will not only foster integration of physics, chemistry, biomechanics, and molecular biology to address intriguing mechanisms involved in plants' fast motions, but also inspire various kinds of biomimetic designs for engineering applications.

# Author Contributions

Q.G., E.D., X. H., and Z.C. generated the idea for the manuscript, performed literature searches, helped to draft and edit it. S.X. and E.C. also performed literature searches, helped to draft and edit it.

# Acknowledgments

The authors thank J. Hanlin, I. Trase, J. Ballard, K. Tan, H. Zheng, and the anonymous reviewers for helpful comments.

# Funding statement

This work is in part supported by the National Science Foundation of China Grant No.11102040, Projects of International Cooperation and Exchanges NSFC (Grant No.11201001044), and Fujian University of Technology Research Fund (Grant No.GY-Z15006). Z.C. acknowledges support from the Society in Science - Branco Weiss fellowship, administered by ETH Z$\ddot{u}$rich.

# Conflict of interests

The authors have no financial or personal relationships with other people or organizations that might inappropriately influence or bias the work.


address:  Dept. of Epidemiology & Biostatistics
  VU University medical center Amsterdam
  Amsterdam
  The Netherlands
author: Carel F.W. Peeters
title: Bayesian Constrained-Model Selection
  for Factor Analytic Modeling

My dissertation (Peeters, 2012a) revolves around Bayesian approaches towards constrained statistical inference in the factor analysis (FA) model. Two interconnected types of restricted-model selection are considered. These types have a natural connection to selection problems in the exploratory FA (EFA) and confirmatory FA (CFA) model and are termed Type I and Type II model selection. Type I constrained-model selection is taken to mean the determination of the appropriate dimensionality of a model. This type of constrained-model selection connects with EFA in the sense of selecting the optimal dimensionality of the latent vector. Type II model selection is taken to mean the determination of appropriate inequality, order or shape restrictions on the parameter space. The dissertation connects Type II constrained-model selection to CFA by focusing on the determination of linear inequality constraints as expressions of the direction and (relative) strength of factor loadings. The figures accompanying this article are taken from the slides of my Division 5 Awards Symposium Invited address at the APA 2015 Annual Convention in Toronto. These slides can be retrieved from <https://github.com/CFWP/ConventionTalk>.

# Summary of Research

Three research aims guide the dissertation. These research aims are motivated by the potential of the connection between Bayesian model selection (by the Bayes factor) and constrained statistical inference for factor analytic modeling. Classical approaches towards analyzing constrained hypotheses are restricted to a limited class of constraints and models. The Bayesian approach is more flexible. To utilize this flexibility, certain problems need to be overcome.

The main problem of the Bayes factor (Kass & Raftery, 1995) for Type I model selection is that both the use of improper noninformative and proper but vague priors will yield indeterminate answers when the models to be compared are of differing dimension (Jeffreys, 1961). This is undesirable as especially under default prior choices the interpretation of the Bayes factor as a weighted likelihood ratio is warranted (Berger & Pericchi, 2004). This problem reflects on Bayesian approaches towards EFA and the selection of the dimensionality of the latent vector. Most approaches use conjugate priors, usually in conjunction with a triangularity condition for rotational determination. Even when the conjugate priors are relatively uninformative (vague), the triangularity conditions induce strong prior information, calling into question the exploratory nature of Bayesian EFA. Research Aim 1 then intends to allow for the use of improper priors in Bayes factor computation for Type I constrained-model selection, allowing for formal dimensionality selection in EFA with all the benefits of the Bayesian machinery but without inducing strong prior information.

**Research Aim 1**. *To construct a conceptually and computationally simple Bayes factor for Type I constrained-model selection that is determinate under usage of improper priors. Subsequently, this Bayes factor is to be embedded within a strategy towards a truly Bayesian EFA concerned with the selection of an optimal dimensionality for the latent vector.*

CFA seeks to incorporate theory into the factor model by imposing certain restrictions. These restrictions are meant to emulate sparse structure and are usually of the fixed-value equality kind with a focus on exclusion restrictions. Such rigidity in parameter specification may pose several problems. First, it implies a loss of information in the sense that more exclusion restrictions are usually applied than is necessary for identification of the FA model. Also, exclusion restrictions may amount to errors of omission, may make the unrealistic assumption that items are factorially pure, and may induce bias in estimates of the free parameters (Ferrando & Lorenzo-Seva, 2000; van Prooijen & van der Kloot, 2001). These issues are intricately connected to the well-known and widespread situation of exploratively obtained factor structures not being confirmed by CFA. Moreover, researchers in substantive fields usually have informed ideas regarding direction and magnitude of parameter effects that cannot be expressed using exclusion restrictions. What is wished for then, is expressions of factor structure, not through usage of exclusion restrictions in the matrix of factor loadings, but by the imposition of inequality constraints. This desire is formulated in the second guiding research aim.

**Research Aim 2**. *To construct a conceptually and computationally simple Bayes factor for Type II constrained-model selection that is geared towards inequalities on regression-type parameters. Subsequently, this Bayes factor is to be embedded within a strategy that specifies factor analytic structure using inequality constraints rather than through exclusion restrictions.*

A gap in FA practice is that there is no unequivocal strategy for integrating EFA and CFA. While the *modi operandi* are often viewed as distinct, it might be fruitful to view EFA and CFA as complementary techniques. For example, in the CFA model all attention regarding misspecification is geared towards the pre-specified pattern of factor loadings. The evaluation of model fit in CFA is then essentially the evaluation of a diffuse hypothesis (Hoyle & Duvall, 2004), as it is unclear in case of misspecification if the pattern of loadings or the factor dimensionality is to blame. The third research aim proposes a factor analytic strategy seeking to integrate EFA and CFA in order to avoid embarking on the evaluation of diffuse hypotheses.

**Research Aim 3**. *To let the provisions from Aims 1 and 2 conjoin in order to develop an integrative factor analytic strategy that articulates the complimentary nature of EFA and CFA.*

Part I of the dissertation covers Research Aims 1 and 2. The focus is on statistical and computational issues regarding formulation, prior selection, parameter estimation and model selection for the constrained Bayesian factor model. Part II comprises Research Aim 3 and purports that researchers often have competing theories that can be translated into inequality-constrained factor analytic models.

# Part I: Statistical & Computational Modeling

Part I starts by reviewing Markov chain Monte Carlo (MCMC) computation of the marginal likelihood. The Bayes factor consists of a ratio of marginal likelihoods. The candidate estimator method for marginal likelihood computation (Besag, 1989; Chib, 1995) is adapted to deal with (i) improper noninformative priors and (ii) the existence of (well-separated) symmetric posterior modes due to permutative invariance over the parameter indices, such that the ensuing Bayes factor is still determinate. Pending certain conditions, the provisions provide for what can be seen as a simulation consistent MCMC implementation of well-known default Bayes factors (Berger & Pericchi, 1996). This automated candidate estimator is subsequently applied to latent factor dimensionality selection in EFA. It is shown that a failure to abide certain regularity conditions on the FA model may result in violation of a crucial regularity condition for simulation consistency of estimates stemming from MCMC sampling. This implies that Bayesian approaches towards factor analytic dimensionality selection may be afflicted by some of the same regularity conditions that hamper classical approaches (see Figure ). The dissertation proposes an assessment strategy that ensures abidance of the regularity condition for simulation consistency such that the automated candidate estimator provides an appropriate stopping rule for factor analytic data compression. In passing, a truly Bayesian EFA is proposed.

Part I subsequently deals with a set of conditions for rotational identification of the oblique factor solution under usage of fixed zero elements in the factor loadings matrix (see also Peeters, 2012). It is shown that the well-known conditions for the oblique factor correlation structure (Jöreskog, 1979) need to be amended in order to obtain global rotational uniqueness. The amended condition set provides a way to design an unrestricted solution to the (Bayesian) CFA model. Unrestricted solutions correspond to EFA: An unrestricted confirmatory factor model (UCFM) is a FA model that places only minimal restrictions on the model parameters for achieving global rotational uniqueness of the factor solution. The restrictions are (in contrast to EFA) chosen so that they convey preconceived theoretical meaning and thus eliminate the need for post-hoc rotation of the solution for interpretation purposes.

The UCFM is pivotal in designing a Bayesian framework that takes parameter restrictions in the context of CFA beyond exclusion restrictions, by allowing inequality and approximate equality constraints to express substantive theoretical ideas regarding direction and magnitude of effect of factor loadings. This framework first requires the development of a Bayes factor for Type II constrained-model selection. Second, a strategy is developed for the demarcation of competing inequality-constrained formulations of factor analytic correlation structure. The strategy consists of choosing as a base model a UCFM. Substantive theory is then not represented by structural exclusions to express a pre-specified loading pattern, but by imposing inequalities on and between the free parameters in the loadings matrix (see Figure ). It is shown that when (i) proper but noninformative priors are chosen that are flat on the parameter space of the parameters on which inequalities are placed; and (ii) all competing inequality-constrained models are subsets of the UCFM; then the ensuing Bayes factor is determinate, its complexity is well-defined, and its computation is greatly simplified. Under this framework model fit and model complexity are explicitly connected to, respectively, the posterior and prior probability mass satisfying the constraints that define the constrained model (see Figure ).

# Part II: Applications

Part II conjoins the developments from Part I in an alternative factor analytic strategy merging EFA and CFA. This strategy lets EFA precede inequality-constrained CFA efforts by making it part of a total inferential procedure involving the selection of an optimal factor dimension before competing confirmatory structures are assessed. The strategy consists of the following steps:

1.  Embark on evaluating a series of unrestricted (EFA) models with respect to their factor dimensionality;

2.  Once the latent factor dimensionality is settled, specify a UCFM;

3.  Formulate, using the UCFM as a base model, competing inequality-constrained factor structures making use of a system of inequality constraints on and between the free parameters in the loadings matrix;

4.  Compute for each constrained model the Type II constrained-model selection Bayes factor and determine the constrained model best supported by the data.

This alternative strategy is brought to bear on research regarding the metabolic syndrome (see Chapter 5 of the dissertation and Peeters, Dziura, & van Wesel, 2014). The conception metabolic syndrome (MBS) refers to a clustering of interrelated risk factors of metabolic origin (Unwin, 2006). MBS is thought to be a precursor for the development of type 2 diabetes mellitus and atherosclerotic cardiovascular disease. The syndrome has social, behavioral, educational and psychological components, as MBS might result from maladaptive human metabolism in the face of food energy abundance in combination with a sedentary lifestyle (Wilkin & Voss, 2004). A high-profile MBS data set with anthropometric measurements on overweight and obese children and adolescents is reanalyzed using the alternative strategy. The alternative strategy is able to extract information with higher specificity. The findings may give inroads for behavioral and educational approaches towards MBS prevention.

# Contributions and Relevance

The contributions of the dissertation are connected to the research aims stated earlier. With the fulfillment of Aim 1 a generic computational procedure for Type I Bayes factors is available that is determinate under usage of default priors. This subsequently allows for improper prior usage in Bayesian EFA, thus respecting its exploratory nature. The developments under Aim 1 also spur learning on some lesser known indeterminacies in the factor model and their interrelationships with computational approaches towards dimensionality selection. These results also hold importance for non-Bayesian approaches towards factor analytic dimensionality selection. They imply that for informed decisions regarding factor dimensionality, likelihood ratio and information theoretic approaches benefit from a complete exploration of the likelihood, which can be achieved by objective Bayesian methods.

Fulfilling Aim 2 extends Bayesian model selection efforts regarding Type II model selection (e.g., Klugkist & Hoijtink, 2007; Mulder, Hoijtink, & Klugkist, 2010) and adds to the body of literature regarding inequality-constrained inference on regression-type parameters. More specifically, this research aim holds importance for the theory and technique of CFA. The formulation and development of inequality-constrained models allows for theoretically meaningful constrained coefficients beyond those needed simply to identify a model. This allows for higher specificity in model formulation as the direction and (relative) strength of factor loadings can be included in a formal model selection procedure. This may lead one to extract more information from a single analysis. The developments emphasize a break with simple structure models and may be viewed as constituting an alternative take on CFA.

Aim 3 answers the call for the integration of EFA and CFA (Steiger, 1994). The integrative strategy that is proposed connects EFA and CFA as complementary techniques that may be part of a total inferential procedure that aims to avoid evaluating diffuse model hypotheses. While Aims 1 and 2 are mostly of theoretical relevance, Aim 3 is thus of practical importance. The usefulness of the technical developments as merged in the integrative strategy is exemplified by bringing this strategy to bear on a published data set concerning risk factors for MBS. With the integrative strategy, it is shown that more and unexpected information can be extracted from these data.

# Acknowledgements

I would like to thank the APA Division 5 Award Committee for bestowing me with the 2015 Anne Anastasi Distinguished Dissertation Award. I express my gratitude towards my supervisors, Peter G.M. van der Heijden and Herbert Hoijtink. I am also indebted to the people who provided data for the applications in the dissertation: Bill Hopkins, Karin Lasthuizen, Peter Esaiasson, and James Dziura. Lastly, I acknowledge Floryt van Wesel who was, among many things, a collaborator on the application-oriented chapters of the dissertation.

This research was supported by grant NWO-VICI-453-05-002 of the Netherlands Organization for Scientific Research (NWO). It is based on the first chapter of my unpublished PhD dissertation (2012a). This version is a preprint of: Peeters, C.F.W. (2016). Bayesian Constrained-Model Selection for Factor Analytic Modeling. *The Score*, April Issue (<http://www.apadivisions.org/division-5/publications/score/2016/04/index.aspx>). David Herzberg's editorial help in improving this manuscript is gratefully acknowledged.