**Physicists Thriving with Paperless Publishing**

Heath B. O'Connell[^1]
Stanford Linear Accelerator Center
Stanford CA 94309
`firstname.lastname@example.com`
`http://www.slac.stanford.edu/grp/lib/people/hoconnell.html`

The Stanford Linear Accelerator Center (SLAC) and Deutsches Elektronen Synchrotron (DESY) libraries have been comprehensively cataloguing the High Energy Particle Physics (HEP) literature online since 1974. The core database, SPIRES-HEP, now indexes over 400,000 research articles, with almost 50% linked to fulltext electronic versions (this site now has over 15 000 search hits per day). This database motivated the creation of the first site in the United States for the World–Wide Web at SLAC. With this database and the invention of the Los Alamos E-print archives in 1991, the HEP community pioneered the trend to "paperless publishing" and the trend to paperless access; in other words, the "virtual library." We examine the impact this has had both on the way scientists research and on paper-based publishing. The standard of work archived at Los Alamos is very high. 70% of papers are eventually published in journals and another 20% are in conference proceedings. As a service to authors, the SPIRES-HEP collaboration has been ensuring that as much information as possible is included with each bibliographic entry for a paper. Such meta-data can include tables of the experimental data that researchers can easily use to perform their own analyses as well as detailed descriptions of the experiment, citation tracking, and links to full-text documents.

# Introduction

The World Wide Web is now ubiquitous. At a time when every advertisement on television seems to end with a URL, it is of some interest to note that the SLAC Library had the first Web server in North America (indeed the first outside Europe), which it set up to aid the research of High Energy/Particle physicists. In this talk we explore the role technology has played in the organisation and dissemination of information in High Energy Physics (HEP), where paper is now a medium of last resort, and discuss why HEP provided the ideal conditions for the rapid adoption of new technology.

# The Compulsive Communicators

Communication has long played a vital role in physicists' research. Not content to wait for the publication of their work in journals, they, or perhaps their departments, would send out *preprints* of their work to institutions or other researchers they personally knew (once the duplication process became affordable). Naturally because of the cost and effort involved in this, large departments and famous scientists enjoyed a numerical advantage both in sending and receiving. However, for every poor soul in a small, out of the way department lamenting intellectual starvation, a well known researcher would be suffering from an overstuffed mailbox. By the 1960's the sheer number of preprints had somewhat perversely made the communication of research *more* difficult, due to information overload.

Since its foundation in 1962, the SLAC Library, at the request of the director, W.K.H. Panofsky, had been actively collecting new preprints (as a world centre for physics, it attracted some three thousand articles per year). The details of each preprint were entered into a card catalogue, including every author, no matter how many there were on the paper, which for an experiment could easily be in the hundreds (the current record is over 1000). From this catalogue, the library published a weekly list of new preprints for the SLAC community.

In 1967 work began at Stanford University on a computer database system that would be able to handle a practically limitless number of large bibliographic records. It was decided early on that the SLAC Library's preprint card catalogue would be a suitable test subject for this database, and so in March 1968 the SLAC Library began participating in the development and testing of SPIRES, the Stanford Physics Information REtrieval System. SPIRES (later renamed the Stanford Public Information REtrieval System) was envisioned as a fully searchable database, and each bibliographic record would contain the preprint's author list (with each author's institution), title, report number, date.

The remaining bibliographic element, which turned out to be of vital importance, was the paper's reference list, that is the list of earlier papers the present paper has cited . Only references to published journal articles would be included (as these were the most likely to be accurate), and they would be stored by the journal, volume and first page number. The importance of these reference lists was with SPIRES, one could search through them for a given paper, and thus find every preprint that *cited* the paper, thereby establishing (in some way) the impact this paper had in the literature. Naturally, this is of considerable interest to physicists .

However, this activity was not limited to the New World. In 1962, the Deutsches Elektronen-Synchrotron (DESY) in Hamburg, Germany, had begun to publish a list named "High Energy Physics - An Index," (HEPI) of all the research literature it received. Included in this list was both the published literature and unpublished preprints (whereas SLAC concentrated mainly on preprints), amounting to at that time some nine thousand documents per year. Like SLAC, DESY provided the title, author and date, but instead of references, HEPI staff physicists assigned to each paper up to 23 keywords (or, more often, key phrases) chosen from a controlled list.

Naturally it made sense for the SLAC and DESY teams to work together in their coverage of the particle physics literature, and so preparations for this began. By June 1969 specifications for the conversion of the DESY data to SPIRES format were complete.

# SPIRES and the Internet 1969–1990

In 1968 the SLAC librarians, Louise Addis, Bob Gex and Rita Taylor were funded by the the US Atomic Energy Commission (precursor to today's Department of Energy), with the Division of Particles and Fields (DPF) of the American Physical Society as a sponsor, to publish a weekly preprint list for the entire particle physics community . This list of all preprints entered in the past week could be easily generated from the newly commissioned SPIRES-HEP database. In January 1969 the first Preprints in Particles and Fields (PPF) list was sent out to over a thousand eager subscribers. Therefore, to get your work known all around the world (or at least the title of your work seen) all you would have to do would be send your preprint to SLAC, and similarly you could browse the titles of all the week's preprints by just looking at one list. This established the idea of a central repository for preprints.

A sister "Anti-preprints" list (PPA) recorded the preprints that had since been published in journals. This allowed the original preprints to be discarded, and proved quite popular with the journal editors, who were then able to match references to preprints with the published article.

By 1974 both the SLAC and DESY libraries (and in later years other collaborators ) were working on the single HEP database, and comprehensively cataloguing preprints, and by extension, the published literature, in high-energy physics.

This new PPF age, though, was not perfect. Receiving lists of new preprints and journal articles is all very well, but after several years one is left once again with a stack of paper. Trying to find a copy of a particular article or information on a certain subject can be very difficult. The Internet allowed a further development.

The problem of finding just the right preprints from this ever growing body of research was not an issue for the SLAC community. The SPIRES database allowed searches by date, author, institution, report number, title, keywords (assigned by DESY) and a number of other fields. Simply by logging in to their computer, SLAC physicists were able to search through thousands of papers, and the paper PPF list would serve as a newspaper, with a similar longevity. With the rise of the Internet anyone in the world could now access this service. A program, QSPIRES , developed in 1985 by SLAC database systems developer George Crane, allowed physicists to search the SPIRES database using E-mail and other bitnet message protocols (which actually had priority over E-mail and hence were faster). One would send off an E-mail query to the database and a swift reply would follow. This made an enormous difference for people outside SLAC, especially those outside the U.S., who had previously only been able to access SPIRES, through rare and specially issued SLAC accounts. People were able to find all the papers by a particular author or from a certain institution, or find out how many citations their work had . The PPF and PPA lists could also be sent by E-mail.

# TeX and the Single Archive

By 1990, lists of new articles were being sent to physicists by E-mail and people were searching for articles in the SPIRES database via E-mail. Everything was modern and electronic, up to the point where you actually wanted to read a paper. In this case, the reader would have to request that an author send a copy in the mail. Thus, if you were overseas it could take weeks to months before you actually saw the article; this was clearly unacceptable.

Around this time another technical innovation had sufficiently matured to obviate this problem. High energy physics articles usually contain a lot of mathematics, which makes them difficult to write with a standard set of alpha-numeric characters. In the late 1970's, Stanford computer scientist Donald Knuth had begun work on a special typesetting program, TeX (pronounced "tek" since the "x" is actually a Greek $\chi$), that could display mathematics beautifully from source containing no special characters . After 1983 when the current version was released, the language soon proved very popular with physicists and mathematicians as it gave them complete control over the production of their documents.[^2]

The underlying "TeX-file" was simply a normal, *portable* text file, in which the mathematics was written using certain rules which could then be processed *on any computer* into a viewable final output format, such as PostScript (PS), , changing

    \int_0^1\frac{e^x}{2\pi} dx

to $$\int_0^1 \frac{e^x}{2\pi} dx.$$ The PS file can then be printed out and read. Therefore, to send a paper, all you would have to do is send the tex-file via E-mail and the person at the other end could process it and print it out; another step in the communication process had entered the future. By 1990 the use of TeX in the HEP community was almost universal.

With preprint list distribution, SPIRES database searches, requests for papers, and even the transfer of these papers being done through E-mail, all stages of the process were electronic and pretty much immediate … except one. SLAC was still receiving the E-prints as "hard-copy". Some authors had begun sending their preprints as TeX files, but for some researchers this would unfortunately lead to overstuffed electronic mailboxes, echoing a problem of the Paper Days. In 1991, Paul Ginsparg of Los Alamos National Laboratory (LANL) decided to do something about this. If authors could send their papers electronically to a central repository, with author-supplied bibliographic data, a list could be sent out of each *day's* additions and the preprints could be obtained directly from this archive . Here, TeX proved a godsend, as it provided a means of storing rather small files that researchers could request and then process into the viewable PS files at their home institutions. In August 1991, the first paper was sent to the Los Alamos archive. Twenty seven papers were sent that month[^3] and physicists began receiving daily E-mails of the new papers sent to LANL.

Soon there were over 2000 subscribers to this daily E-mail notification and it gradually began to replace SLAC's PPF list (though not the SPIRES database itself) as the HEP community's "newspaper" — its most immediate source of new research information. Far from competing, though, the two services at LANL and SLAC, complemented each other nicely. SLAC continued to process hard copy preprints that were never sent to LANL and journal articles that had never been preprints, adding such information as their references, experiment numbers, authors' institutions, etc. The electronic preprints, *E-prints*, began to account for a good deal of the HEP literature, and Ginsparg's system allowed for considerable automation and efficiency in SLAC's work. A good deal of useful information could be harvested electronically from the files authors supplied to LANL (though there was still much that had to be done by hand).

This evolution from preprints to E-prints had another important facet. The E-prints were assigned a unique number of the form *archive/0007040* (eg 9501251 for the 251st paper in January 1995), which the SLAC Library immediately began to store as a new element in each record in the SPIRES database. In the past "report numbers" had been assigned to preprints by the author's home institution in a less than completely systematic way and with a variety of conventions, which made the tracking of citations too difficult for the SPIRES database. Thus, as mentioned earlier, only citations to published journal articles were recorded, and the many citations made to an author's paper while it was an unpublished preprint were lost. The standardised E-print number changed all that, since these numbers soon started appearing in reference lists. After some work, these citations to unpublished works could be registered. It is worth noting that the most cited HEP article in 1998 spent some ten months as an unpublished E-print . The next step was to ensure that the database recognised that the published article and the E-print were essentially the same paper, and thus, that the true number of citations was the sum of the citations of both. In some instances this lead to double counting, which needed special programming to eliminate, as it became fashionable to include *both* the E-print number and the journal reference when citing an article!

# The Mouse That Clicked

With the SPIRES database and the LANL archive you could find a particular paper and be reading it within minutes merely by sending E-mail, after just a little work to process the TeX file and print it out. A lot had changed in 20 years since the SLAC Library first started sending out a list of all the preprints it had received. However, the guiding philosophy behind computing is that no work a person has to do can ever be "too little." A new technology was about to make the research process even easier.

In September 1991, SLAC physicist Paul Kunz, was visiting CERN, Geneva, where he was shown the infant World Wide Web, which allowed information on different computers to be accessed in a very user-friendly (and now completely familiar) manner. Realising the Web interface would be perfect for searching the SPIRES database he told then librarian and SPIRES database manager, Louise Addis, about it upon his return to SLAC. She was quick to recognise the potential of this and by December SLAC had the first WWW server in the Western Hemisphere.[^4] The E-mail access-system became obsolete as more and more physics departments installed the free software necessary to reach the SPIRES WWW interface . A group was formed at SLAC, the WWWizards to provide help with WWW technology .

The key feature of the WWW is linking. Therefore, when displaying search results in the SPIRES database, a number of things could be linked to each record. One of these was the full-text at Los Alamos which was being stored in a minimal fashion as TeX source. As mentioned before, TeX files need to be processed before they can be printed out. In the early stages, many of the source files contained errors or idiosyncrasies that made them difficult to process. In order to enter all the relevant information into the SPIRES database the SLAC Library staff needed to print out each paper, and so would process the TeX files sent to Los Alamos into PS files. At first this was done manually, which soon grew to be very time consuming. Luckily, the SLAC Library was helped by Paul Mende from Brown University who created an automatic procedure. Accompanying figures could be faxed to the SLAC Library, where they would be converted via NextFAX to a PS file. The PS files (of both text and figures) were then stored on a server at SLAC. The record in the SPIRES database would link to both the Los Alamos entry, where the TeX source file was stored, and the SLAC PS server, where the user could view the PS file.

Around the same time a similar script to convert TeX source to PS had been written at Los Alamos by Tanmoy Bhattacharya, a physicist specialising in the heavily computational subfield of Lattice Gauge Theory. Initially for personal use, this PS generating code was built into a powerful tool by Bhattacharya, Rob Hartill and Mark Doyle. By mid 1995, the LANL archive was producing its own PS files thereby freeing the SLAC Library from this task. Ultimately PS generation was incorporated into the LANL submission process, and it came with an ultimatum: if your paper couldn't be processed, it wouldn't be accepted. Many authors were shocked to find that the TeX files they had previously sent to Los Alamos could not be easily processed, thus making their work less accessible, and welcomed this checking procedure. With the ability to generate PS from directly from TeX source, PS that had been generated from TeX was no longer acceptable, and by late 1996 submissions of such would be rejected . For whatever reasons, this led to some grumbling, but faced with being denied the opportunity to post to LANL, there was no real alternative to complying.[^5]

With all this in place the SPIRES database was linking to reliable PS files generated and stored at the Los Alamos archive. The LANL archive would also link to both references and citations stored in the SPIRES database. This formed perhaps the first publically available full-text, cross-linked scientific database. Paper had been eliminated from the process and the virtual library was born .

# The Many Hands Interpretation

The model of authors sending their electronic texts to LANL has pushed the effort of publication onto the researchers themselves, who benefit from an inexpensive, immediate, wide-scale dissemination of their work. This idea has led many computer scientists to ponder automated systems for indexing and retrieving these full-text papers. How much of the work of collecting information on the literature can be realistically facilitated by the authors? The combination of the Web, the SLAC Library's automated systems and the LANL archives provides an interesting testing ground for this question.

The goal of the LANL archive is to be as automated as possible, so that it can exist without administrative intervention (as opposed to the SPIRES database). To this end it has a number of checks to ensure that all submissions meet the entry requirements (one we have already discussed ensures the TeX file is successfully processed into PS). The basic bibliographic information (author, title, etc.) data is supplied by the users in a neatly structured format that can be downloaded into the SPIRES database automatically. This raw information requires only minor attention from the Library staff who, among other things, ensure spelling consistencies, and add in the author affiliations from the INSTITUTIONS database.

Far and away the most time consuming part of this is collecting the references of each paper, from which the citation searching is built . As citation results reflect some degree of professional accomplishment, physicists tend to be rather interested in them and a good deal of their correspondence with the SLAC Library concerns omissions or mistakes in the reference lists (which is exacerbated by the posting to LANL of revised versions of E-prints, which often contain additional references). Originally these reference lists were typed into the database, but as another happy spin-off of TeX about 90% of them can now be extracted from the author's file. Unfortunately, this only makes a huge job large, as they still have to be checked by the Library staff, and authors regularly confound the reference extracting program by adopting imaginative new ways of writing the journal-volume-page sequence, or simply making errors in a reference.

Obviously the LANL model suggested there should be some way to place the burden of constructing and checking the reference list on the authors themselves (as they do with bibliographic information), before they send their paper out to the world. Two things really stood in the way of this. First, the original program that extracted the references was not definitive enough for this sort of user-side checking, as there was no real way of specifying which references the code would extract, and which it would not. The second problem was that with physicists being rather busy, any new system would have to result in less, not more, work for it to be widely adopted.

Once again, TeX helped with the solution. The SPIRES team developed and offered a new service: the database would display records in TeX format, with the exact information needed to construct the reference as an additional tag. The authors then had a very simple process of cut and paste to create their reference lists, and the tag would sit invisibly in the TeX file (not appearing in the final PS file) waiting to be extracted at the SLAC library. The strict structure of the CITATION tag and the simplicity of its extraction allowed for the creation of a checking program that authors could use, for example:

    \bibitem{O'Connell:1997}
    H.~B.~O'Connell,
    ``Recent developments in $\rho\!-\!omega$ mixing,''
    Austral.\ J.\ Phys.\ {\bf 50}, 255 (1997)
    [hep-ph/9604375].
    %%CITATION = HEP-PH 9604375;%%

which becomes:
H. B. O'Connell, "Recent developments in $\rho\!-\!\omega$ mixing," Austral. J. Phys. **50**, 255 (1997) \[hep-ph/9604375\].
So far over five hundred papers have been written using this system.

A second way we've been experimenting with author–supplied data is asking them to help us with another thing that is of great interest to them: missing papers in our database. Traditionally adding them in was another time consuming exercise that would require Library staff to get the journal from the shelves and type in all the relevant information. Using a web form that authors can fill out themselves allows a paper to be added to the database simply by cutting and pasting on the part of the Library staff.

By inviting our users to help us in maintaining the database and automating as many processes as possible, we have been able to accomplish a lot of additional work without spending significant amounts of extra time. This would not be possible, though, without the eagerness and attention to detail that characterises the Physics community.

# All the news that's fit to link

The WWW created a unique opportunity for the SPIRES database by providing a system that would conveniently attach to any record all the related information both inside the database and around the world in a very compact manner. Over the past six years we have worked on making this process as efficient as possible both in terms of computer programming and the output display our users see . I shall now discuss how this works.

A single record in the Literature database might have certain *basic* elements such as

-   Title

-   Author (and author's institution) or corporate author

-   Date

-   Publication note

but the article could also be *hyperlinked* to

-   the appropriate record in the EXPERIMENTS database

-   full text of article at the journal server or the E-print server

-   the references of the article

-   other articles that cite it

-   experimental data in the REACTIONS database,

-   the record for the conference at which the paper was given

allowing the researcher to selectively explore any particular related item. Thus an original record only ten lines long becomes a gateway to information stored around the globe.

The construction of these links requires special care and co-operation with outside services. One particular case that followed an evolutionary process was linking to the published version of the document on journal home pages. Optimally, this URL would point to a unique "abstract page," rendered in HTML so that the link to the journal server can be fast. Once there the user can be presented with all the Journal services such as full text in a variety of formats (the most common being PS or the newer Portable Document File (PDF)). From the point of view of the publisher, this means that a URL has to be found for every article. From our point of view this URL should be calculable from the information we already have about the record, as it then permits us to run in the URL automatically. Articles have always been cited through the journal-volume-page (JVP) convention, so it made sense that a number of major publishers, including the American Physical Society and Elsevier adopted a URL scheme based on these three elements. Setting up such a system does present something of a technical challenge for the journals, but is well worth it in terms of presenting the simplest possible interface to the outside world and providing reliable access to their wares.

The journals also link back into our database, mainly for the references of the paper (though publishers have proven less eager to link to our database record of the actual paper). In an effort which I hope will become more widely adopted, we have worked closely with the American Physical Society's *Physical Review* to share the bibliographic data (including reference lists) and either update an existing record in our database, or add a new record, when they publish new articles.

Where possible we have also tried to link to other literature databases. Some, like the CERN Library's database , have a significant overlap with our system (there are currently over fifty thousand two-way links between SPIRES and CERN). Others, such as Harvard's *Astrophysics Data System* and the American Mathematical Society's *Mathematical Reviews* (AMSMR) cover a much smaller set of papers in our database, but from a different perspective and connect the paper to other academic disciplines. Here, once again, we need to work with the other databases to ensure maximal efficiency and reliability in this linking. Our links to the Harvard database use a JVP scheme, while those to CERN and AMSMR databases use unique record keys. SPIRES can be linked to either way. There is also the consideration that the link should bring additional information, rather than just repeating what record in the SPIRES database (or else why bother?). Therefore, we apply these links selectively.

# So why HEP?

Why did the High Energy Physics community provide such a fertile ground for this particular aspect of the so-called "Information Revolution"? We have touched on the reasons throughout this talk, but it is useful, perhaps, to summarise them to understand how this model might be more widely adopted (or why it might *not* be).

The first thing is that it has not been a revolution so much as an evolution. Particle physicists have always been compulsive communicators. Generally free of commercial or governmental restrictions, the habit of sending out advanced copies of work goes back to the days of carbon paper. They also have a long standing tradition of international collaboration, so any new advances in communication technology are eagerly adopted, and the high technological literacy facilitates this. This is especially true of the experimental community whose collaborations can have hundreds of members all over the world and whose experiments depend on high speed computer network connections. It is worth noting that the two SLAC physicists most closely involved in spinning the web at SLAC, Paul Kunz and Tony Johnson are experimentalists, though I should also point out that Paul Ginsparg of Los Alamos is a theorist.

In this environment the sheer quantity of existing literature necessitated creating the organisation begun by the SLAC Library in the late 1960's. Through using this PPF list service, physicists became accustomed to the idea of a central "clearing house" for preprints. As the use of E-mail spread through the Physics community, particularly the High Energy sector, the SPIRES E-mail interface was introduced in the mid 80's and before long was being used by thousands of researchers in over forty countries. The TeX typesetting language, favoured by those needing to write complex mathematics, allowed one to send simple ASCII files using E-mail, thus giving rise to the electronic distribution of preprints, which was centralised by Ginsparg's E-print archive. The World Wide Web then laid the ground for a highly powerful way to integrate the various facets of research communication.

It should be acknowledged that these are rather special conditions. It is instructive to note that in other scientific communities where TeX is less heavily used, the adoption of the LANL archives has been noticeably slower, though in principle there is nothing to stop them (despite being built with TeX in mind, the archive is not dependent on it, and will accept submissions in a variety of formats). The TeX fluent mathematics community, on the other hand, though quick to adopt the new system, perhaps lacked the detailed computing knowledge to play a leading role.

For different research communities there are other obstacles. Faced with the incredibly rapid and wholehearted adoption of this E-print scheme by high energy physicists, the traditional journal publishers found it best to not to offer any real resistance. They had nothing to fear, really, as over seventy percent of HEP papers sent to Los Alamos are eventually published in journals and another twenty percent appear in conference proceedings. In other fields, with newer archives, this is still something of a touchy subject, such as in the biological and medical research communities, as can be seen in the Biochemical Society's response to the NIH's proposal for a preprint server for the life sciences (further reflections on eprints and the life sciences can be found at ).

# Conclusion

In this paper I have described the evolution of paperless publishing in a particular academic community and how special circumstances conspired to make this happen faster and more comprehensively than in any other field. In doing this I have hoped to convey the sense that this was indeed a process of evolution rather than revolution.

It may be interesting to speculate what the future holds. Will more and more of what used to be called "bibliographic data," become author–supplied "metadata" which annotates each article with various related information? Can this metadata be shared to create cross-database, and cross-discipline, linking ? Our experience shows that if you make the procedure simple enough, and offer an advantage to the authors for doing it (such as the ability to check their reference list), they *will* use it. The challenge now is to bring this to other scholarly communities .

We have begun in the SPIRES-HEP database to exploit Web technology to handle information beyond the traditional bibliographic record or what one could obtain from a print version (such as the full text and reference list). These other facets range from the paper's citation list, which is constantly growing (if the author is fortunate), to the home-page of the experimental collaboration that wrote the paper, to what information researchers might really want out of a paper, in a form far more useful than could ever be delivered in paper (such as a computer file of experimental data or 3-D computer images created by the author).

The fluidity of this new, electronic means of publication also has implications for the traditional research process. At what point does an author "give the final answer" and the paper become "set in stone"? The LANL archive, has neatly addressed this issue in an appropriate manner, by archiving each version of an E-print with a date stamp, documenting any modifications in the paper. Once the article is published in a journal, it is a finished piece of work, with any subsequent alterations handled via errata. Other authors may *comment* on the published paper, to which the author may *reply*, but the comment and reply are treated as new articles, and the original is left unchanged. One journal, however, has broken away from the traditional model. *Living Reviews in Relativity* exists solely on the Web, and allows authors to constantly revise their articles. In part this is due to the pedagogical nature of the journal; it seeks to provide reviews that aid learning, rather than publish original research.

Well aware of the difficulties of predicting anything in the Internet world, I have merely tried, in concluding, to state some of the current trends in electronic publication.

**Acknowledgements**

I would like to thank Louise Addis, Paul Ginsparg, Marek Karliner, Pat Kreitz, Michael Peskin, Hartmut Preissner and Ann Redfield for reading this manuscript and offering helpful comments and suggestions. I have also enjoyed conversations on the topics mentioned here with Richard Dominiak, Mark Doyle, Paul Ginsparg, Tony Johnson, John Jowett and Paul Kunz.

[^1]: Work supported by the US Department of Energy contract DE-AC03-76SF00515.
    *Invited talk at the AAAS Annual Meeting and Science Innovation Exposition, February 17-22, 2000 in Washington, DC*

[^2]: This paper has been written using a derivative of TeX

[^3]: In August 1999, over 800 particle physics papers were sent to the LANL archive.

[^4]: There was no Initial Public Offering.

[^5]: An interesting and oftentimes amusing timeline of the development of the Los Alamos archive can be found at .


author: Mariafelicia De Laurentis$^\dag$, Salvatore Capozziello$^\dag$; Luca Izzo$^\dag$ $^\ddag$
title: Stochastic background of gravitational waves "tuned" by $f(R)$ gravity

The cosmological background of gravitational waves can be tuned by Extended Theories of Gravity. In particular, it can be shown that assuming a generic function $f(R)$ of the Ricci scalar $R$ gives a parametric approach to control the evolution and the production mechanism of gravitational waves in the early Universe.

# Introduction

In the last thirty years several shortcomings came out in Einstein General Relativity (GR) and people began to investigate whether it is the only theory capable of explaining gravitational interactions. Such issues sprang up in Cosmology and Quantum Field Theory. In the first case, the Big Bang singularity, the flatness and horizon problems led to the result that the Standard Cosmological Model is inadequate to describe the Universe in extreme regimes. Besides, GR is a *classical* theory which does not work as a fundamental theory. Due to these facts and to the lack of a self-consistent Quantum Gravity theory, alternative theories have been pursued. A fruitful approach is that of Extended Theories of Gravity (ETGs) which have become a sort of paradigm based on corrections and enlargements of GR.

These theories have acquired interest in Cosmology owing to the fact that they "naturally" exhibit inflationary behaviors . Recently, ETGs are playing an interesting role in describing the today observed Universe. In fact, the amount of good quality data of last decade has made it possible to shed new light on the effective picture of the Universe. In particular, the *Concordance Model* predicts that baryons contribute only for $\sim4\%$ of the total matter - energy budget, while the exotic *cold dark matter* represents the bulk of the matter content ($\sim25\%$) and the cosmological constant $\Lambda$ plays the role of the so called "dark energy" ($\sim70\%$). Although being the best fit to a wide range of data, the $\Lambda$CDM model is severely affected by strong theoretical shortcomings that have motivated the search for alternative models. Dark energy models mainly rely on the implicit assumption that GR is the correct theory of gravity indeed. Nevertheless, its validity on the larger astrophysical and cosmological scales has never been tested, and it is therefore conceivable that both cosmic speed up and dark matter represent signals of a breakdown in GR .

From an astrophysical viewpoint, ETGs do not require to find out candidates for dark energy and dark matter at fundamental level; the approach starts from taking into account only the "observed" ingredients (i.e. gravity, radiation and baryonic matter); this is in agreement with the early spirit of GR which could not act in the same way at all scales. In fact, GR has been definitively probed in the weak field limit and up to Solar System scales. However, a comprehensive effective theory of gravity, acting consistently at any scale, is far, up to now, to be found, and this demands an improvement of observational datasets and the search for experimentally testable theories. A pragmatic point of view could be to "reconstruct" the suitable theory of gravity starting from data. The main issues of this "inverse " approach is matching consistently observations at different scales and taking into account wide classes of gravitational theories where "ad hoc" hypotheses are avoided. In principle, the most popular dark energy models can be achieved by considering $f(R)$ theories of gravity and the same track can be followed to match galactic dynamics . This philosophy can be taken into account also for the cosmological stochastic background of gravitational waves (GW) which, together with CMBR, would carry, if detected, a huge amount of information on the early stages of the Universe evolution.

# Gravitational waves in f(R)-gravity

In this paper, we face the problem to match a generic $f(R)$ theory with the cosmological background of GWs. GWs are perturbations $h_{\mu\nu}$ of the metric $g_{\mu\nu}$ which transform as 3-tensors. The GW-equations in the transverse-traceless gauge are $$\square h_{i}^{j}=0\label{eq: 1},$$ Latin indexes run from 1 to 3. Our task is now to derive the analog of Eqs. () from a generic $f(R)$ given by the action $$\mathcal{A}=\frac{1}{2k}\int d^{4}x\sqrt{-g}f(R)\label{eq:2}\,.$$ From conformal transformation, the extra degrees of freedom related to higher order gravity can be recast into a scalar field $$\widetilde{g}_{\mu\nu}=e^{2\Phi}g_{\mu\nu}\qquad \mbox{with}
\qquad e^{2\Phi}=f'(R)\,.\label{eq:3}$$ Prime indicates the derivative with respect to $R$. The conformally equivalent Hilbert-Einstein action is $$\mathcal{\widetilde{A}}=\frac{1}{2k}\int\sqrt{-\tilde{g}}d^{4}x\left[\widetilde{R}+
\mathcal{L}\left(\Phi\mbox{,}\Phi_{\mbox{;}\mu}\right)\right]\label{eq:4}$$ where $\mathcal{L}\left(\Phi\mbox{,}\Phi_{\mbox{;}\mu}\right)$ is the scalar field Lagrangian derived from $$\widetilde{R}=e^{-2\Phi}\left(R-6\square\Phi-6\Phi_{;\delta}\Phi^{;\delta}\right)\,.\label{eq:6}$$ The GW-equation is now $$\widetilde{\square}\tilde{h}_{i}^{j}=0\label{eq:7}$$ where $$\widetilde{\square}=e^{-2\Phi}\left(\square+2\Phi^{;\lambda}\partial_{;\lambda}\right)\label{eq:9}\,.$$ Since no scalar perturbation couples to the tensor part of gravitational waves, we have $$\widetilde{h}_{i}^{j}=\widetilde{g}^{lj}\delta\widetilde{g}_{il}=e^{-2\Phi}g^{lj}e^{2\Phi}\delta
g_{il}=h_{i}^{j}\label{eq:8}$$ which means that $h_{i}^{j}$ is a conformal invariant.

As a consequence, the plane-wave amplitudes $h_{i}^{j}(t)=h(t)e_{i}^{j}\exp(ik_{m}x^{m}),$ where $e_{i}^{j}$ is the polarization tensor, are the same in both metrics. This fact will assume a key role in the following discussion.

# Cosmological stochastic background of gravitational waves

In a FRW background, Eq.() becomes $$\ddot{h}+\left(3H+2\dot{\Phi}\right)\dot{h}+k^{2}a^{-2}h=0\label{eq:10}$$ being $a(t)$ the scale factor, $k$ the wave number and $h$ the GW amplitude. Solutions are combinations of Bessel's functions. Several mechanisms can be considered for the production of cosmological GWs. In principle, we could seek for contributions due to every high-energy process in the early phases of the Universe.

In the case of inflation, GW-stochastic background is strictly related to dynamics of cosmological model. This is the case we are considering here. In particular, one can assume that the main contribution to the stochastic background comes from the amplification of vacuum fluctuations at the transition between the inflationary phase and the radiation era. However, we can assume that the GWs generated as zero-point fluctuations during the inflation undergo adiabatically damped oscillations $(\sim 1/a)$ until they reach the Hubble radius $H^{-1}$. This is the particle horizon for the growth of perturbations. Besides, any previous fluctuation is smoothed away by the inflationary expansion. The GWs freeze out for $a/k\gg H^{-1}$ and reenter the $H^{-1}$ radius after the reheating. The reenter in the Friedmann era depends on the scale of the GW. After the reenter, GWs can be detected by their Sachs-Wolfe effect on the temperature anisotropy $\bigtriangleup T/T$ at the decoupling. If $\Phi$ acts as the inflaton, we have $\dot{\Phi}\ll H$ during the inflation. Adopting the conformal time $d\eta=dt/a$, Eq. () reads $$h''+2\frac{\chi'}{\chi}h'+k^{2}h=0\label{eq:16}$$ where $\chi=ae^{\Phi}$. The derivation is now with respect to $\eta$. Inside the radius $H^{-1}$, we have $k\eta\gg 1.$ Considering the absence of gravitons in the initial vacuum state, we have only negative-frequency modes and then the solution of () is $$h=k^{1/2}\sqrt{2/\pi}\frac{1}{aH}C\exp(-ik\eta)\,.\label{eq:18}$$ $C$ is the amplitude parameter. At the first horizon crossing $(aH=k)$ the averaged amplitude $A_{h}=(k/2\pi)^{3/2}\left|h\right|$ of the perturbation is $$A_{h}=\frac{1}{2\pi^{2}}C\,.\label{eq:19}$$ When the scale $a/k$ becomes larger than the Hubble radius $H^{-1}$, the growing mode of evolution is constant, i.e. it is frozen. It can be shown that $\bigtriangleup T/T\lesssim A_{h}$, as an upper limit to $A_{h}$, since other effects can contribute to the background anisotropy. From this consideration, it is clear that the only relevant quantity is the initial amplitude $C$ in Eq. (), which is conserved until the reenter. Such an amplitude depends on the fundamental mechanism generating perturbations. Inflation gives rise to processes capable of producing perturbations as zero-point energy fluctuations. Such a mechanism depends on the gravitational interaction and then $(\bigtriangleup T/T)$ could constitute a further constraint to select a suitable theory of gravity. Considering a single graviton in the form of a monochromatic wave, its zero-point amplitude is derived through the commutation relations: $$\left[h(t,x),\,\pi_{h}(t,y)\right]=i\delta^{3}(x-y)\label{eq:20}$$ calculated at a fixed time $t$, where the amplitude $h$ is the field and $\pi_{h}$ is the conjugate momentum operator. Writing the Lagrangian for $h$ $$\widetilde{\mathcal{L}}=\frac{1}{2}\sqrt{-\widetilde{g}}\widetilde{g}^{\mu\nu}h_{;\mu}h{}_{;\nu}\label{eq:21}$$ in the conformal FRW metric $\widetilde{g}_{\mu\nu}$, where the amplitude $h$ is conformally invariant, we obtain $$\pi_{h}=\frac{\partial\widetilde{\mathcal{L}}}{\partial\dot{h}}=e^{2\Phi}a^{3}\dot{h}\label{eq:22}$$

Eq. () becomes $$\left[h(t,x),\,\dot{h}(y,y)\right]=i\frac{\delta^{3}(x-y)}{a^{3}e^{2\Phi}}\label{eq:23}$$ and the fields $h$ and $\dot{h}$ can be expanded in terms of creation and annihilation operators. The commutation relations in conformal time are $$\left[hh'^{*}-h^{*}h'\right]=\frac{i(2\pi)^{3}}{a^{3}e^{2\Phi}}\,.\label{eq:26}$$ From () and (), we obtain $C=\sqrt{2}\pi^{2}He^{-\Phi}$, where $H$ and $\Phi$ are calculated at the first horizon-crossing and, being $e^{2\Phi}=f'(R)$, the relation $$A_{h}=\frac{H}{\sqrt{2f'(R)}}\,, \label{eq:27}$$ holds for a generic $f(R)$ theory.

# Conclusion

The central result (eq. ()) of this paper deserves some discussion. Clearly the amplitude of GWs produced during inflation depends on the theory of gravity which, if different from GR, gives extra degrees of freedom. On the other hand, the Sachs-Wolfe effect could constitute a test for gravity at early epochs. This probe could give further constraints on the GW-stochastic background, if ETGs are independently probed at other scales.


abstract: We present a method of detecting sequence defects by supercoiling DNA with magnetic tweezers. The method is sensitive to a single mismatched base pair in a DNA sequence of several thousand base pairs. We systematically compare DNA molecules with 0 to 16 adjacent mismatches at 1 $M$ monovalent salt and 3.5 pN force and show that, under these conditions, a single plectoneme forms and is stably pinned at the defect. We use these measurements to estimate the energy and degree of end-loop kinking at defects. From this, we calculate the relative probability of plectoneme pinning at the mismatch under physiologically relevant conditions. Based on this estimate, we propose that DNA supercoiling could contribute to mismatch and damage sensing *in vivo*.
author: Andrew Dittmore$^1$, Sumitabha Brahmachari$^2$, Yasuhara Takagi$^1$, John F. Marko$^{2,3}$, and Keir C. Neuman
date: 2024-09-29
title: Supercoiling DNA locates mismatches

With increasing imposed torsion, a thin elastic rod builds torque until it reaches a limit point, and then abruptly loops into a self-contacting, plied structure called a plectoneme . Although in a macroscopic system this abrupt buckling is expected to occur at a site of elastic discontinuity , the analogous situation of a base-pair mismatch defect in a single DNA molecule is less clear: Thermal fluctuations could mask the effects of the defect in a DNA molecule of several thousand base-pairs. Alternatively, the bending energy decrease could overcome entropy, in which case the position of the defect would determine where the plectoneme forms. In addition to being a unique physical system in which to explore the role of defects in the torsional buckling of a fluctuating elastic rod, understanding how defects influence DNA supercoiling has potential biological implications. Whereas DNA supercoiling is known to be a major determinant of the large-scale organization of cellular DNA , we propose that torque may also act globally to facilitate the detection of defect sites among millions of normal base pairs .

Supercoiling potentially provides a mechanism of defect sensing by localizing a defect at a plectoneme tip, where the duplex will be further destabilized by bending stress. Previous work provides precedent for the idea of plectoneme pinning at DNA defects. Brutzer *et al.* showed overwound DNA buckles and forms a plectoneme that is pinned at a permanent kink . Computations by Matek *et al.* demonstrated that energy is minimized by the colocalization of plectoneme end-loops and regions of unpaired bases as "tip bubbles." Ganji *et al.* used an intercalator based method for creating fluorescent plectonemes and showed a preference for plectoneme pinning at the location of a 10 bp mismatch.

Here we develop a framework to quantify how defects influence DNA supercoiling. By isolating individual plectonemes in DNA constrained by force and ionic screening, we find that even a single mismatch site specifies the location of a plectoneme. Statistical-mechanical calculations allow us to extend these results to physiologically relevant conditions.

Our experimental method is based on a magnetic tweezers supercoiling assay and makes use of the abrupt, discontinuous buckling of DNA at the onset of plectoneme formation (Fig. ). If a single pinned plectoneme forms at a defect positioned near a surface, it is prevented from lengthening as it impinges on the surface. This results in a second, abrupt buckling transition that occurs at double the distance of the defect from the surface. Thus the DNA extension at re-buckling reports the position of the defect.

To demonstrate this measurement approach, we tested torsionally constrained, 6 kb DNA molecules with $n=0$, 1, 2, 4, or 16 adjacent mismatches (Fig. 2). We positioned the mismatch site roughly 8$\%$ from one end of the DNA using a cassette based single-strand nicking template generated by PCR and simultaneously ligated a mismatch-containing oligonucleotide and handles for torsional constraint . Upon supercoiling we observed re-buckling at the expected DNA extension $2\times8\%=16\%$ below the maximum. This confirms a single plectoneme formed and was pinned by the defect for all defect sizes, even a single mismatch ($n$=1). We do not observe re-buckling in DNA molecules with a centered defect, or in intact DNA ($n$=0). The re-buckling signal also disappears below the salt-force phase boundary demarcating the edge of the single-plectoneme regime (Fig. ); at lower forces or ionic strengths, multiple plectonemes are expected , and the pinned plectoneme may begin to exchange length with other plectonemes or diffuse . All curves shown in Fig. 2 were collected at a fixed force of $f=3.6\pm0.4$ pN and 1 $M$ monovalent salt; these conditions favor a single plectoneme , which we find to be stably pinned at the mismatch.

*Defects cause kinking of the pinned plectoneme end-loop* – We focus on the initial plectoneme formation, which occurs at the mismatch whenever re-buckling is observed (Fig. 2). The abrupt extension drop, $\Delta X$, upon initial buckling decreases quadratically with defect size, $n$ (Fig. 2a,c). This is consistent with the expectation that the plectoneme end-loop decreases in size and is sharply bent at its apex to a degree that depends on the extent of the mismatch (Fig. 1). Crucially, defects reduce the critical buckling torque $\Gamma_{\rm c}$ occurring at linking number ${\rm Lk}^\dagger$ .

To quantify the critical buckling torque $\Gamma_{\rm c}$ as a function of defect size, we define $\Gamma_{\rm c}$ relative to an internal reference torque, $\Gamma^*$, which is constant in the DNA after re-buckling and estimated using the analytic formula derived by Clauvelin, Audoly, and Neukirch (Fig. 1): $\Gamma^*=21.0\pm0.5$ pN$\cdot$nm. This analysis also provides estimates of the structural parameters of the plectoneme : $\alpha=23.8\pm0.4^{\circ}$ is the ply angle and $r=1.7\pm0.1$ nm is its radius, consistent with a high degree of electrostatic screening .

For intact DNA ($n=0$) we assume that $\Gamma_{\rm c}-\Gamma^*$ is well approximated by the discrete torque drop $\Delta \Gamma_0$, which can be derived from the equilibrium buckling probability obtained from analysis of the extension fluctuations (Fig. 2a). Using measurement of $K=P_{unbuckled}/P_{buckled}$ from the relative probabilities of the unbuckled and buckled states (Fig. 2a), we estimate the free energy difference $\Delta F=-k_BT\ln{K}$ and note that torque is the partial derivative of free energy with respect to rotation angle, $2\pi {\rm Lk}$: $\Delta\Gamma=-k_BT/(2\pi)(\partial{\ln{K}}/\partial{{\rm Lk}})$. Our measurements at controlled ${\rm Lk}$ with all other differential quantities in the potential fixed permit straightforward evaluation of the partial derivative and model-free estimation of $\Delta\Gamma$. Under the assumption that $\Gamma_{\rm c}-\Gamma^*$ is proportional to the linking number difference, ${\rm Lk}^\dagger-{\rm Lk}^*$ (Fig. 1, Fig. 2d), we calculate $\Gamma_{\rm c}(n)\approx\Gamma^*+\Delta\Gamma_0({\rm Lk}^\dagger_n-{\rm Lk}^*_n)/ ({\rm Lk}^\dagger_0-{\rm Lk}^*_0)$. Together, these estimates provide the critical torque over the range of abrupt buckling ($n<5$) (Fig. 2d). For larger defects ($n>5$) abrupt buckling is not observed. We conclude that for more than five adjacent mismatches there is no end-loop and abrupt buckling is replaced by a kinetic pathway without a transition barrier. Consistent with this, the $n=16$ curve (Fig. 2b) shows a continuous transition from the buckling point, ${\rm Lk}^*$; *i.e.*, the plectoneme forms once $\Gamma\approx\Gamma^*$.

For $n=16$, $\Delta\Gamma=0$. Surprisingly, $\Delta \Gamma$ for each of the smaller defects is nearly constant with an average value of $\Delta\Gamma=2.92\pm 0.08$ pN$\cdot$nm for $n=0$ to $n=4$. This indicates that the torque initially drops below $\Gamma^*$ but approaches $\Gamma^*$ with increasing ${\rm Lk}$. This undershoot is consistent with a more compact initial plectoneme structure due to the decreased bending energy of the end-loop.

*Kinking lowers the end-loop energy* – These measurements establish a relationship between end-loop energy and defect size. The end-loop energy can be related to geometry in a simple planar loop model . The size of the end-loop is obtained via minimizing the elastic energy, $$\beta\mathcal{E}_\gamma=\epsilon\frac{A}{\gamma}+\mu\frac{\gamma}{A},$$ where $\beta$ is the inverse thermal energy, $A$ is the bending persistence length, $\mu\equiv\beta A f$, and $\epsilon$ is a numerical prefactor related to the elastic energy of a loop of size $\gamma=A\sqrt{\epsilon/\mu}$. Notably, $\epsilon\approx16$ for a "teardrop" shape, and $\epsilon\approx4$ for a loop with its apex kinked by 90$^\circ$ (Fig. 4a) . Within the context of this model, the end-loop energy is related to the extension change $\Delta X(n)$ (Fig. 2c) through $\gamma(n)$. The data in Fig. 2 correspond to plectonemes with end-loops ranging from zero kinking at $n=0$ to a collapsed loop for large defects ($n>5$), which we assume to be kinked to a limiting interior apex angle of $2\alpha$. By fixing these limits and fitting to $\Delta{X(n)}$ (Fig. 2c), we estimated $\epsilon$ vs $n$ for our experimental data set (Fig 4b).

*End-loop kinking model predicts a preference for pinned plectonemes* – Outside the regime where there is a single pinned plectoneme (Fig. 3) we lose the re-buckling signal and our experimental method cannot unambiguously determine the position or number of plectonemes: Plectoneme pinning at the mismatch becomes probabilistic rather than deterministic . To obtain these probabilities, we turn to theory, employing the energy-scaling parameter $\epsilon$ obtained from experiments (Fig. 4b).

Building on a detailed mesoscopic model described previously , we consider a DNA molecule of contour length $L$ with defect positioned at $L$'. We sum over all possible sizes and numbers of plectonemes to construct a canonical partition function that includes $m\in[0,1,2\dots]$ plectonemes with non-kinked loops and the $m^\dagger\in[0,1]$ pinned plectoneme with its end-loop kinked by the defect.

In calculating the partition sum, both energetic and entropic contributions are accounted for in the fixed force and fixed linking number ensemble. The DNA is partitioned into a force-extended fraction and plectonemic fraction. The extended fraction includes the free energy of twisting as well as the total energy associated with force-extension and lateral fluctuations. The plectonemic fraction includes twist energy, superhelical bending, electrostatic repulsion, and the elastic energy contribution of plectoneme end-loops ($\epsilon<16$ at a defect and $\epsilon_0=16$ elsewhere); in addition, we explicitly include contributions to the free energy of fluctuations in the plectoneme , and entropy correction factors to account for plectonemes that are mobile and exchange length. We minimize free energy to calculate the equilibrium values of the number of mobile plectoneme domains $\langle  m\rangle$, and occupancy of the pinned plectoneme ($0\leq\langle m^\dagger\rangle\leq1$) at the defect. We will describe the details of this calculation and its range of predictive results in a forthcoming article.

In essence, the decrease in bending energy at the defect competes with the loss of entropy associated with a defect-pinned plectoneme. The pinned plectoneme is expected to be present in the fluctuating system with a fractional occupancy that approaches 1 only in the limit of high tension and ionic screening (Fig. 3). Although at lower force and ionic strength, multiple plectonemes may be present and the defect causes only a small change to the global probability of plectoneme formation, a plectoneme pinned at the defect represents a large local change, which we can now calculate, even if its fractional occupancy is $\ll1$.

To assess whether supercoiling could influence DNA damage sensing and repair, we focus on conditions most relevant to DNA in the bacterial cell and calculate $\langle m\rangle$ and $\langle m^\dagger\rangle$ for an average $10$ kb topological domain at supercoiling density 0.05 and $0.2~M$ monovalent salt. Whereas the $\langle m\rangle$ diffusing plectonemes are randomly distributed across the $N=10^4$ base pairs, the pinned plectoneme with fractional occupancy $\langle m^\dagger\rangle$ occurs only at the defect ($\sim1$ base pair). We therefore compare the ratio of probabilities, $\langle m^\dagger\rangle/(\langle m\rangle/N)$, yielding the relative enhancement of plectoneme occupancy at the defect as a function of $\epsilon$ (Fig. 4c). We restrict our attention to small defects ($n<2$) and note that the values ($8<\epsilon<16$) are expected to cover a range of different defects on the scale of 1 bp, including mismatches, abasic sites, and other lesions . Relative to a random position on the DNA, kinking generally enhances plectoneme occupancy at the defect by a factor that depends on force and the end-loop kinking energy (Fig. 4c). The experimental phase diagram (Fig. 3) is broadly consistent with the calculated force dependence of defect-enhanced buckling in Fig. 4. At 3.5 pN $\langle m^\dagger\rangle/\langle m\rangle>1$, consistent with the observation of re-buckling at this force (Fig. 3). At the low force of 0.25 pN and for end-loop energies corresponding to a single base-pair mismatch (Fig. 4b; $\epsilon\approx13.5$ for $n=1$), physiologically relevant supercoiling promotes buckling at the defect by a factor $\sim10^2-10^3$ (Fig. 4c).

*Repair of mismatches and damaged bases begins with sensing* – The single-molecule data and theory presented here establishes that supercoiling localizes mismatches to the tips of kinked plectoneme end-loops. From a biological perspective, our data suggest a role for supercoiling in defect sensing and repair. Since plectonemes project outward from the dense bacterial nucleoid, presentation of defects at plectoneme tips could promote access to lesion sites and accelerate the diffusive search by proteins that bind and initiate repair . Furthermore, binding at the lesion is accompanied by DNA kinking and stabilization of base-flipped configurations . Since bending stress at the plectoneme end-loop favors these protein-bound structural states, tighter protein binding is expected from an energetic standpoint. The hundred-fold or greater relative enhancement of plectoneme occupancy at a mismatch under physiological supercoiling conditions could therefore promote mismatch detection by a similar or possibly larger factor.

This work was supported by in part by the Intramural Research Program of the National Heart, Lung, and Blood Institute, National Institutes of Health. Work at NU was supported by the NIH through grants R01-GM105847, U54-CA193419 (CR-PS-OC) and a subcontract to grant U54-DK107980, and by the NSF through grants MCB-1022117 and DMR-1206868.

# Supporting Information

## Supercoiling DNA locates mismatches

*DNA preparation methods* – Using the Bam-HI restriction site, we inserted the forward sequence GATCGCTGAGGCGCAGCTTCCGACTGCAGCCTGACGCCAGGGCTGAGGT after position 198 in the pET-28b plasmid. This introduced both a central Pst-I restriction site (CTGCAG) and two sites for the nicking enzyme Nt.BbvCI. Successful insertion also eliminated the Bam-HI site to allow for clean selection of the recombinant DNA prior to transfection. We amplified a region of this plasmid DNA with Phusion polymerase using the forward primer -GAACCATCACCCTAATCAAG and reverse primer -GAACAACACTCAACCCTATC. These were prepended by the overhang sequences GCTGGGTCTCGCAAC- and GCTGGGTCTCGACCA-, respectively, to introduce nonpalendromic Bsa-I cleavage sites . After purification of the PCR product, we simultaneously nicked with Nt.BbvCI and cut with Bsa-I. We then annealed a 10-fold excess of the sequence GGCGCAGCTTCCGACTGCAGCCTGACGCCAGGGCTGA to competitively hybridize with the excised DNA fragment. Upon second purification, this produced DNA with an internal 37-nt gap and unique 4-nt sticky overhangs. We ligated 0.5 kb handles containing the complementary 4-nt overhangs and functionalized with either multiple biotin or digoxigenin moieties and each of the following sequences to introduce a 1, 2, 4, or 16 bp mismatch (underscore):

-   $n=1$:   /5'Phos/TCAGCCCTGGCGTCAGGCAGCAGTCGGAAGCTGCGCC

-   $n=2$:   /5'Phos/TCAGCCCTGGCGTCAGGCACCAGTCGGAAGCTGCGCC

-   $n=4$:   /5'Phos/TCAGCCCTGGCGTCAGGGACGAGTCGGAAGCTGCGCC

-   $n=16$:  /5'Phos/TCAGCCCTGGGCAGTCCGACGTCAGCGAAGCTGCGCC

Successful insertion of the mismatch-containing sequence destroys the Pst-I site. This allowed us to finish the product with Pst-I digestion, ensuring a pure sample for testing in magnetic tweezers. Intact DNA was prepared similarly but without Nt.BbvCI or Pst-I.
*Isolation and testing of single DNA molecules using magnetic tweezers* – Each DNA molecule ($2\pm0.2~\mu$m contour length) was isolated and torsionally constrained between an antibody-coated glass surface and a $\sim1~\mu$m streptavidin-coated paramagnetic bead (Dynal) using multiple biotin or digoxigenin moieties at either extremity. A magnetic field gradient, controlled by the proximity of permanent magnets, pulled the bead upward and extended the DNA. At constant force, magnet rotation synchronously rotated the bead to add turns to the DNA, thereby controlling its topology. We performed all measurements in Fig. 2 in $10$ m$M$ Tris buffer, pH 7.5 with 1 $M$ NaCl and 0.1$\%$ Tween-20. For the data in Fig. 3, we used the same buffer conditions but varied the salt concentration.
*Comparison to traditional magnetic tweezers measurements* – In the absence of the re-buckling signal (Fig. 1), traditional magnetic tweezers measurements comparing DNA molecules with $n=0$ and $n=2$ mismatches show no apparent differences (Fig. (a,b)). In contrast, analysis of the densely sampled data in Fig. 2 shows the defect causes buckling to occur at a lower critical torque, $\Gamma_{\rm c}(n)$, which, ignoring small differences in torsional elasticity, is assumed to linearly increase from $\Gamma^*$ in proportion to ${({\rm Lk}_n^\dagger-{\rm Lk}_n^*)}$ (Fig. 1, Fig. 2d). Compensating trends account for differences in $\Gamma_{\rm c}(n=0)$ and $\Gamma_{\rm c}(n=2)$ even though ${\rm Lk}_0^\dagger$ and ${\rm Lk}_2^\dagger$ are similar: As (${\rm Lk}_n^\dagger-{\rm Lk}_n^*$) is reduced by *bending* compliance, *torsional* compliance of the defect causes ${\rm Lk}_n^*$ to increase. The enhanced torsional compliance is noticeable for a large defect (Fig. (c,d)). These trends are repeatable across different constant forces (tabulated in Fig. e).


abstract: We study the minimal flavor violation in the context of string effective field theory. Stringy selection rules indicate that $n$-point couplings among fermionic zero-modes and lightest scalar modes in the string effective action are given by a product of Yukawa couplings which are regarded as spurion fields of stringy and geometrical symmetries. Hence, Yukawa couplings determine the dynamics of flavor and CP violations. This observation strongly supports the hypothesis of minimal flavor violation in the Standard Model effective field theory.
author: Tatsuo Kobayashi${}^{1}$; Hajime Otsuka${}^{2}$
title: On Stringy Origin of Minimal Flavor Violation

EPHOU-21-011, KEK-TH-2340

# Introduction

Precision measurements of flavor and CP violations give us not only a precision test of the Standard Model (SM) but also clues of new physics beyond the SM. Current experimental searches put a strong constraint on the scale of new physics, namely the cutoff scale larger than the TeV scale in the context of low-energy effective field theory (EFT).

The SMEFT approach is a powerful tool to capture signatures of new physics beyond the SM, and the flavor- and CP-violating processes are sensitive to higher-dimensional operators justified below the cutoff scale. Recalling the observation that the SM has an $U(3)^5$ flavor symmetry in the limit of vanishing Yukawa couplings, it is natural to assume that the SM Yukawa couplings are the only source of flavor and CP violations. In this hypothesis called Minimal Flavor Violation (MFV), Yukawa couplings are considered spurion fields under, e.g., the $U(3)^5$ maximal flavor symmetry ; the higher-order couplings in the SMEFT are controlled by a product of Yukawa couplings to be invariant under the SM flavor and gauge symmetries. From the viewpoint of high-energy dynamics, the theory beyond the SM has to possess the non-trivial structure such that higher-dimensional operators generated by integrating out heavy states around the cutoff scale satisfy the criterion of the MFV hypothesis.

In the higher-dimensional theory, such as the string theory, only a limited class of $n$-point couplings is allowed in the four-dimensional (4D) EFT due to the existence of higher-dimensional gauge and Lorentz symmetries, and stringy selection rules. Explicit forms of $n$-point couplings were calculated in the string EFT on toroidal backgrounds, complex projective spaces, and some classes of Calabi-Yau manifolds with background gauge fields by utilizing conformal field theory (CFT) techniques and field theoretical methods. It turned out that higher-order couplings of matter zero-modes were constructed by the product of Yukawa couplings. For earlier attempts in the calculation of higher-order couplings, see Ref. for Type IIA intersecting D-brane models, Ref. for Type IIB magnetized D-brane models, Ref. for heterotic orbifold models, and Ref. for ${\cal N}=2$ topological string theory. Furthermore, the zero-mode solutions constitute a complete orthonormal system on general curved manifolds with background gauge fields; the $n$-point couplings of fermionic zero-modes and/or lightest scalar modes are determined by their $(n-1)$-point couplings in a similar way .

The phenomenological aspects of these higher-order couplings in the string EFT have not been fully investigated. The purpose of this paper is to explore the MFV hypothesis from the ultraviolet point of view by utilizing the structure of $n$-point couplings calculated in the string EFT. We find that various classes of string EFTs satisfy the criterion of MFV hypothesis, and flavor- and CP-violating processes are controlled by stringy selection rules.

# Higher-order couplings in string EFT

We investigate the basic structure of $n$-point couplings in the string EFT framework from the MFV point of view.

## $n$-point couplings

Let us begin with a class of 4D ${\cal{N}}=1$ EFTs arising from superstring theory compactified on toroidal as well as Calabi-Yau backgrounds.

On 4D ${\cal{N}}=1$ supersymmetric backgrounds, it has been known that the 4D ${\cal{N}}=1$ superpotential of matter zero-modes are expressed by $$\begin{aligned}
\begin{split}
W = y_{i_1 i_2 i_3}\Phi_{i_1} \Phi_{i_2} \Phi_{i_3} + \sum_{n=1} \frac{y_{i_1i_2\cdots i_{n+3}}}{\Lambda^n}\Phi_{i_1}\cdots \Phi_{i_{n+3}},
\label{eq:W}
\end{split}
\end{aligned}$$ where $\Phi_{i_n}$ denote matter chiral superfields originating from ten- and higher-dimensional gauge boson/gauginos for heterotic string theory and Type II string theory with D-branes, respectively. Here, $\Lambda$ represents the cutoff scale depending on the compactification scale.

The higher-order couplings are explicitly calculated in toroidal orbifolds in the CFT as well as the 4D EFT. It is remarkable that $n$-point couplings of matter zero-modes take the form $$\begin{aligned}
y_{i_1i_2\cdots i_{n}} = \sum_m y_{i_1i_2\cdots i_{n-2}m} \cdot y_{\bar{m}i_{n-1}i_{n}},
\label{eq:formula}
\end{aligned}$$ for intersecting D6-brane models , the T-dual magnetized D7-brane models , and heterotic orbifolds . Here, we sum over the gauge-invariant modes $\Phi_m$ which do not necessarily correspond to zero-modes but include heavy states. This expression is inherited from operator product expansions in the two-dimensional CFT: $$\begin{aligned}
    V_i (z) V_j (0) \sim \sum_k \frac{y_{ijk}}{z^{h_k -h_i -h_j}}V_k(0),
\end{aligned}$$ where $h_i$ denote the conformal dimensions of vertex operators $V_i$, corresponding to massless modes. Here, we rely on the perturbative calculation for $n$-point couplings, but non-perturbative effects such as D-brane instanton effects in Type IIA string theory lead to a similar structure of $n$-point couplings protected by non-anomalous symmetries , where the coefficient includes the instanton action factor $e^{-S_{\rm instanton}}$.

The explicit forms of $n$-point couplings are known on toroidal backgrounds, but it is challenging to derive $n$-point couplings on more general backgrounds such as Calabi-Yau threefolds. One of the difficulties in the calculation is due to the lack of analytical metric of entire Calabi-Yau threefolds and their submanifolds. However, in the heterotic string theory with standard embedding, $n$-point couplings of matter zero-modes can be extracted from the corresponding moduli couplings due to the identification between $SU(3)$ gauge bundle of heterotic $E_8$ gauge group and the tangent bundle of Calabi-Yau threefolds. It results in the same formula of $n$-point couplings (), where the $n$-point moduli couplings reduce to the derivative of Yukawa couplings with respect to the moduli fields . Furthermore, the equations of motion for matter zero-modes on general curved manifolds with background gauge fields also have a structure that $n$-point couplings among fermionic zero-modes and lightest scalar modes obey the form () by utilizing the complete orthonormal system of Dirac and Klein-Gordon operators. The orthogonality of matter zero-modes $\Phi_i$ indicates that the lightest eigenstates of operators $D_i = \nabla_i -i A_i$ with background gauge fields $A_i$ satisfy the relation: $$\begin{aligned}
\Phi_i \cdot \Phi_j \sim \sum_k y_{ijk}\cdot \Phi_k,
\end{aligned}$$ which results in an underlying structure of the formula () in a way similar to the operator product expansions in the CFT.

Note that the selection rule of $n$-point couplings () will be applicable to the non-supersymmetric (SUSY) background. It means that the non-holomorphic $n$-point couplings of massless modes also obey the formula () as confirmed in e.g., Ref. . For instance, four-fermion operators $y_{\bar{i}j\bar{k}l}\bar{\psi}_i\psi_j\bar{\psi}_k\bar{\psi}_l$ are generated by integrating out heavy states such as string modes as well as superpartners of SM fields.

In this way, the selection rule of $n$-point couplings () was explicitly demonstrated on toroidal and some classes of Calabi-Yau threefolds with background gauge fields and was also supported at the level of equations of motion . The structure of the formula () will be understood from the viewpoint of string scattering amplitudes. For the tree-level string amplitude, the vertex operators inserted at the position of $i_1,...,i_{n-2}$ on Riemann surfaces correspond to massless modes, which join to give a massless or massive mode ($\Phi_m$) propagating in the Riemann surface and emit two massless modes $\Phi_{i_{n-1}}$ and $\Phi_{i_n}$. These observations strongly suggest that the decomposition of $n$-point couplings into Yukawa couplings holds on general grounds if they admit the CFT description.

## Minimal flavor violation

We will consider the phenomenological interpretation of Eq. () predicted in the string EFT. More concrete estimations of $n$-point couplings in the low-energy EFT will be presented in the next section.

In the MFV hypothesis, Yukawa couplings are supposed to be spurion fields under the SM flavor symmetries, e.g., $U(3)^5$, and the SMEFT is constructed from the SM fields and Yukawa couplings, keeping the flavor and/or CP symmetries of quarks and leptons. The flavor-violating processes are then determined by the structure of Yukawa couplings. Remarkably, this treatment is realized in a certain class of string EFTs, since all the $n$-point couplings are decomposed into the product of Yukawa couplings, thereby yielding the magnitude of all the flavor- and CP-violating processes. Furthermore, $n$-point couplings are constrained by stringy selection rules and geometrical symmetries of the internal manifold in addition to the SM gauge symmetries, which allows us to restrict the unknown flavor symmetry of Yukawa couplings in the SMEFT prescription.

In typical string models, matter zero-modes have representations under continuous and discrete symmetries associated with stringy selection rules and geometrical symmetries of the internal manifold, and the corresponding Yukawa couplings also have some representations under them. In the string EFT, Yukawa couplings, as well as other couplings, are functions of moduli. That corresponds to the important aspect of MFV; Yukawa couplings are spurion fields. Since these moduli transform non-trivially under stringy and geometrical symmetries, couplings in the EFT and matter modes transform non-trivially.

One example is the modular symmetry $SL(2,\mathbb{Z})$ on toroidal (orbifold) backgrounds. Matter fields and their couplings transform non-trivially under the modular symmetry, but their Lagrangian is invariant. In particular, zero-mode wavefunctions of matter fields are described by the Jacobi theta function with modular weight 1/2 for D-brane models and 1 for heterotic orbifold models. Since $n$-point couplings are determined by the overlap integrals of $n$ number of matter zero-modes, these couplings are also described by modular forms as a function of the complex structure of the torus, which results in, e.g., $A^4$ or $S^4$ non-Abelian discrete flavor symmetries. (For more details, see Ref.  for heterotic orbifold models, and Ref.  for magnetized D-brane models.) In multi-moduli theories, the modular symmetry is enlarged to the symplectic modular symmetry $Sp(2h,\mathbb{Z})$, which governs the flavor and CP structures. For instance, in the heterotic string theory with standard embedding on Calabi-Yau threefolds, Yukawa and higher-order couplings computed by taking an appropriate number of derivatives of the prepotential are tensor representations of a subgroup of $Sp(2h,\mathbb{Z})$ in the Calabi-Yau moduli space (such as $S_4$ flavor symmetry ). These couplings depend on the complex structure and Kähler moduli for the fundamental and anti-fundamental representations of $E_6$ gauge group, respectively.

Another example is given by heterotic orbifold models, where stringy selection rules and geometrical symmetries lead to flavor symmetries such as $D_4$ and $\Delta(54)$ . Twisted modes are non-trivial representations of these symmetries, while untwisted modes are trivial singlet. Yukawa couplings, as well as other couplings, are also non-trivial representations when twisted modes develop their vacuum expectation values corresponding to the blow-up of orbifold singularities.

Starting from the supersymmetric constructions, the SUSY-breaking terms will cause the flavor violating processes in the low-energy EFT. In the typical string models, the rank of Yukawa couplings is 1 in the perturbative regime, i.e., the large volume and weak string coupling regime . It corresponds to the $U(2)$ flavor symmetry acting on first-two light generations. In this way, the flavor violations induced by the soft SUSY-breaking terms are approximately suppressed. Also, we can study more generic models, where Yukawa matrices correspond to rank-3 matrices. The soft scalar masses $m_{\Phi_{i}}^2$ and the $A$-terms induced by the SUSY-breaking fields $X$ are written by , $$\begin{aligned}
    m_{\Phi_{i}}^2 &= m_{3/2}^2 - \sum_X |F^X|^2 \partial_X\partial_{\bar{X}}\ln \partial_{\Phi_{i}} \partial_{\bar{\Phi}_i}K_{\rm matter},
    \nonumber\\
    A_{ijk} &= A_{i} + A_{j} +A_{k} - \sum_X \frac{F^X}{y_{ijk}} \partial_X (y_{ijk}),
\end{aligned}$$ with $$\begin{aligned}
    A_i = \sum_X F^X \partial_X \ln \left(e^{-K_{\rm mod}} \partial_{\Phi_{i}} \partial_{\bar{\Phi}_i}K_{\rm matter}\right),
\end{aligned}$$ where $m_{3/2}$ is the gravitino mass, $K_{\rm mod}$ and $K_{\rm matter}$ denote the moduli and matter Kähler potentials, respectively. Here, we assume that the matter Kähler metric is a diagonal form. When $X$ are moduli fields, soft terms are also invariant under geometrical symmetries of the internal manifold, such as the modular symmetries, since moduli fields and their $F$-terms are the same representations under the modular symmetries.

In this way, various classes of string EFTs satisfy the requirement of MFV hypothesis, and the flavor structure of higher-order couplings is determined by Yukawa couplings reflecting geometrical symmetries of the internal manifold. The representations of higher-order couplings under the geometrical symmetries are given by the tensor product of Yukawa couplings in a way similar to the MFV. Furthermore, the magnitudes of CP violation in Yukawa and higher-order couplings have a common origin in string axions associated with the complex structure or Kähler moduli of toroidal and Calabi-Yau geometries.

Even if string EFTs satisfy the MFV at the compactification scale, various physical stages may occur between the compactification scale and low energy one, (i) some modes become massive, and (ii) some scalar fields develop their vacuum expectation values. At the stage (i), we integrate out massive modes $\Phi_m$, and then new operators appear after the integration keeping the rule (). At the stage (ii), new couplings are induced by vacuum expectation values of scalar fields. For example, suppose that there is the coupling, $y_{ijk\ell}\Phi_i \Phi_j \Phi_k \Phi_\ell$ and $\Phi_i$ develops its vacuum expectation value. Then, a new coupling appears as $y'_{jk\ell} \Phi_j \Phi_k \Phi_\ell$ with $y'_{jk \ell}=y_{ijk\ell}\langle \Phi_i \rangle$. Since the coupling constant $y_{ijk\ell}$ is the spurion depending on moduli, the new coupling constant $y'_{jk \ell}$ is also the spurion as functions of $\langle \Phi_i \rangle$ as well as moduli. The transformation behavior of $y'_{jk \ell}$ under (flavor and CP) symmetries is inherited from one of $y_{ijk\ell}\langle \Phi_i \rangle$. Thus, the MFV framework is not violated and remains in the low-energy EFT, although the cutoff scale may shift from the compactification scale to the scale including $\langle \Phi_i \rangle$ as well as mass scale of massive modes. This statement would hold true for some fermion condensations. Thus, the MFV framework can be realized in various classes of string EFTs, and remain in the low-energy EFT.

# The structure of higher-order couplings in the low-energy EFT

We discuss the phenomenological implications of the formula () into higher-dimensional operators relevant to flavor physics in the low-energy EFT.

## Dimension-5 operators

Here, we begin with a dimension-5 operator of the 4D ${\cal{N}}=1$ effective action realizing the Minimal Supersymmetric Standard Model (MSSM) at the low-energy scale. The dimension-5 operator invariant under $SU(2)_L\times U(1)_Y$ is the Weinberg operator: $$\begin{aligned}
 \frac{y_{L_i HL_j H}}{\Lambda}(L_i H)(L_j H).
\end{aligned}$$ The formula () indicates that the four-point couplings $y_{L_i HL_j H}$ are determined by the Yukawa couplings including the lepton sector: $y_{L_i HL_j H} \simeq \sum_m y_{L_i \Phi_m H}\times y_{L_j \bar{\Phi}_m H}$, independently on whether right-handed neutrinos appear below the compactification scale or not.

The SM or its supersymmetric extension can be engineered using several stacks of D-branes, e.g., intersecting D-brane models. In the context of heterotic string theory, the model building of grand unified theories (GUTs) has been based on $SU(5)$, $SO(10)$, and $E_6$ GUTs at the intermediate scale, which is subsequently broken to the SM gauge group by a Wilson line. The selection rules of $n$-point couplings constrain the higher-dimensional operators not only in the SM but also in GUTs. For instance, the dimension-5 operators for fundamental representations in the SUSY $E_6$ GUT $$\begin{aligned}
\frac{a_{ijkl}}{\Lambda}27_i 27_j 27_k 27_l
\end{aligned}$$ are given by $a_{ijkl}=\partial_i \partial_j \partial_k \partial_l {\cal F}$ with ${\cal F}$ being the prepotential of Calabi-Yau threefolds, where $a_{ijkl}$ are generated by the contribution of massive modes in the formula (). In the large complex structure regime, the prepotential takes the form: $$\begin{aligned}
{\cal F} = {\cal F}_{\rm poly} +{\cal F}_{\rm inst},
\end{aligned}$$ where ${\cal F}_{\rm poly}$ denotes the cubic polynomial functions with respect to the complex structure moduli, and ${\cal F}_{\rm inst}$ represents quantum corrections corresponding to worldsheet instanton effects in the topological A-model. Since higher-dimensional operators have an origin in the quantum corrections, these exponentially suppressed corrections will not be significant for the flavor- and CP-violating processes.

The couplings of dimension-5 proton decay operators in the SUSY $SU(5)$ GUT $$\begin{aligned}
    \frac{b_{ijk\bar{l}}}{\Lambda}10_i 10_j 10_k \bar{5}_l
\end{aligned}$$ are also written by the product of Yukawa couplings: $$\begin{aligned}
    b_{ijk\bar{l}} = \sum_m y_{10_i10_j \Phi_m} \times y_{10_k \bar{5}_l \bar{\Phi}_m},
\end{aligned}$$ where $\Phi_m$ will include GUT Higgs fields or other heavy states. If these operators are prohibited by the existence of hidden $U(1)$ symmetries, $SU(5)$ singlet fields appear in the proton decay operators in a gauge-invariant way. Even in this case, the couplings of proton decay operators would be restricted to satisfy the relation (). On the toroidal backgrounds, the four-point couplings are given by the Jacobi theta function with respect to the complex structure of the torus $\tau$. Since such an elliptic function behaves as $e^{2\pi i \tau}$, these higher-order couplings can also be suppressed in the large complex structure regime ${\rm Im}(\tau)\gg 1$.

## Dimension-6 flavor-violating operators

Let us study the dimension-6 flavor-violating operators in the SMEFT following Ref. , including the 1350 CP-even and 1149 CP-odd operators in the absence of any flavor symmetry . Concerning the recent flavor anomalies such as $B$ physics, we focus on four-fermion operators: the $\Delta F=2$ quark operators $$\begin{aligned}
    \frac{c_{\bar{i}j\bar{k}l}}{\Lambda^2}(\bar{Q}_i \gamma_\mu P_{L,R} Q_j)(\bar{Q}_k \gamma^\mu P_{L,R} Q_l),
\end{aligned}$$ and the $\Delta F=1$ semileptonic operators $$\begin{aligned}
    &\frac{d_{\bar{i}j\bar{k}l}}{\Lambda^2}(\bar{Q}_i \gamma_\mu P_{L,R} Q_j)(\bar{L}_k \gamma^\mu P_{L,R} L_l),
\end{aligned}$$ with the 4D Gamma matrix $\gamma_\mu$, and the projective operators $P_{L,R}$, where we pick up specific four-fermion operators for illustrative purposes. Given the superpotential (), the four-fermion operators are generated by the spontaneous SUSY breaking or the contribution of stringy states, and couplings $\{c_{\bar{i}j\bar{k}l}, d_{\bar{i}j\bar{k}l}\}$ are decomposed into the product of moduli-dependent three-point couplings: $\sum_m y_{\bar{i}jm}y_{\bar{m}\bar{k}l}$ where the virtual mode $\phi_m$ corresponds to heavy superpartners or stringy modes appearing in the gauge-invariant way. (For the derivation of dimension-6 operators in the case of matching the MSSM onto the SMEFT, e.g., Ref. .) Here, the cutoff scale determined by the SUSY-breaking scale or the string scale should be larger than at least a few TeV scale due to recent flavor experiments, e.g., various meson mixing in $\Delta F=2$ processes and the lepton-flavor universality violation in the semileptonic decay processes of the $B$-meson . (For more details about the constraints on dimension-6 SMEFT operators, see, e.g., Ref. and references therein.) Furthermore, the rank 1 Yukawa couplings are typically realized at the leading order of string EFTs , corresponding to the $U(2)$ flavor symmetry acting on first-two light generations. This setup provides the useful EFT approach for the recent flavor anomalies .

When the four-fermion operators are generated by heavy superpartner contributions, the Yukawa couplings of quarks and leptons are constrained by a realization of their mass hierarchies and suppression of dangerous flavor and CP violations. Such a hierarchical structure would be possible at a specific point of moduli space such as the large complex structure regime, and the vacuum expectation values of moduli and axion fields appearing in Yukawa couplings will be fixed by the stabilization mechanism of the internal manifold such as flux compactifications with or without non-perturbative effects developed in Type IIB string theory .

## Derivative couplings

In addition to the SM flavor and gauge symmetries, the gauge and higher-dimensional Lorentz symmetries constrain the form of $n$-point couplings with derivatives. The derivative couplings in the higher-dimensional theory such as $\psi \bar{\psi}|D_M\phi|^2$ lead to the 4D couplings through the compactification: $\chi \bar{\chi}|D_\mu \xi|^2$, where $\chi$ and $\xi$ are 4D fermions and scalars with the 4D coordinates $x^\mu$, originating from higher-dimensional fermions $\psi$ and scalars $\phi$ with the coordinate in higher-dimensional spacetime $x^M$, respectively. Similarly, the coefficient of $\chi \bar{\chi}|D_\mu \xi|^2$ is the same form of $\chi \bar{\chi}|\xi|^2$ thorough the compactification of the internal space. When we identify $\xi$ with the Higgs fields, the coupling of $\chi \bar{\chi}|D_\mu \xi|^2$ operator corresponds to the four-point coupling of fermionic zero-modes and Higgs fields. In this way, the derivative couplings are also governed by $n$-point couplings and satisfy the criterion of the MFV.

Note that these analyses would be applicable to other derivative couplings with or without multiple scalar fields such as multi-Higgs models, which are often predicted in the string effective action. The effective action, as well as the flavor structure of the multiple Higgs fields, will also be constrained by the formula () (see, e.g., Ref. ).

# Conclusions

It has been known that the MFV hypothesis is a powerful tool to capture signatures of flavor and CP violations, but its high-energy behavior has not been established. In this paper, we focused on the properties of higher-order couplings in a certain class of string EFTs. We have shown that the stringy structure naturally satisfies the criterion of MFV hypothesis; all the flavor- and CP-violating processes are controlled by the product of Yukawa couplings. Furthermore, moduli fields appearing in Yukawa couplings transform non-trivially under stringy and geometrical symmetries, implying that Yukawa couplings are regarded as spurion fields under them. The MFV scenario also works in the low-energy EFT below the compactification scale.

In the context of string EFT, the appearance of flavor symmetries depends on the structure of Yukawa couplings obtained through the compactification of extra-dimensional spaces. One of the interesting symmetries is the modular symmetry. For the toroidal backgrounds, the modular symmetry $SL(2,\mathbb{Z})$ and its finite subgroups such as $A_4$ and $S_4$ play an important role in the flavor structure of quarks and leptons, where the Yukawa couplings have non-trivial representations under the modular symmetry. A similar phenomenon appears in the heterotic string theory with standard embedding on Calabi-Yau threefolds. Other interesting symmetries can be provided by stringy selection rules and geometrical symmetries of compact spaces. It would be fascinating to examine in detail the flavor and CP violations induced by higher-order couplings such as the semileptonic decays of $B$-mesons on concrete models and exemplify the realization of MFV hypothesis on other corners of string models.

Furthermore, we have only focused on the classical $n$-point couplings of matter fields. It would be interesting to incorporate quantum corrections into the matter couplings and check the validity of the formula () against these corrections.