Significance Lossless transport is the Holy Grail of energy science in general and superconductivity research in particular. The main obstacle is the dissipative motion of Abrikosov vortices, which can be reduced or eliminated by pinning at nonsuperconducting defects. Pinning effectiveness nontrivially depends on various factors such as the shape, concentration, and spatial distribution of defects, rendering the optimization of the vortex pinning landscape highly difficult. Here, we use concepts from biological evolution to develop an efficient strategy for vortex pinning improvement. We replace natural selection with targeted selection, where only pinning configurations with better vortex immobilization survive. In combination with high-performance numerical algorithms, it allows us to dynamically evolve the defect landscape into the best possible pinning configuration with maximal lossless current.

Abstract The ability of type II superconductors to carry large amounts of current at high magnetic fields is a key requirement for future design innovations in high-field magnets for accelerators and compact fusion reactors, and largely depends on the vortex pinning landscape comprised of material defects. The complex interaction of vortices with defects that can be grown chemically, e.g., self-assembled nanoparticles and nanorods, or introduced by postsynthesis particle irradiation precludes a priori prediction of the critical current and can result in highly nontrivial effects on the critical current. Here, we borrow concepts from biological evolution to create a vortex pinning genome based on a genetic algorithm, naturally evolving the pinning landscape to accommodate vortex pinning and determine the best possible configuration of inclusions for two different scenarios: a natural evolution process initiating from a pristine system and one starting with preexisting defects to demonstrate the potential for a postprocessing approach to enhance critical currents. Furthermore, the presented approach is even more general and can be adapted to address various other targeted material optimization problems.

Life has undergone tremendous changes due to natural selection—from relatively simple molecules with replication capability to complex organisms, whose understanding is still far beyond present contemplation. Modern computer systems have enabled the effective exploitation of the idea of natural selection for practical purposes. The underlying genetic algorithms are widely used in electromagnetic and mechanical design, financial mathematics, energy applications, scheduling problems, circuit design, image processing, medicine, etc. Within this approach, one only needs to specify the direction of positive mutations to find optimal or beneficial characteristics of the system of interest, i.e., replace natural evolution by targeted evolution, which is especially effective in complex systems with a large number of degrees of freedom.

A key science aspect to advance the deployment of high-temperature superconductors (HTSs) is the discovery of novel materials that can carry large currents without dissipation at high-magnetic fields (1). These materials are especially desirable for high-performance applications (2) such as superconducting motors, generators, magnets, and power lines in urban areas (3⇓–5). Low dissipation is also very important for superconducting cavities for particle accelerators (6), undulators for x-ray synchrotrons (7), and compact fusion reactors (8). The main challenge is to suppress the dissipation in these systems caused by the motion of quantized elastic magnetic flux tubes or vortices, which appear in type II superconductors in magnetic fields above the first critical field (9). Since most applied superconductors are of type II, the study of efficient pinning mechanisms benefits a majority of superconducting technologies. Vortices can be pinned (or trapped) by inhomogeneities in the material, usually in the form of nonsuperconducting defects (10). Examples are point-like pinning centers (impurities, vacancies, inclusions), 1D defects (dislocations, irradiation tracks), or 2D defects (twin boundaries, stacking faults). Although extensive knowledge has been gained in the pursuit of high critical currents (the highest current the system can carry without dissipation) (11⇓⇓⇓–15), the fundamental solution to the dynamics of interacting vortices in disordered media is still unknown. Only recently more systematic, computer-assisted approaches were developed (16, 17), leading to the critical-current-by-design methodology (18).

While sophisticated numerical optimization methods (19) and corresponding experiments can guide the design of superconductors with enhanced critical current densities ( J c ), the problem requires defining the general geometry of the vortex pinning landscape (or pinscape) a priori. This works well if only a certain type of pinning defects is present—in other words, pinscapes defined by only a few parameters. Hence, the overall best pinscape for the highest J c cannot be determined by these approaches. To address this question, one needs to study all possible combinations of defects, resulting in highly mixed pinscapes. Each of the individual defects are described by numerous material and geometrical parameters, resulting in an extremely high-dimensional parameter space for the pinscape. This is where evolutionary concepts can be used.

In this work, we borrow concepts from biological evolution to create a vortex pinning genome with targeted evolution for predicting high in-field J c . We focus on the geometrical aspect of the defects to produce the best pinscape for a given system. In particular, we evolve the pinscape by changing the shapes of individual defects (see sketch of targeted evolution in Fig. 1), thereby including the possibility of all major defect types such as columnar and spherical defects, which can be experimentally realized. Moreover, our approach can also be adapted for many different materials optimization/design problems. Here, we demonstrate its power for (i) numerically determining the maximum possible J c in superconductors with nonmagnetic normal inclusions and (ii) developing a universal postprocessing strategy for enhancing the performance of superconductors with preexisting pinning landscapes such as in commercial HTS wires and superconductors in alternating or nonhomogeneous magnetic fields.

Fig. 1. Sketch of a targeted evolution of the pinning landscape. We start with generation 0, which contains a single configuration without defects. Each defect has elliptical shape and is characterized by three independent diameters. The evolution process mutates the pinning landscape by adding/removing, translating, scaling, and reshaping particles. These mutations create the next generation. We accept the pinning landscape with maximal critical current density ( J c ) and discard all others. The evolution ends at some generation N with configuration having maximal J c (shown in red).

Targeted Evolution An essential ingredient for our approach is to obtain the critical current for a given evolved pinscape. Here, we describe the collective dynamics and pinning of vortices by the time-dependent Ginzburg–Landau equation (TDGLE), allowing us to determine J c (16, 17)—the fitness function. The TDGLE yields J c as function of shapes, sizes, and positions of pinning defects (see details in SI Appendix). As a quite general model, we consider pinning landscapes containing D ellipsoidal metallic pinning centers with principle axes ( a i , b i , c i ), aligned in the x, y, and z directions, with center positions ( x i , y i , z i ), where i = 1 , … , D . These ellipsoidal defects can describe a large variety of defect geometries in superconductors such as precipitates, point defects, dislocations, grain boundaries, and stacking faults, as well as particle irradiation-induced columnar or spherical defects. For example, point defects can be modeled by small spherical inclusions, grain boundaries by flattened spheroids and columnar defects by spheroids with one of the diameters larger than the system size. To find pinscapes with ellipsoidal defects that yield the highest J c , we use an evolution-based algorithm with three distinct stages: mutations and targeted selection (stage 1), extrapolation and analysis (stage 2), and verification (stage 3), described below. Stage 1: Mutations and Targeted Selection. This step implements the evolutionary paradigm, during which the shape and position of individual inclusions is altered independently (mutation) and J c is calculated. A set of random mutations produces a new generation. Each new pinscape or successor may contain one (typical) or more sequential mutations (rare). Each pinscape is evaluated, and the one with the largest J c is chosen for further evolution (see sketch in Fig. 1). The initial pinscape usually depends on the problem to be studied. For a general situation (discussed in the next section), one can initiate the targeted evolution algorithm with an empty pinscape, the 0th generation, which represents an infinite homogeneous system with zero J c . Mutations have random type, strength, direction, and number of affected inclusions, namely: (i) copying of existing inclusions or adding new inclusions of random shape; (ii) removing inclusions; (iii) changing the inclusion principle axes a i , b i , and c i ; (iv) changing the inclusion position ( x i , y i , z i ); (v) repelling/attracting pairs of inclusions, i.e., increasing/decreasing the distance between randomly chosen inclusions i and j; (vi) squishing inclusions, i.e., changing the inclusion’s axes a i , b i , and c i while maintaining its volume; (vii) splitting inclusions, i.e., creating a pair of inclusions with the same volume as the original one; and (viii) merging pairs of inclusions. Mutation types (vi) to (viii) preserve the volume of the affected inclusions. Note, that if we start the mutation process with an empty pinscape, the only possible mutation is the addition of defects. We calculate J c for each pinscape in a generation. These are then compared with the maximal J c of the previous generation. In case none of the mutations increase J c , we repeat the mutation procedure and expand the population in the current generation until at least one pinscape produces a J c larger than the maximal critical current of the previous generation or a maximum population is reached. This stage is implemented to work in parallel. If a configuration with larger J c is found within a generation, we select the pinscape with the largest J c as the seed configuration for the following generation and then apply the mutation procedure again. Repeating this protocol produces subsequent generations of pinscapes with even higher J c . We stop in generation N if no further J c enhancement is found (the cutoff population size is 2,048). The evolution approach provides us with the types and parameters of defects that ensure maximum vortex pinning and, consequently, maximum J c . The results are obtained without any assumptions of the pinscape structure and only depend on external parameters such as magnetic field and temperature. In some application relevant situations, the initial pinscape and the type of possible mutations may have some constraints, which we address below. Stage 2: Extrapolation and Analysis. Stage 1 provides information regarding the distribution of the particle sizes and, in some cases, their spatial distribution. We can model/extrapolate these distributions with only a few parameters such as the size and typical distances between defects. In other words, one can use the general knowledge of the defect shapes obtained by the evolutionary approach and characterize the corresponding pinscape with a simplified global parameter set. For example, if the optimal pinscape consists of randomly distributed spherical defects of similar diameters, the configuration can be characterized by two parameters: concentration and diameter of the defects (16). Based on the simplified global parameter set, near-optimal pinscapes can be fine-tuned using conventional optimization methods (19). Furthermore, one can sample J c for near-optimal parameter sets to determine the robustness of the configuration and compare them to analytical results (10, 20). Stage 3: Verification. To test the model obtained in stage 2, we restart the evolution process with the best model configuration and change the positions and sizes of each inclusion individually. The model is verified, if subsequent evolution cannot further increase J c by a significant amount (we typically use a threshold of 3% within 2,048 mutations). Stages 2 and 3 are in a sense optional, as they elucidate the underlying mechanism for the optimal pinscape, extract a model, and show the stability of the process. Stage 1 alone can determine the general optimal pinscapes. Pinning Landscape for Maximal Critical Current Density in Fixed Field. Starting with empty pinscapes and allowing almost any possible mutation is usually not very relevant for applications. However, it is instructive to study this case as it ultimately yields the best pinning configurations for given external parameters. Consider the exemplary situation of a fixed magnetic field applied along the z axis (or c axis in HTSs) and current flowing along the x direction. Naïvely, the optimal pinning landscape should mimic the vortex configuration for zero applied current at the given field, namely the Abrikosov vortex lattice. Hence, the pinscape should be a hexagonal array of columnar defects, with each column trapping a single vortex. However, the evolutionary approach yields an even better pinscape: a periodic array of planar pinning defects (walls) that are aligned with the current and field direction (here parallel to the x z plane). In the simulation, we apply a constant external magnetic field B = 0.1 H c 2 at low temperatures, corresponding to nearly zero noise (reduced temperature T f = 1 0 − 5 ; see ref. 21 for details). Inclusions are modeled by a nonsuperconducting material with zero critical temperature, T c , i , resulting in a suppressed order parameter, ψ ( r ) , inside the defects. Here, H c 2 is the upper critical field at given temperature. Following the evolution approach described above, an actual evolution process for of an initially empty pinscape is illustrated in Fig. 2. Note that J c rises faster in early generations; improvements in later generations require more mutations and lead to a smaller gain in J c . The evolution terminates with the 37th generation and results in a set of almost equidistant planar defects oriented in the direction of applied current and having a thickness on the order of a coherence length (Fig. 2B). The distance between planar defects roughly corresponds to positions of vortex rows in a perfect hexagonal lattice (blue circles) generated by the external magnetic field. The full evolution tree has 37 generations and 6,331 pinning configurations (Fig. 2C). The best landscapes in each generation are numbered and have color ranges from blue with almost zero J c to orange with maximal J c = 0.40 J d p , where J d p is the depairing current. Each numbered configuration has at least 20 successors: (i) the successor with maximum J c becomes the numbered seed for the next generation; (ii) all other successors are shown by small colored circles. These configurations have smaller J c values than the seed and are discarded by our targeted selection. The final 37th configuration in the center has 2,040 mutations with smaller J c values. The vast majority (90%) of these dead mutations lead to marginal decreases of J c (from 1 to 15%, shown in green). Therefore, the determined configuration, shown in Fig. 2B, is rather stable with respect to mutation. Fig. 2. Evolution history. (A) The evolution process starts with a superconductor without inclusions shown in the left panel. The following panels show pinning landscapes having highest J c in first, third, fourth, sixth, ninth, 12th, and 21st generation, correspondingly. In the first generation, the maximal J c is achieved with the configuration containing a single nearly spherical inclusion. In second and third generations, this inclusion evolves to a flattened ellipsoid lying in the plane spanned by the current and magnetic field. In subsequent generations, this ellipsoid is copied multiple times to enhance the total pinning. The remaining generations of the evolution process fine tunes the landscape by copying, removing, moving, and slightly deforming successors of the seed inclusion. (B) The final pinning landscape consisting of a periodic array of almost planar defects has the best possible J c in the framework of our model. The positions of pinned vortices are shown schematically by blue circles. (C) The evolution tree. The numbered circles represent configurations with the maximum J c per generation. Dead mutations are indicated by dots. All J c values are color-coded. The used parameters produce a rather large critical current, J c ( T ) = 0.40 J d p ( T ) , at almost zero temperature. For a larger noise level ( T f = 0.28 ) corresponding to a temperature T ∼ 77 K, J c ( T ) reduces to 0.34 J d p ( T ) . Weaker metallic pinning centers with higher critical temperature, e.g., T c , i = 2 T − T c , b ( T c , b is the critical temperature in the bulk superconductor) can produce a maximal J c = 0.31 J d p at zero noise. In all these cases, we can easily extrapolate a model for the optimal pinning configuration with only two parameters: the thickness of the planar defects and their separation. Based on the optimal pinscape, we studied the critical dynamics close to J c (see SI Appendix, Fig. S1 and Movies S1–S3 for corresponding order parameter and supercurrent densities). Namely, the depinning process involves the collective motion of vortices, which effectively increases the pinning force of the system. The same collective behavior occurs for other types of J c -optimized pinning landscapes, e.g., for ordered defects (22) or disordered nanorods extended along the direction of the applied magnetic field (18). A similar but somewhat less pronounced effect was observed for randomly placed spherical particles (16). We conclude that the collective depinning in very-large- J c systems leads to a much more pronounced and sharp transition to the dissipative state than individual vortex depinning in suboptimal pinscapes, confirmed by current–voltage curves (SI Appendix, Fig. S2 and Movies S3 and S4). Our targeted evolution derived pinning landscape with J c = 0.40 J d p can be compared with other typical pinscapes with potentially high J c at the same magnetic field B = 0.1 H c 2 and low thermal noise: (i) randomly placed spherical defects with optimal diameter and concentration have a maximum possible J c of 0.061 J d p (16); (ii) field-aligned, randomly placed columnar inclusions with the best diameter and concentration lead to J c = 0.091 J d p (19); (iii) hexagonally ordered, field-aligned columnar defects with optimal size and concentration generate a significantly larger critical current J c = 0.32 J d p but still smaller than for planar defects. Next, we compare the properties of hexagonally ordered columnar with that of arrays of planar defects—the idealized model derived from the genetic approach. Planar vs. Columnar Defects. A systematic comparison of a hexagonal lattice of columnar defects to the extrapolated model of a periodic array of planar defects requires comparable parameters. The natural parameters for columnar defects are the matching field B Φ (the hypothetic magnetic field producing an Abrikosov vortex lattice with the same density as the lattice of columns) and their diameter d. For arrays of planar defects, one can use the same matching field B Φ and place the defects along one of the main axes (parallel to the current) of the hexagonal lattice (and along the field) (Fig. 3 A and B, Insets). The distance between the planar defects is then h = 3 1 / 4 ξ ( π H c 2 / B Φ ) 1 / 2 , i.e., 3 1 / 2 / 2 less than for columnar defects. The second parameter is the thickness of the planar defects. In both cases, the maximum J c is reached when B = B Φ (Fig. 3 C and D). A main difference is that the planar array is more robust against changes in B Φ than the discrete columnar defects structure, i.e., small changes in B Φ (or h) result in very small changes in the optimal J c . Fig. 3. Critical current as a function of magnetic field, J c ( B ) , for different landscapes and matching fields B Φ . (A) Planar defects with fixed thickness of b = 0.5 ξ and range of B Φ from 0.025 to 0.4 H c 2 . The star shows the targeted evolution result for B = 0.1 H c 2 . The envelope curve (black line with open circles) shows the maximal possible J c , m a x ( B ) at a given field B. Inset shows the corresponding optimal matching field B Φ , m a x , i.e., the distance needed to achieve this maximum. (B) Hexagonal pattern of columnar defects with diameter d = 3 ξ . Inset shows J c as a function of hexagonal lattice rotation angle α with respect to the applied current. It is π/3-periodic and maximal if the current is aligned with the lattice axes. (C) J c as a function of matching field for planar defects in applied field B = 0.1 H c 2 for different wall thicknesses, b. (D) The same for columnar defects ordered in a hexagonal pattern for different diameters, d, of the columns. Fig. 3A shows the J c ( B ) dependence for the planar array with fixed thickness, b = 0.5 ξ , and different B Φ ranging from 0.025 H c 2 to 0.4 H c 2 . All curves display a relatively smooth behavior. The most representative is the green curve simulated for B Φ = 0.1 H c 2 ; the open star on the green curve corresponds to J c associated with the pinning landscape shown in Fig. 2B obtained by targeted evolution for B = 0.1 H c 2 . For the green curve J c ( B ′ ) ⩾ J c ( 0.1 H c 2 ) for all B ′ ⩽ 0.1 H c 2 ; this property persists for all reasonable wall-pattern parameters with optimized B Φ and b at a given field B, i.e. J c ( B ′ ) ⩾ J c ( B ) for B ′ ⩽ B . The envelope curve (black line with circles) shows J c , m a x ( B ) , the critical current for optimized landscapes for each fixed B with optimal wall thickness b ( B ) and matching field B Φ , m a x ( B ) . The deviation of B Φ , m a x ( B ) from a simple linear dependence B Φ = B (Inset) is due to different b ( B ) , ranging from ∼ 0.5 ξ at low fields to ∼ 0.1 ξ at higher fields. Fig. 3B shows the J c ( B ) dependence for hexagonal-patterned columnar defects with fixed diameter d = 3 ξ for different B Φ from 0.025 H c 2 to 0.2 H c 2 . The green curve shows a peak at the first matching field with J c = 0.32 J d p , which coincides with the maximal J c of the hexagonal lattice at B = 0.1 H c 2 . A rotation of the hexagonal pattern can reduce this value (it is maximum if a main axis of the lattice is aligned with the current; see angular dependence in Inset). Fig. 3C depicts the J c ( B Φ ) dependence for arrays of planar defects with different wall thickness b at fixed applied field B = 0.1 H c 2 . This sampling shows a single robust optimum near B Φ = 0.1 H c 2 and b = 0.5 ξ . A similar sampling for the hexagonal columnar defect pattern presented in Fig. 3D shows significantly sharper peaks in the vicinity of the matching field, resulting in less robust behavior against small changes of the parameters. Samplings for other parameters are shown in SI Appendix, Figs. S4–S8.

Application-Relevant Examples of Targeted Evolution A recent report on doubling J c of commercial HTS wire by additional particle irradiation (23) highlights the importance and advantages of a postsynthesis approach to enhance J c , while leaving the wire synthesis process untouched. Our targeted evolution approach can also be applied to systems with preexisting defects. Fig. 4 demonstrates results of targeted evolution in different environments, defined by either preexisting pinscapes or different external parameters. In Fig. 4 A and B, we apply the evolutionary algorithm to pinscape with fixed preexisting defects. These defects partially block the left-to-right current flow and, thus, dramatically change the result of the targeted evolution described above. Mainly, the evolution leaves some defect-free regions in the superconducting matrix to allow for a supercurrent path. In the case of preexisting tilted walls in Fig. 4A, the total current I c = J c w t through the system was increased by evolution from I c = 56 J d p ξ 2 ( J c = 0.11 J d p ) to I c = 147 J d p ξ 2 ( J c = 0.29 J d p ) in applied field B = 0.1 H c 2 , where w = 64 ξ and t = 8 ξ are the system’s width and thickness, respectively. In the case of the preexisting two half-ellipses shown in Fig. 4B, the critical current rises, from I c = 35 J d p ξ 2 ( J c = 0.068 J d p ) to I c = 104 J d p ξ 2 ( J c = 0.20 J d p ) upon evolution of added defects. Fig. 4. Pinscape evolution for predefined environments. (A) Current is applied from left to right, and magnetic field is fixed at B = 0.1 H c 2 perpendicular to the figure plane as in Fig. 2. The difference is in the preexisting pinscapes containing tilted planar defects shown in gray. These plates redirect the supercurrent flow (boundary conditions are periodic in the figure plane) and make the optimal pinscape shown in Fig. 2 inefficient. The evolutionary approach generates smaller planar defects (or flat cylinders) along the current between the preexisting inclusions. (B) In this scenario, two flattened half-cylinders block the current with open (no-current) boundary conditions at the top and bottom surfaces. Generated inclusions are more cylindrical between inclusions to avoid blocking supercurrents. (C) In this scenario, the current can be applied in left–right and bottom–up directions, and the largest critical current is defined by the minimum critical current in either direction. The pinscape evolves to hyperuniformly placed columnar defects. In Fig. 4C, we apply the current both in the horizontal and vertical directions and consider the fitness function J c , u = min { J c → , J c ↑ } , where J c → is left-to-right J c and J c ↑ is bottom-to-up J c , rather than only J c → as before. J c , u approximately models arbitrary directions of applied currents. The resulting pinscape consists of columnar defects along the magnetic field arranged in a hyperuniform pattern (22, 24). The corresponding critical current density, J c = 0.27 J d p , is 5% less than the J c for a hexagonal lattice oriented in the wrong way (rotated π / 6 from the main axes; see the angular dependence in Fig. 3B, Inset). In all of the simulations above, we intentionally did not limit the size, shape, or placement of the mutated defects. However, it is possible to limit the defect morphology to mimic the limitations of practical postprocessing procedures.

Discussion and Conclusions In this paper, we introduced an evolutionary approach for the optimization of pinscapes in type II superconductors. This approach utilizes the idea of targeted selection inspired by biological natural selection. We demonstrated that it can be applied to enhance the current-carrying capacity of superconductors in a magnetic field. We discovered that certain patterns of defects composed of metallic inclusions can maximize the critical current up to 40% of J d p for fixed direction of the current perpendicular to the magnetic field at 10% of H c 2 . We numerically demonstrated that no other mixture of different defect shapes can reach this level of J c . The discussed pining structure may arise in niobium titanium wires, in which a sequence of heating/drawing steps result in a microstructure composed of nanometer-scale metallic and almost parallel α-titanium lamellae embedded in the niobium titanium matrix (25). Furthermore, the layered structure of cuprate HTSs give rise to intrinsic pinning of similar nature. In contrast to conventional optimization techniques such as coordinate descent, where one varies only a few parameters characterizing the entire sample (e.g., size and concentration of defects), our targeted evolution approach allows us to vary each defect individually without any a priori assumptions about the defects configuration. This flexibility outweighs its higher computational cost. The considered optimization problem has basically infinite degrees of freedom, prompting one to ask why the evolution method convergences relatively quickly. One reason is that there are a lot of configurations with J c quite close to the maximum possible one, which are in practice, indistinguishable from each other. The evolutionary approach just allows us to find one such configuration. Typically, larger regions of near-optimum configurations correspond to a broader maximum of J c as a function of a set of appropriate parameters, e.g., the system in Fig. 3C evolutionally adapts faster than the system in Fig. 3D. We also demonstrated the enhancement of J c for two cases of preexisting defects, found in commercial HTSs. Our approach provides a computer-assisted route to rational enhancement of the critical current in applied superconductors. It can be used to define a postsynthesis optimization step for existing state-of-the-art HTS wires for high-field magnet applications by modeling the actual geometry of the wire within the magnet and taking into account external magnetic field distributions and self-fields. This can be done by coupling transport simulations with Maxwell equations and initiating the simulation with a preexisting defect distribution in the wire. Finally, we note that the described evolutionary algorithm is a local method and thus can easily get stuck in a local maximum. An analog in biological evolution is the extreme detour of a giraffe’s recurrent laryngeal nerves (26), which became trapped under the aortic arch in the thorax. However, in contrast to natural selection, targeted evolution can be performed multiple times. Namely, a comparison of the resultant pinscapes and corresponding J c values allows us to estimate how close they are to the best possible pinscape, making targeted evolution global. Moreover, by finding different near-maximum points, it is possible to understand which parameters are important for large J c and which ones are not. An experimental analog in organic systems is the process of in vitro selection (27). A particular example is the selection of RNA molecules being able to bind to specific ligands (28): it was shown that evolved molecules bind stronger than those of the first generation and an a priori guess of the best binding RNA sequence would not have been possible. In conclusion, our methodology of using targeted evolutionary concepts to improve the intrinsic properties of condensed matter systems is a promising path toward the design of tailored functional materials. It can be applied to a large variety of different physical systems and has demonstrated its usefulness in the enhancement of superconducting critical currents. Furthermore, its ability to take existing environments into account allows for optimization by postprocessing.

Materials and Methods The evolutionary algorithm was implemented in Python, and the TDGLE simulations were implemented for high-performance computers with general-purpose graphics processing unit coprocessors; see details and used parameters in SI Appendix.

Acknowledgments We thank L. Civale, R. Willa, and I. S. Aranson for numerous useful comments. I.A.S., A.E.K., and A.G. were supported by the Scientific Discovery through Advanced Computing program, funded by US Department of Energy (DOE), Office of Science, Advanced Scientific Computing Research and Basic Energy Science, Division of Materials Science and Engineering. U.W. and W.-K.K. were supported by the Center for Emergent Superconductivity, an Energy Frontier Research Center funded by the US DOE, Office of Basic Energy Sciences. Simulations were performed at the Oak Ridge Leadership Computing Facility (LCF) (DOE Contract DE-AC05-00OR22725), the Argonne LCF (DOE Contract DE-AC02-06CH11357), and the Computing Facility at Northern Illinois University.

Footnotes Author contributions: I.A.S., W.-K.K., U.W., and A.G. designed research; I.A.S., A.E.K., and A.G. performed research; I.A.S. and A.G. contributed new reagents/analytic tools; I.A.S. and A.G. analyzed data; and I.A.S., A.E.K., W.-K.K., U.W., and A.G. wrote the paper.

The authors declare no conflict of interest.

This article is a PNAS Direct Submission.

See Commentary on page 10201

This article contains supporting information online at www.pnas.org/lookup/suppl/doi:10.1073/pnas.1817417116/-/DCSupplemental.