Search

For $q = p^n$ with $p$ an odd prime, the projective linear group $PGL(2,q)$ can be seen as the stabilizer of a conic $O$ in a projective plane $\pi = PG(2,q)$. In that setting, involutions of $PGL(2,q)$ correspond bijectively to points of $\pi$ not in $O$. Triples of involutions $\{ \alpha_P,\alpha_Q,\alpha_R \}$ of $PGL(2,q)$ can then be seen also as triples of points $\{P,Q,R\}$ of $\pi$. We investigate the interplay between algebraic properties of the group $H = \langle \alpha_P,\alpha_Q,\alpha_R \rangle$ generated by three involutions and geometric properties of the triple of points $\{P,Q,R\}$. In particular, we show that the coset geometry $\Gamma = \Gamma(H,(H_0,H_1,H_2))$, where $H_0 = \langle \alpha_Q,\alpha_R \rangle, H_1 = \langle \alpha_P,\alpha_R \rangle$ and $H_2 = \langle \alpha_P,\alpha_Q \rangle$ is a regular hypertope if and only if $\{P,Q,R\}$ is a strongly non self-polar triangle, a terminology we introduce. This entirely characterizes hypertopes of rank $3$ with automorphism group a subgroup of $PGL(2,q)$. As a corollary, we obtain the existence of hypertopes of rank $3$ with non linear diagrams and with automorphism group $PGL(2,q)$, for any $q = p^n$ with $p$ an odd prime. We also study in more details the case where the triangle $\{P,Q,R\}$ is tangent to $O$., Comment: 17 pages

Given a residually connected incidence geometry $\Gamma$ that satisfies two conditions, denoted $(B_1)$ and $(B_2)$, we construct a new geometry $H(\Gamma)$ with properties similar to those of $\Gamma$. This new geometry $H(\Gamma)$ is inspired by a construction of Percsy, Percsy and Leemans [1]. We show how $H(\Gamma)$ relates to the classical halving operation on polytopes, allowing us to generalize the halving operation to a broader class of geometries, that we call non-degenerate leaf hypertopes. Finally, we apply this generalization to cubic toroids in order to generate new examples of regular hypertopes., Comment: 25 pages

Recently, Leemans and Stokes constructed an infinite family of incidence geometries admitting trialities but no dualities from the groups PSL(2,q) (where $q=p^{3n}$ with $p$ a prime and $n>0$ a positive integer). Unfortunately these geometries are not flag transitive. In this paper, we construct the first infinite family of incidence geometries of rank three that are flag transitive and have trialities but no dualities. These geometries are constructed using chamber systems of Suzuki groups Sz(q) (where $q=2^{2e+1}$ with $e$ a positive integer and $2e+1$ is divisible by 3) and the trialities come from field automorphisms. We also construct an infinite family of regular hypermaps with automorphism group Sz(q) that admit trialities but no dualities., Comment: 14 pages

Dynamic compression of iron to Earth-core conditions is one of the few ways to gather important elastic and transport properties needed to uncover key mechanisms surrounding the geodynamo effect. Herein a new machine-learned ab-initio derived molecular-spin dynamics (MSD) methodology with explicit treatment for longitudinal spin-fluctuations is utilized to probe the dynamic phase-diagram of iron. This framework uniquely enables an accurate resolution of the phase-transition kinetics and Earth-core elastic properties, as highlighted by compressional wave velocity and adiabatic bulk moduli measurements. In addition, a unique coupling of MSD with time-dependent density functional theory enables gauging electronic transport properties, critically important for resolving geodynamo dynamics., Comment: 3 Figures in main document, 8 Figures in the supplemental information

We present a study on the transport and materials properties of aluminum spanning from ambient to warm dense matter conditions using a machine-learned interatomic potential (ML-IAP). Prior research has utilized ML-IAPs to simulate phenomena in warm dense matter, but these potentials have often been calibrated for a narrow range of temperature and pressures. In contrast, we train a single ML-IAP over a wide range of temperatures, using density functional theory molecular dynamics (DFT-MD) data. Our approach overcomes computational limitations of DFT-MD simulations, enabling us to study transport and materials properties of matter at higher temperatures and longer time scales. We demonstrate the ML-IAP transferability across a wide range of temperatures using molecular-dynamics (MD) by examining the thermal conductivity, diffusion coefficient, viscosity, sound velocity, and ion-ion structure factor of aluminum up to about 60,000 K, where we find good agreement with previous theoretical data.

We introduce the notion of moving absolute geometry of a geometry with triality and show that, in the classical case where the triality is of type $(I_\sigma)$ and the absolute geometry is a generalized hexagon, the moving absolute geometry also gives interesting flag-transitive geometries with Buekenhout diagram with parameters $(d_p, g, d_L) = (5, 3, 6)$ for the groups $G_2(k)$ and $^3D_4(k)$, for any integer $k \geq 2$. We also classify the classical absolute geometries for geometries with trialities but no dualities coming from maps of Class III with automorphism group $L_2(q^3)$, where $q$ is a power of a prime. We then investigate the moving absolute geometries for these geometries, illustrating their interest in this case.

We determine the Thurston unit ball of a family of $n$-chained link, denoted by $C(n,p)$, where $n$ is the number of link components and $p$ is the number of twists. When $p$ is strictly positive, we prove that the Thurston unit ball for $C(n,p)$ is an $n$-dimensional cocube, for arbitrary $n$. Moreover, we clarify the condition for which $C(n,p)$ is fibered and find at least one fibered face for any $p$. Finally we provide the Teichm\"uller polynomial for the face of Thurston unit ball of $C(n, -2)$ with $n\geq 3$., Comment: 25 pages, 15 figures. Comments are welcome

Tungsten (W) is a material of choice for the divertor material due to its high melting temperature, thermal conductivity, and sputtering threshold. However, W has a very high brittle-to-ductile transition temperature and at fusion reactor temperatures ($\geq$1000K) may undergo recrystallization and grain growth. Dispersion-strengthening W with zirconium carbide (ZrC) can improve ductility and limit grain growth, but much of the effects of the dispersoids on microstructural evolution and thermomechanical properties at high temperature are still unknown. We present a machine learned Spectral Neighbor Analysis Potential (SNAP) for W-ZrC that can now be used to study these materials. In order to construct a potential suitable for large-scale atomistic simulations at fusion reactor temperatures, it is necessary to train on ab initio data generated for a diverse set of structures, chemical environments, and temperatures. Further accuracy and stability tests of the potential were achieved using objective functions for both material properties and high temperature stability. Validation of lattice parameters, surface energies, bulk moduli, and thermal expansion is confirmed on the optimized potential. Tensile tests of W/ZrC bicrystals show that while the W(110)-ZrC(111) C-terminated bicrystal has the highest ultimate tensile strength (UTS) at room temperature, observed strength decreases with increasing temperature. At 2500K, the terminating C layer diffuses into the W, resulting in a weaker W-Zr interface. Meanwhile, the W(110)-ZrC(111) Zr-terminated bicrystal has the highest UTS at 2500K.

We use spin-lattice dynamics simulations to study the possibility of modeling the magnetic hysteresis behavior of a ferromagnetic material. The temporal evolution of the magnetic and mechanical degrees of freedom is obtained through a set of two coupled Langevin equations. Hysteresis loops are calculated for different angles between the external field and the magnetocrystalline anisotropy axes. The influence of several relevant parameters is studied, including the field frequency, magnetic damping, magnetic anisotropy (magnitude and type), magnetic exchange, and system size. The role played by a moving lattice is also discussed. For a perfect bulk ferromagnetic system we find that, at low temperatures, the exchange and lattice dynamics barely affect the loops, while the field frequency and magnetic damping have a large effect on it. The influence of the anisotropy magnitude and symmetry are found to follow the expected behavior. We show that a careful choice of simulation parameters allows for an excellent agreement between the spin-lattice dynamics measurements and the paradigmatic Stoner-Wohlfarth model. Furthermore, we extend this analysis to intermediate and high temperatures for the perfect bulk system and for spherical nanoparticles, with and without defects, reaching values close to the Curie temperature. In this temperature range, we find that lattice dynamics has a greater role on the magnetic behavior, especially in the evolution of the defective samples. The present study opens the possibility for more accurate inclusion of lattice defects and thermal effects in hysteresis simulations

In this work, we leverage atomistic spin-lattice simulations to examine how magnetic interactions impact the propagation of sound waves through a ferromagnetic material. To achieve this, we characterize the sound wave velocity in BCC iron, a prototypical ferromagnetic material, using three different approaches that are based on the oscillations of kinetic energy, finite-displacement derived forces, and corrections to the elastic constants, respectively. Successfully applying these methods within the spin-lattice framework, we find good agreement with the Simon effect including high order terms. In analogy to experiments, morphic coefficients associated with the transverse and longitudinal waves propagating along the [001] direction are extracted from fits to the fractional change in velocity data. The present efforts represent an advancement in magnetoelastic modelling capabilities which can expedite the design of future magneto-acoustic devices.

We investigate the relation between the topological entropy of pseudo-Anosov maps on surfaces with punctures and the rank of the first homology of their mapping tori. On the surface $S$ of genus $g$ with $n$ punctures, we show that the entropy of a pseudo-Anosov map is bounded from above by $\dfrac{(k+1)\log(k+3)}{|\chi(S)|}$ up to a constant multiple when the rank of the first homology of the mapping torus is $k+1$ and $k, g, n$ satisfy a certain assumption. This is a partial generalization of precedent works of Tsai and Agol-Leininger-Margalit., Comment: 20 pages, 8 figures, comments are welcome

Given a regular covering map $\varphi:\Lambda \to \Gamma$ of graphs, we investigate the subgroup $\operatorname{LAut}(\varphi)$ of the automorphism group $\operatorname{Aut}(A_\Gamma)$ of the right-angled Artin group $A_\Gamma$. This subgroup comprises all automorphisms that can be lifted to automorphisms of $A_\Lambda$. We first show that $\operatorname{LAut}(\varphi)$ is generated by a finite subset of Laurence's elementary automorphisms. For the subgroup $\operatorname{FAut}(\varphi)$ of $\operatorname{Aut}(A_\Lambda)$, which consists of lifts of automorphisms in $\operatorname{LAut}(\varphi)$, there exists a natural homomorphism $\operatorname{FAut}(\varphi)\to\operatorname{LAut}(\varphi)$ induced by $\varphi$. We then show that the kernel of this homomorphism is virtually a subgroup of the Torelli subgroup $\operatorname{IA}(A_\Lambda)$ and deduce a short exact sequence reminiscent of results from the Birman--Hilden theory for surfaces., Comment: 29 pages, 4 figures

We show that the automorphism groups of right-angled Artin groups whose defining graphs have at least 3 vertices are not relatively hyperbolic. We then show that the outer automorphism groups are not relatively hyperbolic, if they are not virtually isomorphic to a right-angled Artin group whose defining graph is either a single vertex or disconnected., Comment: We correct the theorem of relative hyperbolicity of outer automorphism groups of right-angled Artin groups. The corrigendum will appear in Bulletin of the Australian Mathematical Society

A data-driven framework is presented for building magneto-elastic machine-learning interatomic potentials (ML-IAPs) for large-scale spin-lattice dynamics simulations. The magneto-elastic ML-IAPs are constructed by coupling a collective atomic spin model with an ML-IAP. Together they represent a potential energy surface from which the mechanical forces on the atoms and the precession dynamics of the atomic spins are computed. Both the atomic spin model and the ML-IAP are parametrized on data from first-principles calculations. We demonstrate the efficacy of our data-driven framework across magneto-structural phase transitions by generating a magneto-elastic ML-IAP for {\alpha}-iron. The combined potential energy surface yields excellent agreement with first-principles magneto-elastic calculations and quantitative predictions of diverse materials properties including bulk modulus, magnetization, and specific heat across the ferromagnetic-paramagnetic phase transition., Comment: 15 pages, 6 figures, 2 tables

We present a methodology based on the N\'{e}el model to build a classical spin-lattice Hamiltonian for cubic crystals capable of describing magnetic properties induced by the spin-orbit coupling like magnetocrystalline anisotropy and anisotropic magnetostriction, as well as exchange magnetostriction. Taking advantage of the analytical solutions of the N\'{e}el model, we derive theoretical expressions for the parameterization of the exchange integrals and N\'{e}el dipole and quadrupole terms that link them to the magnetic properties of the material. This approach allows to build accurate spin-lattice models with the desire magnetoelastic properties. We also explore a possible way to model the volume dependence of magnetic moment based on the Landau energy. This new feature can allow to consider the effects of hydrostatic pressure on the saturation magnetization. We apply this method to develop a spin-lattice model for BCC Fe and FCC Ni, and we show that it accurately reproduces the experimental elastic tensor, magnetocrystalline anisotropy under pressure, anisotropic magnetostrictive coefficients, volume magnetostriction and saturation magnetization under pressure at zero-temperature. This work could constitute a step towards large-scale modeling of magnetoelastic phenomena.

Putman introduced a notion of a partitioned surface which is a surface with boundary with decoration restricting how the surface can be embedded into larger surfaces, and defined the Torelli group of the partitioned surfaces. In this paper, we study some topological and dynamical aspects of the Torelli groups of partitioned surfaces. More precisely, first we obtain upper and lower bounds on the cohomological dimension of Torelli groups of partitioned surfaces and show that those two bounds coincide when at most three boundary components are grouped together in the partition of the boundary. Second, we study the asymptotic translation lengths of Torelli groups of partitioned surfaces on the corresponding curve complexes. We show that the minimal asymptotic translation length asymptotically behaves almost like the reciprocal of the Euler characteristic of the surface. This generalizes the previous result of the first and second authors on Torelli groups for closed surfaces., Comment: 19 pages, 4 figures. A gap in the previous version has been fixed and a new result on cohomological dimension has been added

The magnetic behavior of bcc iron nanoclusters, with diameters between 2 and 8 nm, is investigated by means of spin dynamics (SD) simulations coupled to molecular dynamics (MD-SD), using a distance-dependent exchange interaction. Finite-size effects in the total magnetization as well as the influence of the free surface and the surface/core proportion of the nanoclusters are analyzed in detail for a wide temperature range, going beyond the cluster and bulk Curie temperatures. Comparison is made with experimental data and with theoretical models based on the mean-field Ising model adapted to small clusters, and taking into account the influence of low coordinated spins at free surfaces. Our results for the temperature dependence of the average magnetization per atom M(T), including the thermalization of the transnational lattice degrees of freedom, are in very good agreement with available experimental measurements on small Fe nanoclusters. In contrast, significant discrepancies with experiment are observed if the translational degrees of freedom are artificially frozen. The finite-size effects on M(T) are found to be particularly important near the cluster Curie temperature. Simulated magnetization above the Curie temperature scales with cluster size as predicted by models assuming short-range magnetic ordering (SRMO). Analytical approximations to the magnetization as a function of temperature and size are proposed.

We present a scale-bridging approach based on a multi-fidelity (MF) machine-learning (ML) framework leveraging Gaussian processes (GP) to fuse atomistic computational model predictions across multiple levels of fidelity. Through the posterior variance of the MFGP, our framework naturally enables uncertainty quantification, providing estimates of confidence in the predictions. We used Density Functional Theory as high-fidelity prediction, while a ML interatomic potential is used as the low-fidelity prediction. Practical materials design efficiency is demonstrated by reproducing the ternary composition dependence of a quantity of interest (bulk modulus) across the full aluminum-niobium-titanium ternary random alloy composition space. The MFGP is then coupled to a Bayesian optimization procedure and the computational efficiency of this approach is demonstrated by performing an on-the-fly search for the global optimum of bulk modulus in the ternary composition space. The framework presented in this manuscript is the first application of MFGP to atomistic materials simulations fusing predictions between Density Functional Theory and classical interatomic potential calculations.

Efficient algorithms for the calculation of minimum energy paths of magnetic transitions are implemented within the geodesic nudged elastic band (GNEB) approach. While an objective function is not available for GNEB and a traditional line search can, therefore, not be performed, the use of limited memory Broyden-Fletcher-Goldfarb-Shanno (LBFGS) and conjugate gradient algorithms in conjunction with orthogonal spin optimization (OSO) approach is shown to greatly outperform the previously used velocity projection and dissipative Landau-Lifschitz dynamics optimization methods. The implementation makes use of energy weighted springs for the distribution of the discretization points along the path and this is found to improve performance significantly. The various methods are applied to several test problems using a Heisenberg-type Hamiltonian, extended in some cases to include Dzyaloshinskii-Moriya and exchange interactions beyond nearest neighbours. Minimum energy paths are found for magnetization reversals in a nano-island, collapse of skyrmions in two-dimensional layers and annihilation of a chiral bobber near the surface of a three-dimensional magnet. The LBFGS-OSO method is found to outperform the dynamics based approaches by up to a factor of 8 in some cases.

We present a new algorithm to compute all the chiral polytopes that have a given group $G$ as full automorphism group. This algorithm uses a new set of generators that characterize the group, all of them except one being involutions. It permits to compute all chiral polytopes of groups that were previously unreachable by other known algorithms.

We show that a chiral coset geometry constructed from a $C^+$-group necessarily satisfies residual connectedness and is therefore a hypertope.

A temperature-dependent approach involving Green-Kubo equilibrium atomic and spin dynamics (GKEASD) is reported to assess phonon and magnon thermal transport processes accounting for phonon-magnon interactions. Using body-center cubic (BCC) iron as a case study, GKEASD successfully reproduces its characteristic temperature-dependent spiral and lattice thermal conductivities. The non-electronic thermal conductivity, i.e., the sum of phonon and magnon thermal conductivities, calculated using GKEASD for BCC Fe agrees well with experimental measurements. Spectral energy analysis reveals that high-frequency phonon-magnon scattering rates are one order of magnitude larger than those at low frequencies due to energy scattering conservation rules and high densities of states. Higher temperatures further accentuate this phenomenon. This new framework fills existing gaps in simulating thermal transport across the ferro- to para-magnetic transition. Future application of this methodology to phonon- and magnon-dominant insulators and semiconductors will enhance understanding of emerging thermoelectric, spin caloritronic and superconducting materials.

A parallel implementation of coupled spin-lattice dynamics in the LAMMPS molecular dynamics package is presented. The equations of motion for both spin only and coupled spin-lattice dynamics are first reviewed, including a detailed account of how magneto-mechanical potentials can be used to perform a proper coupling between spin and lattice degrees of freedom. A symplectic numerical integration algorithm is then presented which combines the Suzuki-Trotter decomposition for non-commuting variables and conserves the geometric properties of the equations of motion. The numerical accuracy of the serial implementation was assessed by verifying that it conserves the total energy and the norm of the total magnetization up to second order in the timestep size. Finally, a very general parallel algorithm is proposed that allows large spin-lattice systems to be efficiently simulated on large numbers of processors without degrading its mathematical accuracy. Its correctness as well as scaling efficiency were tested for realistic coupled spin-lattice systems, confirming that the new parallel algorithm is both accurate and efficient., Comment: 15 pages, 10 figures

When nano-magnets are coupled to random external sources, their magnetization becomes a random variable, whose properties are defined by an induced probability density, that can be reconstructed from its moments, using the Langevin equation, for mapping the noise to the dynamical degrees of freedom. When the spin dynamics is discretized in time, a general fluctuation-dissipation theorem, valid for non-Markovian noise, can be established, even when zero modes are present. We discuss the subtleties that arise, when Gilbert damping is present and the mapping between noise and spin degrees of freedom is non--linear., Comment: 3 pages LaTeX2e. Presentation for MMM2016, New Orleans

The Fokker--Planck equation describes the evolution of a probability distribution towards equilibrium--the flow parameter is the equilibration time. Assuming the distribution remains normalizable for all times, it is equivalent to an open hierarchy of equations for the moments. Ways of closing this hierarchy have been proposed; ways of explicitly solving the hierarchy equations have received much less attention. In this paper we show that much insight can be gained by mapping the Fokker--Planck equation to a Schr\"odinger equation, where Planck's constant is identified with the diffusion coefficient., Comment: 7 pages, LaTeX2e, 2 EPS figures. Uses jpconf style file. Contribution to the ICM-SQUARE conference, 23-26 May 2016, Athens, Greece. Accepted for publication in the proceedings

A functional calculus approach is applied to the derivation of evolution equations for the moments of the magnetization dynamics of systems subject to stochastic fields. It allows us to derive a general framework for obtaining the master equation for the stochastic magnetization dynamics, that is applied to both, Markovian and non-Markovian dynamics. The formalism is applied for studying different kinds of interactions, that are of practical relevance and hierarchies of evolution equations for the moments of the distribution of the magnetization are obtained. In each case, assumptions are spelled out, in order to close the hierarchies. These closure assumptions are tested by extensive numerical studies, that probe the validity of Gaussian or non--Gaussian closure Ans\"atze., Comment: 17 pages, 5 figures

Statistical averaging theorems allow us to derive a set of equations for the averaged magnetization dynamics in the presence of colored (non-Markovian) noise. The non-Markovian character of the noise is described by a finite auto-correlation time, tau, that can be identified with the finite response time of the thermal bath to the system of interest. Hitherto, this model was only tested for the case of weakly correlated noise (when tau is equivalent or smaller than the integration timestep). In order to probe its validity for a broader range of auto-correlation times, a non-Markovian integration model, based on the stochastic Landau-Lifshitz-Gilbert equation is presented. Comparisons between the two models are discussed, and these provide evidence that both formalisms remain equivalent, even for strongly correlated noise (i.e. tau much larger than the integration timestep)., Comment: 4 pages LaTeX2e, 3 EPS figures; uses IEEEtran.cls

A stochastic approach for the description of the time evolution of the magnetization of nanomagnets is proposed, that interpolates between the Landau-Lifshitz-Gilbert and the Landau-Lifshitz-Bloch approximations, by varying the strength of the noise. Its finite autocorrelation time, i.e. when it may be described as colored, rather than white, is, also, taken into account and the consequences, on the scale of the response of the magnetization are investigated. It is shown that the hierarchy for the moments of the magnetization can be closed, by introducing a suitable truncation scheme, whose validity is tested by direct numerical solution of the moment equations and compared to the averages obtained from a numerical solution of the corresponding colored stochastic Langevin equation. This comparison is performed on magnetic systems subject to both an external uniform magnetic field and an internal one-site uniaxial anisotropy., Comment: 4 pages, 3 figures

Frequency-dependent permeability tensor for unsaturated polycrystalline ferrites is derived through an effective medium approximation that combines both domain-wall motion and rotation of domains in a single consistent scattering framework. Thus derived permeability tensor is averaged on a distribution function of the free energy that encodes paramagnetic states for anhysteretic loops. The initial permeability is computed and frequency spectra are given by varying macroscopic remanent field., Comment: 24 pages, 3 figures

We propose a stochastic approach for the description of the time evolution of the magnetization of nanomagnets, that interpolates between the Landau--Lifshitz--Gilbert and the Landau--Lifshitz--Bloch approximations, by varying the strength of the noise. In addition, we take into account the autocorrelation time of the noise and explore the consequences, when it is finite, on the scale of the response of the magnetization, i.e. when it may be described as colored, rather than white, noise and non-Markovian features become relevant. We close the hierarchy for the moments of the magnetization, by introducing a suitable truncation scheme, whose validity is tested by direct numerical solution of the moment equations and compared to the average deduced from a numerical solution of the corresponding stochastic Langevin equation. In this way we establish a general framework, that allows both coarse-graining simulations and faster calculations beyond the truncation approximation used here., Comment: 5 pages LaTeX2e, 2 EPS figures; uses elsarticle.cls. Presented at HMM 2015, 10th International Symposium on Hysteresis Modeling and Micromagnetics

The Landau--Lifshitz--Gilbert equations for the evolution of the magnetization, in presence of an external torque, can be cast in the form of the Lorenz equations and, thus, can describe chaotic fluctuations. To study quantum effects, we describe the magnetization by matrices, that take values in a Lie algebra. The finite dimensionality of the representation encodes the quantum fluctuations, while the non-linear nature of the equations can describe chaotic fluctuations. We identify a criterion, for the appearance of such non-linear terms. This depends on whether an invariant, symmetric tensor of the algebra can vanish or not. This proposal is studied in detail for the fundamental representation of $\mathfrak{u}(2)=\mathfrak{u}(1)\times\mathfrak{su}(2)$. We find a knotted structure for the attractor, a bimodal distribution for the largest Lyapunov exponent and that the dynamics takes place within the Cartan subalgebra, that does not contain only the identity matrix, thereby can describe the quantum fluctuations., Comment: 5 pages, 3 figures. Uses jpconf style. Presented at the ICM-SQUARE 4 conference, Madrid, August 2014. The topic is a special case of the content of 1404.7774, currently under revision

A generalization of the Lorenz equations is proposed where the variables take values in a Lie algebra. The finite dimensionality of the representation encodes the quantum fluctuations, while the non-linear nature of the equations can describe chaotic fluctuations. We identify a criterion, for the appearance of such non-linear terms. This depends on whether an invariant, symmetric tensor of the algebra can vanish or not. This proposal is studied in detail for the fundamental representation of $\mathfrak{u}(2)$. We find a knotted structure for the attractor, a bimodal distribution for the largest Lyapunov exponent and that the dynamics takes place within the Cartan subalgebra, that does not contain only the identity matrix, thereby can describe the quantum fluctuations., Comment: 10 pages Revtex, 3 figures